Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 469,883 A→C 5.7% N338T (AAT→ACT)  mdlA → putative multidrug ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913469,8830AC 5.7% 194.7 / 2.5 70N338T (AAT→ACT) mdlAputative multidrug ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base A (28/38);  new base C (2/2);  total (30/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.46e-01

TGGCCAATGCTGGCGCTGGCATGGATGTTTAACATTGTGGAACGTGGTAGTGCTGCGTACAGCCGTATTCGCGCGATGCTGGCGGAAGCGCCGGTGGTGAACGATGGTAGCGAACCCGTGCCGGAAGGGCGTGGCGAACTGGATGTAAATATTCACCAGTTCACGTATCCGCAGACTGACCATCCTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGAC  >  NC_000913/469735‑470031
                                                                                                                                                    |                                                                                                                                                    
gggCCAATGCTGGCGCTGGCATGGATGTTTAACATTGTGGAACGTGGTAGTGCTGCGTACAGCCGTATTCGCGCGATGCTGGCGGAAGCGCCGGTGGTGAACGATGGTAGCGAACCCGTGCCGGAAGGGCGTGGCGAACTGGATGTaaa                                                                                                                                                      <  2:23084/148‑1 (MQ=255)
 ggCCAATGCTGGCGCTGGCATGGATGTTTAACATTGTGGAACGTGGTAGTGCTGCGTACAGCCGTATTCGCGCGATGCTGGCGGAAGCGCCGGTGGTGAACGATGGTAGCGAACCCGTGCCGGAAGGGCGTGGCGAACTGGATGTAAat                                                                                                                                                     <  1:119811/149‑1 (MQ=255)
   ccAATGCTGGCGCTGGCATGGATGTTTAACATTGTGGAACGTGGTAGTGCTGCGTACAGCCGTATTCGCGCGATGCTGGCGGAAGCGCCGGTGGTGAACGATGGTAGCGAACCCGTGCCGGAAGGGCGTGGCGAACTGGATGTAAatat                                                                                                                                                   >  2:1087545/1‑149 (MQ=255)
    cAATGCTGGCGCTGGCATGGATGTTTAACATTGTGGAACGTGGTAGTGCTGCGTACAGCCGTATTCGCGCGATGCTGGCGGAAGCGCCGGTGGTGAACGATGGTAGCGAACCCGTGCCGGGAGGGCGTGGCGAACTGGATGTAAATAtt                                                                                                                                                  >  2:227033/1‑149 (MQ=255)
     aaggctggcgctgggATGGGTGTTTAACATTGGGGAACGGGGTAGTGCTGCGGACAGCCGTATTCGCGCGATGCTGGCGGAAGCGCCGGTGGTGAACGATGGTAGCGAACCCGTGCCGGAAGGGCGTGGCGAACTGGATGTAAATATTc                                                                                                                                                 <  1:257663/146‑1 (MQ=255)
      aTGCTGGCGCTGGCATGGATGTTTAACATTGTGGAACGTGGTAGTGCTGCGTACAGCCGTATTCGCGCGATGCTGGCGGAAGCGCCGGTGGTGAACGATGGTAGCGAACCCGTGCCGGTAGGGCGCGGCGAACTGGATGTAAATATTCa                                                                                                                                                <  1:709375/149‑1 (MQ=255)
       tGCTGGCGCTGGCATGGATATTTAACATTGTGGCACGTGGTAGTGCTGCGTATAGCCGTATTCGCGCGATGCTGGCGGAATCGCCGGTGGTGAACGATGGTAGCGAACCCGTGCCGGAAGGGCGTGGCGAACTGGATGTAAATATTCAc                                                                                                                                               <  1:1000919/149‑1 (MQ=255)
                tggcATGGATGTTTAACATTGTGGAACGTGGTAGTGCTGCGTACAGCCGTATTCGCGCGATGCTGGCGGAAGCGCCGGTGGTGAACGATGGTAGCGAACCCGTGCCGGAAGGGCGTGGCGAACTGGATGTAAATATTCACCAGTTCACg                                                                                                                                      >  2:764655/1‑149 (MQ=255)
                tggcATGGATGTTTAACATTGTGGAACGTGGTAGTGCTGCGTACAGCCGTATTCGCGCGATGCTGGCGGAAGCGCCGGTGGTGAACGATGGTAGCGAACCCGTGCCGGAAGGGCGTGGCGAACTGGATGTAAATATTCACCAGTTCACg                                                                                                                                      >  2:764695/1‑149 (MQ=255)
                   cATGGATGTTTAACATTGTGGAACGTGGTAGTGCTGCGTACAGCCGTATTCGCGCGATGCTGGCGGAAGCGCCGGTGGTGAACGATGGTAGCGAACCCGTGCCGGAAGGGCGTGGCGAACTGGATGTAAATATTCACCAGTTCACGTAt                                                                                                                                   <  2:587343/149‑1 (MQ=255)
                    aTGGATGTTTAACATTGTGGAACGTGGTAGTGCTGCGTACAGCCGTATTCGCGCGATGCTGGCGGAAGCGCCGGTGGTGAACGATGGTAGCGAACCCGTGCCGGAAGGGCGTGGCGGACTGGATGTAAATATTCACCAGTTCACGTagc                                                                                                                                  >  2:1276474/1‑147 (MQ=255)
                        ttGTTTAGCATTGTGGAACGGGGTGGGGCTGGCGACAGGCGGGTTGGCGGGATGCGGGCGGAAGCGGGGGGGGTGAACGATGGTAGCGAACCTGTGCCGGAAGGGCGTGGCGCACTGGATGTAAATATTCACCTGTTCACGTATCCGcc                                                                                                                              <  2:128074/148‑2 (MQ=255)
                           tttAACATTGTGGAACGTGGTAGTGCTGCGTACAGCCGTATTCGCGCGATGCTGGCGGAAGCGCCGGTGGTGAACGATGGTAGCGAACCCGTGCCGGAAGGGCGTGGCGAACTGGAGGTAAATATTCACCAGGTCACGTATCCGCAGAc                                                                                                                           >  1:754098/1‑149 (MQ=255)
                                 aTTGTGGGACGTGGTAGTGCTGCGTACGGCCGTATTCGCGCGATGCTGGCGGAAGCGCCGGTGGTGAACGATGGTAGCGAACCCGTGCCGGAAGGGCGTGGCGAACTGGATGTAAATATTCACCAGTTCACGTATCCGCAGACTGACCa                                                                                                                     <  2:68605/149‑1 (MQ=255)
                                  ttGTGGAACGTGGTAGTGCTGCGTACAGCCGTATCGGCGCGGTGCAGGCGGAGACGCAGGAGGTGAACGGTGGGATCGCACCCGAGCCGGAAGGGACTGCCGCGCAGGAGAGAAAAAATGACCAGTTCAAGTATCCGCAGACTGAcgac                                                                                                                    >  2:425576/1‑146 (MQ=255)
                                     tgGAACGTGGTAGTGCTGCGTACAGCCGTATTCGCGCGATGCTGGCGGAAGCGCCGGTGGTGAACGATGGTAGCGAACCCGTGCCGGAAGGGCGTGGCGAACTGGATGTAAATATTCACCAGTTCACGTATCCGCAGACTGACCATCCt                                                                                                                 <  2:1155171/149‑1 (MQ=255)
                                       gAACGTGGTAGTGCTGCGTACAGCCGTATTCGCGCGATGCTGGCGGAAGCGCCGGTGGTGAACGATGGTAGCGAACCCGTGCCGGAAGGGCGTGGCGAACTGGATGTAAATATTCACCAGTTCACGTATCCGCAGACTGACCATCCTgc                                                                                                               >  1:495813/1‑149 (MQ=255)
                                        aaCGTGGTAGTGCTGCGTACAGCCGTATTCGCGCGATGCTGGCGGAAGCGCCGGTGGTGAACGATGGTAGCGCACCCGTGCCGGAAGGGCGTGGCGAGCTGGATGTAAATATTCACCAGTTCACGTATCCGCAGACTGACCGTCCTgag                                                                                                              >  2:372118/1‑147 (MQ=255)
                                          cGTGGTAGTGCTGCGTACAGCCGTATTCGCGCGATGCTGGCGGAAGCGCCGGTGGTGAACGATGGTAGCGAACCCGTGCCGGAAGGGCGTGGCGAACTGGATGTAAATATTCACCAGTTCACGTATCCGCAGACTGACCATCCTGCGCt                                                                                                            >  2:1128517/1‑149 (MQ=255)
                                            ttgTAGTGCTGCGTACAGCCGTATGCGCGCGATGCTGGCGGAAGCGCCGGTGGTGAACGATGGTAGCGAACCCGTGCCGGAAGGGCGTGGCGAACTTGATGTAAATATTCACCAGTTCACGTATCCGCAGACTGACCATCCTGCGCTgg                                                                                                          <  2:1192550/147‑1 (MQ=255)
                                                aGTGCTGCGTACAGCCGTATTCGCGCGATGCTGGCGGAAGCGCCGGTGGTGAACGATGGTAGCGAACCCGTGCCGGAAGGGCGTGGCGAACTGGATGTAAATATTCACCAGTTCACGTATCCGCAGACTGACCATCCTGCGCTGGaaaa                                                                                                      <  1:492078/149‑1 (MQ=255)
                                                   gctgcGTACAGCCGTATTCGCGCGATGCTGGCGGAAGCGCCGGTGGTGAACGATGGTAGCGAACCCGTGCCGGAAGGGCGTGGCGAACTTGATGTAAATATTCACCAGTTCACGTATCCGCAGACTGACCATCCTGCGCTGGAAAACGt                                                                                                   <  2:999711/149‑1 (MQ=255)
                                                          aCAGCCGTATTCGCGCGATGCTGGCGGAAGCGCCGGTGGTGAACGATGGTAGCGAACCCGTGCCGGAAGGGCGTGGCGAACTGGATGTAAATATTCACCAGTTCACGTATCCGCAGACTGACCATCCTGCGCTGGAACACGTCAATTTc                                                                                            >  1:103271/1‑149 (MQ=255)
                                                           cAGCCGGATTCGCGCGATGCGGGCGGAAGCGCCTGTGGTGAACGATGGTAGCGAACCCGTGCCGGAAGGGCGTGGCGAACTGGATGTAAATATTCACCAGTTCACGTATCCGCAGACTGACCATCCTGCGCTGGAAAACGTCAATTTCg                                                                                           <  1:130114/149‑1 (MQ=255)
                                                            aGCCGTATTCGCGCGATGCTGGCGGAAGCGCCGGTGGTGAACGATGGTAGCGAACCCGTGCCGGAAGGGCGTGGCGAACTGGATGTAAATATTCACCAGTTCACGTATCCGCAGACTGACCATCCTGCGCTGGAAAACGTCAATTTCGc                                                                                          <  2:524426/149‑1 (MQ=255)
                                                             gCCGTATTCGCGCGATGCTGGCGGAAGCGCCGGTGGTGAACGATGGTAGCGAACCCGTGCCGGAAGGGCGTGGCGAACTGGATGTAAATATTCACCAGTTCACGTATCCGCAGACTGACCATCCTGCGCTGGAAAACGTCAATTTCGcc                                                                                         <  2:883482/149‑1 (MQ=255)
                                                                      gcgcgATGCTGGCGGAAGCGCCGGTGGTGAACGATGGTAGCGAACCCGTGCCGGAAGGGCGTGGCGAACTGGATGTAAATATTCACCAGTTCACGTATCCGCAGACTGACCATCCTGCGCTGGAAAACGTCAATTTCGCCCTGAAAccc                                                                                >  1:892766/1‑149 (MQ=255)
                                                                       cgcgATGCTGGCGGAAGCGCCGGTGGTGAACGATGGTAGCGAACCCGTGCCGGAAGGGCGTGGCGAACTGGATGTAAATATTCACCAGTTCACGTATCCGCAGACTGACCATCCTGCGCTGGAAAACGTCACTTTCGCCCCGAACCCCg                                                                               >  2:330470/1‑149 (MQ=255)
                                                                              cTGGCGGAAGCGGCGGTGGTGAACGATGGTAGCGAACCCGTGCCGGACGGGCGTGGCGAACTGGATGTAAATATTCACCAGTTCACGTATCCGCAGACTGACCATCCTGCGCTGGAAAACGTCAATTTCGCCCTGAACCCCGGTCAGAt                                                                        <  2:1311929/149‑1 (MQ=255)
                                                                                  gggAAGGGCCGGTGGTGACCGATGGTAGCGAACCCGTGCCGGAAGGGCGTGGCGAACTGGATGTAAATATTCACCAGTTCACGTATCCGCAGACTGACCATCCTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTg                                                                    <  2:1137055/148‑1 (MQ=255)
                                                                                   ggAAGCGCCGGTGGTGAACGATGGTAGCGAACCCGTGCCGGAAGGGCGTGGCGAACTGGATGTAAATATTCACCAGTTCACGTATCCGCAGACTCACCATCCTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTgg                                                                   <  1:691934/149‑1 (MQ=255)
                                                                                   ggAAGCGCCGGTGGTGAACGATGGTAGCGAACCCGGGCCGGAAGGGCGTGGCGAACTGGATGTAAATATTCACCAGTTCACGTATCCGCAGACTGACCATCCTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTgg                                                                   <  2:1030039/149‑1 (MQ=255)
                                                                                   ggAAGCGCCGGTGGGGAACGATGGTAGCGAACCCGTGCCGGAAGGGCGTGGCGAACTGGATGTAAATATTCACCAGTTCACGTATCCGCAGACTGACCATCCTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTgg                                                                   <  1:58878/149‑1 (MQ=255)
                                                                                     aaGCGCCGGTGGTGAACGATGGTAGCGAACCCGTGCCGGAAGGGCGTGGCGAACTGGATGTAAATATTCACCAGTTCACGTATCCGCAGACTGACCATCCTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGt                                                                 >  1:944004/1‑149 (MQ=255)
                                                                                       gccccGGCGGGGAACGAGGGCAGCGAACCCGTGCCGGAAGGGCGTGGCGAACTGGATGTAAATATTCACCAGTTCACGTATCCGCAGACTGACCATCCTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTAt                                                               <  1:437642/146‑1 (MQ=255)
                                                                                         gcCGGTGGTGAACGATGGTAGCGAACCCGTGCCGGAAGGGCGTGGCGAACTGGATGTAAATATTCACCAGTTCACGTATCCGCAGGCTGACCATCCTGCGCTGGAAACCGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCt                                                             >  1:1153409/1‑149 (MQ=255)
                                                                                           cGGTGGTGAACGATGGTAGCGAACCCGTGCCGGAAGGGCGTGGCGAACTGGATGTAAATATTCACCAGTTCACGTATCCGCAGACTGACCATCCTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGc                                                           <  1:772415/149‑1 (MQ=255)
                                                                                            ggtggtGAACGAAGGTAGAGAACACGTGCCGGAAGGGCGTGGCGAACTGGATGTAAATATTCACCAGTTCACGTATCCGCAGACTGACCATCCTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCGGCg                                                          >  2:1288514/1‑149 (MQ=255)
                                                                                             gtgcgaaaCGATGGTAGCGAACCCGTGCCGGAAGGGCGTGCCGAACTGGATGTAAAGATTCACCAGTTCACGTATCCGCAGACTGACCATCCTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCgg                                                         <  2:34788/143‑1 (MQ=255)
                                                                                                gtgAACGATGGTAGCGAACCCGTGCCGGAGGGGCGTGGCGAACTGGATGTAAATATTCACCAGTTCACGTATCCGCAGACTGACCATCCTGCGCTGGAAAACGACAAAGTAGCCCAGGAACCCGGTCAGATGCTGGGTAGCTGCGGGcc                                                      >  2:283858/1‑149 (MQ=255)
                                                                                                gtgAACGATGGTACGGAACCCGTGCCGGAAGGGCGTGGCGAACTGGATGTAAATATTCACCAGTTCACGTATCCGCAGACTGACCATCCTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGcc                                                      <  1:357507/149‑1 (MQ=255)
                                                                                                   aaCGATGGTAGCGAACCCGTGCCGGAAGGGCGTGGCGAACTGGATGTAAATATTCACCAGTTCACGTATCCGCAGACTGACCATCCTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGAc                                                   >  1:100344/1‑149 (MQ=255)
                                                                                                     cGATGGTAGCGAACCCGTGCCGGAAGGGCGTGGCGAACTGGATGTAAATATTCACCAGTTCACGTATCCGCAGACTGACCATCCTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTg                                                 >  1:137370/1‑149 (MQ=255)
                                                                                                      gATGGTAGCGAACCCGTGCCGGAAGGGCGTGGCGAACTGGATGTAAATATTCACCAGTCCACGTATCCGCAGACTGTCCACCCTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTgg                                                <  1:1276474/149‑1 (MQ=255)
                                                                                                       gtGGTAGCGAACCCGGGCCGGAAGGGCGTGGCGAACTGGATGTAAATATTCACCAGTTCACGTATCCGCAGACTGACCATCCTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGt                                               <  2:384445/148‑1 (MQ=255)
                                                                                                         gggaccGAACCCGTGCCGGAAGGGCGTGGCGAACTCGATGTAAATATTCACCAGTTCACGTATCCGCAGACTGACCATCCTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTc                                             <  1:749407/144‑1 (MQ=255)
                                                                                                         ggTAGCGAACCCGTGCCGGAAGGGCGTGGCGAACTGGATGTAAATATTCACCAGTTCACGTATCCGCAGACTGACCATCCTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTc                                             <  1:632081/149‑1 (MQ=255)
                                                                                                           tAGCGAACCCGTGCCGGAAGGGCGTGGCGAACTGGATGTAAATATTCACCAGTTCACGTATCCGCAGACTGACCATCCTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCg                                           <  1:764695/149‑1 (MQ=255)
                                                                                                           tAGCGAACCCGTGCCGGAAGGGCGTGGCGAACTGGATGTAAATATTCACCAGTTCACGTATCCGCAGACTGACCAGCCTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCg                                           <  1:764655/149‑1 (MQ=255)
                                                                                                            aGCGAACCCGTGCCGGAAGGGCGTGGCGAACTGGATGTAAATATTCACCAGTTCACGTATCCGCAGACTGACCATCCTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCAGCGGGCCGACTGGTTCCgg                                          >  2:549430/1‑149 (MQ=255)
                                                                                                              cGAACCCGTTCCGGAAGGGCGTGGCGAACTGTATGTAAATATTCACCAGTTCACGTATCCGCAGACTGACCATCCTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCa                                        <  1:1222432/149‑1 (MQ=255)
                                                                                                              cGAACCCGTGCCGGAAGGGCGTGGCGAACTGGATGTAAATATTCACCAGTTCACGTATCCGCAGACTGACCATCCTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCa                                        >  2:1352751/1‑149 (MQ=255)
                                                                                                                   ccGTGCCGGAAGGGCGTGGCGAACTGGATGTAAATATTCACCAGTTCACGTATCCGCAGACTGACCATCCTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGt                                   >  2:1159769/1‑149 (MQ=255)
                                                                                                                    cGTGCCGGAAGGGCGTGGCGAACTGGATGTAAATATTCACCAGTTCACGTATCCGCAGACTGACCATCCTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTa                                  >  2:431933/1‑149 (MQ=255)
                                                                                                                            aaGGGCGTGGCGAACTGGATGTAAATATTCACCAGTTCACGTATCCGCAGACTGACCATCCTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCtgttg                          <  2:495813/149‑1 (MQ=255)
                                                                                                                             aGGGCCTGGCTGACGGCATGTGAATAGTCACCAGTTCACGTAGCCGCAGACTGACCACTCTGAGCTGGAATCCGTCAATTGCGCCCTGACACCCGGTCAGATGCTGGGTATCTGCCGGCCGACTGGTTTCGGCAAAAGTACCCtgttgt                         <  2:1094909/149‑1 (MQ=255)
                                                                                                                              gggCGTGGCGAACTGGATGTAAATAATCACCAGTTCACGTATCCGCAGACTGACCATCCTGCGCTGGAAAACGTCAACTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTc                        <  1:180316/149‑1 (MQ=255)
                                                                                                                                 cGTGGCGAACTGGATGTAACTATTCACCAGTTCACGTATCCGCAGACTGACCATCCTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCtgttgt                         >  2:113681/1‑145 (MQ=255)
                                                                                                                                 cGTGGCGAACTGGATGTAACTATTCACCAGTTCACGTATCCGCAGACTGACCATCCTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCtgttgt                         <  1:113692/145‑1 (MQ=255)
                                                                                                                                 cGTGGCGAACTGGATGTAACTATTCACCAGTTCACGTATCCGCAGACTGACCATCCTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCtgttgt                         <  1:113681/145‑1 (MQ=255)
                                                                                                                                 cGTGGCGAACTGGATGTAACTATTCACCAGTTCACGTATCCGCAGACTGACCATCCTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCtgttgt                         >  2:113692/1‑145 (MQ=255)
                                                                                                                                 cGTGGCGAACTGGATGTAAATATTCACCAGTTCACGTATCCGCAGACTGACCATCCTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTAGCt                     >  1:839461/1‑149 (MQ=255)
                                                                                                                                  gcggCGCACTGGACGTAAATATTCACCAGTTCACGTATCCGCAGACTGACCACCCTGCGCTGGAACCCGTCAATTTCGCCCCGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTc                    <  2:1084510/147‑1 (MQ=255)
                                                                                                                                       gAACTGGATGTAAATATTCACCAGTTCACGTATCCGCAGACTGACCATCCTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCa               >  1:347371/1‑149 (MQ=255)
                                                                                                                                       gAACTGGATGTAAATATTCACCAGTTCACGTATCCGCAGACTGACCATCCTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCa               >  1:766272/1‑149 (MQ=255)
                                                                                                                                       gAACTGGATGGAAATATTCACCAGTTCACGTATCCGCAGACTGACCATCCTGCGCTGGAAAACGTCAATTGCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCa               <  2:759586/149‑1 (MQ=255)
                                                                                                                                         aCTGGATGTAAATATTCACCAGTTCACGTATCCGCAGACTGACCATCCTGCGCTGGAAAACGTCAATTTCGCCCTGAACCCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGc             <  1:1164440/149‑1 (MQ=255)
                                                                                                                                          cTGGATGTAAATATTCACCAGTTCACGTATCCGCAGACTGACCATCCTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCcgcg            >  2:91028/1‑149 (MQ=255)
                                                                                                                                          cTGGATGTAAATATTCACCAGTTCACGTATCCGCAGACTGACCATCCTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCg            <  1:751309/149‑1 (MQ=255)
                                                                                                                                            ggATGTAAATATTCACCAGTTCACGTATCCGCAGACTGACCATCCTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTc          <  1:1342514/149‑1 (MQ=255)
                                                                                                                                                gTAAATATTCACCAGTTCACGTATCCGCAGACTGACCATCCTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCAttt      >  1:252054/1‑149 (MQ=255)
                                                                                                                                                gTAAATATTCACCAGTTCACGTATCCGCAGACTGACCATCCTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCAttt      >  1:577609/1‑149 (MQ=255)
                                                                                                                                                 tAAATATTCACCAGTTCACGTATCCGCAGACTGACCATCCTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTc     <  2:152513/149‑1 (MQ=255)
                                                                                                                                                    atatTCACCAGTTCACGTATCCGCAGACTGACCATCCTGCGCTGGAAAACGTCAATTTCGCCCTGAACCCCGGCCAGATGCCGGGTATCTGCGGGCCGACTGGTTCCCGCAAAAGTACCCTGTTGCCGCTCATTGCGCGTCCTTTCgcc  >  2:334343/1‑147 (MQ=255)
                                                                                                                                                    atatTCACCAGTTCACGTATCCGCAGACTGACCATCCTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGCCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGAc  >  2:654126/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TGGCCAATGCTGGCGCTGGCATGGATGTTTAACATTGTGGAACGTGGTAGTGCTGCGTACAGCCGTATTCGCGCGATGCTGGCGGAAGCGCCGGTGGTGAACGATGGTAGCGAACCCGTGCCGGAAGGGCGTGGCGAACTGGATGTAAATATTCACCAGTTCACGTATCCGCAGACTGACCATCCTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGAC  >  NC_000913/469735‑470031

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: