Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,021,257 A→C 5.4% intergenic (‑80/+65) dadX ← / ← dadA‑II alanine racemase/D‑amino acid:quinone oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,021,2570AC 5.4% 948.8 / 20.6 308intergenic (‑80/+65)dadX/dadA‑IIalanine racemase/D‑amino acid:quinone oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (158/127);  new base C (9/11);  total (167/141)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 330 total aligned reads displayed.

TCAGTTCACGGGCCAGACGGTAGTTGTGGCGCAGGGCCTGGAGGTCGATCAGGGCGCGGGCGGGACGCATAGCGGCTGCCTTATGGTGGTTCGGGTTGAAATCGGTGGTGGATTATTCGCGGGCTTGCCCGCTCCCACAGGTCGTGCACAGCTTTCAAAAGCGGTGGTGTTCCTGTGGGAGCGGGCAAGCCCGCAAATGGGCCGGCACAGGCTGTTACTGCTGGTGAGCCGGCGCGCTTTGCCCGTGCTTGGCCACTTCCCGGCTATTGCCATAACGGGAGATATCCAGGCCTTCGGCA  >  NC_002947/6021108‑6021406
                                                                                                                                                     |                                                                                                                                                     
tCAGTTCACGGGCCAGACGGTAGTTGTGGCGCAGGGCCTGGAGGTCGATCAGGGCGCGGGCGGGGCGCATAGCGGCTGCCTTATGGTGGTTCGGGTTGAAATCGGTGGTGGATTATTCGCGGGCTTGCCCGCTCCCACAGGTCGTGcaca                                                                                                                                                       <  1:5224983/150‑1 (MQ=255)
tCAGTTCACGGGCCAGACGGTAGTTGTGGCGCAGGGCCTGGAGGTCGATCAGGGCGCGGGCGGGACGCATAGCGGCTGCCTTATGGTGGTTCGGGTTGAAATCGGTGGTGGATTATTCGCGGGCTTGCCCGCTCCCACAGGTCGTGcaca                                                                                                                                                       <  2:5618177/150‑1 (MQ=255)
tCAGTTCACGGGCCAGACGGTAGTTGTGGCGCAGGGCCTGGAGGTCGATCAGGGCGCGGGCGGGACGCATAGCGGCTGCCTTATGGTGGTTCGGGTTGAAATCGGTGGTGGATTATTCGCGGGCTTGCCCGCTCCCACAGGTCGTGcaca                                                                                                                                                       <  2:7443815/150‑1 (MQ=255)
tCAGTTCACGGGCCAGACGGTAGTTGTGGCGCAGGGCCTGGAGGTCGATCAGGGCGCGGGCGGGACGCATAGCGGCTGCCTTATGGTGGTTCGGGTTGAAATCGGTGGTGGATTATTCGCGGGCTTGCCCGCTCCCACAGGTCGTGcaca                                                                                                                                                       <  2:3546530/150‑1 (MQ=255)
   gtgcACGGGCCAGACGGTACAAGTGCCGCACGGCATGTAGGTCGATCAGGGCGCCGGCGGGACGCATAGCGCCTGCCTAATGGTGGTTCGGGATGAAATCGGTGGGGGATTATTCGCGGGCTTGCCCGCTCCCACAGGTCGTGCACAGCt                                                                                                                                                    >  2:8071941/4‑150 (MQ=255)
   gTTCACGGGCCAGACGGTAGTTGTGGCGCAGGGCCTGGAGGTCGATCAGGGCGCGGGCGGGACGCATAGCGGCTGCCTTATGGTGGTTCGGGTTGAAATCGGTGGTGGATTATTCGCGGGCTTGCCCGCTCCCACAGGTCGGGCACAGCt                                                                                                                                                    >  2:3185386/1‑150 (MQ=255)
     tCACGGGCCAGACGGTAGTTGTGGCGCAGGGCCTGGAGGTCGATCAGGGCGCGGGCGGGACGCATAGCGGCTGCCTTATGGTGGTTCGGGTTGAAATCGGTGGTGGATTATTCGCGGGCTTGCCCGCTCCCACAGGTCGTGCACAGCttt                                                                                                                                                  <  2:6967541/150‑1 (MQ=255)
        cGGGCCAGACGGTAGTTGTGGCGCAGGGCCTGGAGGTCGATCAGGGCGCGGGCGGGACGCATAGCGGCTGCCTTATGGTGGTTCGGGTTGAAATCGGTGGTGGATTATTCGCGGGCTTGCCCGCTCCCACAGGTCGTGCACAGCTTTCaa                                                                                                                                               <  2:1393908/150‑1 (MQ=255)
           gCCAGACGGTAGTTGTGGCGCAGGGCCTGGAGGTCGATCAGGGCGCGGGCGGGACGCATAGCGGCTGCCTTATGGTGGTTCGGGTTGAAATCGGTGGTGGATTATTCGCGGGCTTGCCCGCTCCCACAGGTCGTGCACAGCTTTCAAAAg                                                                                                                                            >  2:4847054/1‑150 (MQ=255)
            ccAGACGGTAGTTGTGGCGCAGGGCCTGGAGGTCGATCAGGGCGCGGGCGGGACGCATAGCGGCTGCCTTATGGTGGTTCGGGTTGAAATCGGTGGTGGATTATTCGCGGGCTTGCCCGCTCCCACAGGTCGTGCACAGCTTTCAAAAGc                                                                                                                                           >  2:2586158/1‑150 (MQ=255)
             cAGACGGTAGTTGTGGCGCAGGGCCTGGAGGTCGATCAGGGCGCGGGCGGGACGCATAGCGGCTGCCTTATGGTGGTTCGGGTTGAAATCGGTGGTGGATTATTCGCGGGCTTGCCCGCTCCCACAGGTCGTGCACAGCTTTCAAAAGCg                                                                                                                                          <  2:232650/150‑1 (MQ=255)
             cAGACGGTAGTTGTGGCGCAGGGCCTGGAGGTCGATCAGGGCGCGGGCGGGACGCATAGCGGCTGCCTTATGGTGGTTCGGGTTGAAATCGGTGGTGGATTATTCGCGGGCTTGCCCGCTCCCACAGGTCGTGCACAGCTTTCAAAAGCg                                                                                                                                          <  1:1037824/150‑1 (MQ=255)
              aGACGGTAGTTGTGGCGCAGGGCCTGGAGGTCGATCAGGGCGCGGGCGGGACGCATAGCGGCTGCCTTATGGTGGTTCGGGTTGAAATCGGTGGTGGATTATTCGCGGGCTTGCCCGCTCCCACAGGTCGTGCACAGCTTTCAAAAGCgg                                                                                                                                         <  1:6747602/150‑1 (MQ=255)
              aGACGGTAGTTGTGGCGCAGGGCCTGGAGGTCGATCAGGGCGCGGGCGGGACGCATAGCGGCTGCCTTATGGTGGTTCGGGTTGAAATCGGTGGTGGATTATTCGCGGGCTTGCCCGCTCCCACAGGTCGTGCACAGCTTTCAAAAGCgg                                                                                                                                         >  2:2444486/1‑150 (MQ=255)
               gACGGTAGTTGTGGCGCAGGGCCTGGAGGTCGATCAGGGCGCGGGCGGGACGCATAGCGGCTGCCTTATGGTGGTTCGGGTTGAAATCGGTGGTGGATTATTCGCGGGCTTGCCCGCTCCCACAGGTCGTGCACAGCTTTCAAAAGCggt                                                                                                                                        <  2:3147073/150‑1 (MQ=255)
               gACGGTAGTTGTGGCGCAGGGCCTGGAGGTCGATCAGGGCGCGGGCGGGACGCATAGCGGCTGCCTTATGGTGGTTCGGGTTGAAATCGGTGGTGGATTATTCGCGGGCTTGCCCGCTCCCACAGGTCGTGCACAGCTTTCAAAAGCggt                                                                                                                                        <  1:5715542/150‑1 (MQ=255)
               gACGGTAGTTGTGGCGCAGGGCCTGGAGGTCGATCAGGGCGCGGGCGGGACGCATAGCGGCTGCCTTATGGTGGTTCGGGTTGAAATCGGTGGTGGATTATTCGCGGGCTTGCCCGCTCCCACAGGTCGTGCACAGCTTTCAAAAGCggt                                                                                                                                        <  1:1145680/150‑1 (MQ=255)
                aCGGTAGTTGTGGCGCAGGGCCTGGAGGTCGATCAGGGCGCGGGCGGGACGCATAGCGGCTGCCTTATGGTGGTTCGGGTTGAAATCGGTGGTGGATTATTCGCGGGCTTGCCCGCTCCCACAGGTCGTGCACAGCTTTCAAAAGCggtg                                                                                                                                       >  2:6805690/1‑150 (MQ=255)
                aCGGTAGTTGTGGCGCAGGGCCTGGAGGTCGATCAGGGCGCGGGCGGGACGCATAGCGGCTGCCTTATGGTGGTTCGGGTTGAAATCGGTGGTGGATTATTCGCGGGCTTGCCCGCTCCCACAGGTCGTGCACAGCTTTCAAAAGCggtg                                                                                                                                       >  1:2620324/1‑150 (MQ=255)
                 cGGTAGTTGTGGCGCAGGGCCTGGAGGTCGATCAGGGCGCGGGCGGGACGCATAGCGGCTGCCTTATGGTGGTTCGGGTTGAAATCGGTGGTGGATTATTCGCGGGCTTGCCCGCTCCCACAGGTCGTGCACAGCTTTCAAAAGCggtgg                                                                                                                                      <  1:6598704/150‑1 (MQ=255)
                     aGTTGTGGCGCAGGGCCTGGAGGTCGATCAGGGCGCGGGCGGGACGCATAGCGGCTGCCTTATGGTGGTTCGGGTTGAAATCGGTGGTGGATTATTCGCGGGCTTGCCCGCTCCCACAGGTCGTGCACAGCTTTCAAAAGCGGTGGTGtt                                                                                                                                  >  1:4353822/1‑150 (MQ=255)
                         gtgGCGCAGGGCCTGGAGGTCGATCAGGGCGCGGGCGGGACGCATAGCGGCTGCCTTATGGTGGTTCGGGTTGAAATCGGTGGTGGATTATTCGCGGGCTTGCCCGCTCCCACAGGTCGTGCACAGCTTTCAAAAGCGGTGGTGTTcctg                                                                                                                              >  1:362721‑M2/1‑146 (MQ=255)
                          tgGCGCAGGGCCTGGAGGTCGATCAGGGCGCGGGCGGGACGCATAGCGGCTGCCTTATGGTGGTTCGGGTTGAAATCGGTGGTGGATTATTCGCGGGCTTGCCCGCTCCCACAGGTCGTGCACAGCTTTCAAAAGCGGTGGTGTTcctgt                                                                                                                             >  2:1828589‑M2/1‑145 (MQ=255)
                          tgGCGCAGGGCCTGGAGGTCGATCAGGGCGCGGGCGGGACGCATAGCGGCTGCCTTATGGTGGTTCGGGTTGAAATCGGTGGTGGATTATTCGCGGGCTTGCCCGCTCCCACAGGTCGTGCACAGCTTTCAAAAGCGGTGGTGTTcctgt                                                                                                                             >  1:5336506‑M2/1‑145 (MQ=255)
                           ggCGCAGGGCCTGGAGGTCGATCAGGGCGCGGGCGGGACGCATAGCGGCTGCCTTATGGTGGTTCGGGTTGAAATCGGTGGTGGATTATTCGCGGGCTTGCCCGCTCCCACAGGTCGTGCACAGCTTTCAAAAGCGGTGGTGTTcctgtg                                                                                                                            >  2:3498995‑M2/1‑144 (MQ=255)
                            tatctcccgttatggcaatagccgggaagtggccaagcacgggcaaagcgcgccggctcaccagcagtaacagcctgtgccggcccatttGCGGGCTTGCCCGCTCCCACAGGTCGTGCACAGCTTTCAAAAACGGTGGTGTTCCTGTgg                                                                                                                           >  2:8028323‑M1/91‑150 (MQ=255)
                            gcgcAGGGCCTGGAGGTCGATCAGGGCGCGGGCGGGACGCATAGCGGCTGCCTTATGGTGGTTCGGGTTGAAATCGGTGGTGGATTATTCGCGGGCTTGCCCGCTCCCACAGGTCGTGCACAGCTTTCAAAATCGGTGGTGGTcctgtgg                                                                                                                           <  1:8979989‑M2/150‑8 (MQ=255)
                            gcgcAGGGCCTGGAGGTCGATCAGGGCGCGGGCGGGACGCATAGCGGCTGCCTTATGGTGGTTCGGGTTGAAATCGGTGGTGGATTATTCGCGGGCTTGCCCGCTCCCACAGGTCGTGCACAGCTTTCAAAAGCGGTGGTGTTcctgtgg                                                                                                                           >  2:6460357‑M2/1‑143 (MQ=255)
                            gcgcAGGGCCTGGAGGTCGATCAGGGCGCGGGCGGGACGCATAGCGGCTGCCTTATGGTGGTTCGGGTTGAAATCGGTGGTGGATTATTCGCGGGCTTGCCCGCTCCCACAGGTCGTGCACAGCTTTCAAAAGCGGTGGTGTTcctgtgg                                                                                                                           >  1:8440257‑M2/1‑143 (MQ=255)
                            gcgcAGGGCCTGGAGGTCGATCAGGGCGCGGGCGGGACGCATAGCGGCTGCCTTATGGTGGTTCGGGTTGAAATCGGTGGTGGATTATTCGCGGGCTTGCCCGCTCCCACAGGTCGTGCACAGCTTTCAAAAGCGGTGGTGTTcctgtgg                                                                                                                           >  1:42411‑M2/1‑143 (MQ=255)
                            gcgcAGGGCCTGGAGGTCGATCAGGGCGCGGGCGGGACGCATAGCGGCTGCCTTATGGTGGTTCGGGTTGAAATCGGTGGTGGATTATTCGCGGGCTTGCCCGCTCCCACAGGTCGTGCACAGCTTTCAAAAGCGGTGGTGTTcctgtgg                                                                                                                           >  1:2852588‑M2/1‑143 (MQ=255)
                              tctcccgttatggcaatagccgggaagtggccaagcacgggcaaagcgcgccggctcaccagcagtaacagcctgtgccggcccattcGCGGGCTTGCCCGCTCCCACAGGTCGTGCACAGCTTTCAAAAGCGGTGGTGTTCCTGTGGGa                                                                                                                         >  2:7212782‑M1/89‑150 (MQ=255)
                                aGGGCCTGGAGGTCGATCAGGGCGCGGGCGGGACGCATAGCGGCTGCCTTATGGTGGTTCGGGTTGAAATCGGTGGTGGATTATTCGCGGGCTTGCCCGCTCCCACAGGTCGTGCACAGCTTTCAAAAGCGGTGGTGTTcctgtgggagc                                                                                                                       >  1:637377‑M2/1‑139 (MQ=255)
                                  ggCCTGGAGGTCGATCAGGGCGCGGGCGGGACGCATAGCGGCTGCCTTATGGTGGTTCGGGTTGAAATCGGTGGTGGATTATTCGCGGGCTTGCCCGCTCCCACAGGTCGTGCACAGCTTTCAAAAGCGGTGGTGTTcctgtgggagcgg                                                                                                                     >  2:4477625‑M2/1‑137 (MQ=255)
                                  ggCCTGGAGGTCGATCAGGGCGCGGGCGGGACGCATAGCGGCTGCCTTATGGTGGTTCGGGTTGAAATCGGTGGTGGAGTATTCGCGGGCTTGCCCGCTCCCACAGGTCGTGCACAGCTTTCAAAAGCGGTGGTGTTcctgtgggagcgg                                                                                                                     >  2:8398340‑M2/1‑137 (MQ=255)
                                   gCCTGGAGGTCGATCAGGGCGCGGGCGGGACGCATAGCGGCTGCCTTATGGTGGTTCGGGTTGAAATCGGTGGTGGATTATTCGCGGGCTTGCCCGCTCCCACAGGTCGTGCACAGCTTTCAAAAGCGGTGGTGTTcctgtgggagcggg                                                                                                                    >  1:7388047‑M2/1‑136 (MQ=255)
                                      tGGAGGTCGATCAGGGCGCGGGCGGGACGCATAGCGGCTGCCTTATGGTGGTTCGGGTTGAAATCGGTGGTGGATTATTCGCGGGCTTGCCCGCTCCCACAGGTCGTGCACAGCTTTCAAAAGCGGTGGTGTTcctgtgggagcgggcaa                                                                                                                 >  1:5015910‑M2/1‑133 (MQ=255)
                                        gAGGTCGATCAGGGCGCGGGCGGGACGCATAGCGGCTGCCTTATGGTGGTTCGGGTTGAAATCGGTGGTGGATTATTCGCGGGCTTGCCCGCTCCCACAGGTCGTGCACAGCTTTCAAAAGCGGTGGTGTTcctgtgggagcgggcttgc                                                                                                               >  1:3795188‑M2/1‑131 (MQ=255)
                                         aGGTCGATCAGGGCGCGGGCGGGACGCATAGCGGCTGCCTTATGGTGGTTCGGGTTGAAATCGGTGGTGGATTATTCGCGGGCTTGCCCGCTCCCACAGGTCGTGCACAGCTTTCAAAAGCGGTGGTGTTcctgtgggagcgggcaagcc                                                                                                              <  2:4072470‑M2/150‑21 (MQ=255)
                                          ggTCGATCAGGGCGCGGGCGGGACGCATAGCGGCTGCCTTATGGTGGTTCGGGTTGAAATCGGTGGTGGATTATTCGCGGGCTTGCCCGCTCCCACAGGTCGTGCACAGCTTTCAAAAGCGGTGGTGTTcctgtgggagcgggcaagccc                                                                                                             >  2:8445824‑M2/1‑129 (MQ=255)
                                            tCGATCCGGGCGGGGGCGGGACGCATAGGGGCTGCCTTATGGTGGTTCGGGTTGAAATCGGTGGTGGATTATTCGCGGGCTTGCCCGCTCCCACAGGTCGTGCACAGCTTCCAAAAGCGGTGGTGTTcctgtgggagggggcaagcccgc                                                                                                           >  2:7751646‑M2/1‑127 (MQ=255)
                                              gATCAGGGCGCGGGCGGGACGCATAGCGGCTGCCTTATGGTGGTTCGGGTTGAAATCGGTGGTGGATTATTCGCGGGCTTGCCCGCTCCCACAGGTCGTGCACAGCTTTCAAAAGCGGTGGTGTTCCTGTGGGAGCGGGCAAGCCCGCaa                                                                                                         >  1:8263026/1‑150 (MQ=255)
                                                tCAGGGCGCGGGCGGGACGCATAGCGGCTGCCTTATGGTGGTTCGGGTTGAAATCGGTGGTGGATTATTCGCGGGCTTGCCCGCTCCCACAGGTCGTGCACAGCTTTCAAAAGCGGTGGTGTTCCTGTGGGAGCGGGCAAGCCCGCaaag                                                                                                       >  1:1446118/1‑149 (MQ=255)
                                                    ggCGCGGGCGGGACGCATAGCGGCTGCCTTATGGTGGTTCGGGTTGAAATCGGTGGTGGATTATTCGCGGGCTTGCCCGCTCCCACAGGTCGTGCACAGCt                                                                                                                                                    <  2:9043204/101‑1 (MQ=255)
                                                    ggCGCGGGCGGGACGCATAGCGGCTGCCTTATGGTGGTTCGGGTTGAAATCGGTGGTGGATTATTCGCGGGCTTGCCCGCTCCCACAGGTCGTGCACAGCTTTCAAAAGCGGTGGTGTTCCTGTGGGAGCGGGCAAGCCCGCAAATGGGc                                                                                                   >  1:9068316/1‑150 (MQ=255)
                                                                    aTAGCGGCTGCCTTATGGTGGTTCGGGTTGAAATCGGTGGTGGATTATTCGCGGGCTTGCCCGCTCCCACAGGTCGTGCACAGCTTTCAAAAGCGGTGGTGTTCCTGTGGGAGCGGGCAAGCCCGCAAATGGGCCGGCACAGGCTGTTAc                                                                                   >  1:8207691/1‑150 (MQ=255)
                                                                    aTAGCGGCTGCCTTATGGTGGTTCGGGTTGAAATCGGTGGTGGATTATTCGCGGGCTTGCCCGCTCCCACAGGTCGTGCACAGCTTTCAAAAGCGGTGGTGTTCCTGTGGGAGCGGGCAAGCCCGCAAATGGGCCGGCACAGGCTGTTAc                                                                                   >  2:5007951/1‑150 (MQ=255)
                                                                     tAGCGGCTGCCTTATGGTGGTTCGGGTTGAAATCGGTGGTGGATTATTCGCGGGCTTGCCCGCTCCCACAGGTCGTGCACAGCTTTCAAAAGCGGTGGTGTTCCTGTGGGAGCGGGCAAGCCCGCAAATGGGCCGGCACAGGCTGTTAct                                                                                  >  1:4528224/1‑150 (MQ=255)
                                                                        cGGCTGCCTTATGGTGGTTCGGGTTGAAATCGGTGGTGGATTATTCGCGGGCTTGCCCGCTCCCACAGGTCGTGCACAGCTTTAAAAAGCGGTGGTGTTCCTGTGGGAGCGGGCAAGCCCGCAAATGGGCCGGCACAGGCTGTTActgct                                                                               >  1:551896/1‑150 (MQ=255)
                                                                         ggCTGCCTTTTGGTGGTTGGGGTTGAAATCGGTGGTGGATTATTCGCGGGCTGGCCCGCGCCCACTGGTCGTGCACAGCTTTCAAACGCGGTGGTGTTCCTGTGGGCGCGGGCAAGCCCGCAAATGGGCCGGCACAGGCTGTTActgctg                                                                              <  1:2137194/150‑1 (MQ=255)
                                                                         ggCTGCCTTATGGTGGTTCGGGTTGAAATCGGTGGTGGATTATTCGCGGGCTTGCCCGCTCCCACAGGTCGTGCACAGCTTTCAAAAGGGGTGGTGTTCCTGTGGGAGCGGGCAAGCCCGCAAATGGGCCGGCACAGGCTGTTActgctg                                                                              <  2:1156395/150‑1 (MQ=255)
                                                                         ggCTGCCTTAGGGTGGTGCGGGTTGAAAGCGGGGGTGGATTATTCGCGGGCGGGCCCGCTCCCACCGGTCGTGCACCGCTTTCAAACGCGGTGGTGTGCCTGTGGGAGCGGGCAAGCCCGCAAATGGGCCGGCACAGGCTGTTActgctg                                                                              <  1:2972631/150‑1 (MQ=38)
                                                                          gCTGCCTTATGGTGGTTCGGGTTGAAATCGGTGGTGGATTATTCGCGGGCTTGCCCGCTCCCACAGGTCGTGCACAGCTTTCAAAAGCGGTGGTGTTCCTGTGGGAGCGGGCAAGCCCGCAAATGGGCCGGCACAGGCTGTTACTGCTgg                                                                             >  2:2207088/1‑150 (MQ=255)
                                                                                 cggctcaccagcagtaacagcctgtgccggcccatttGCGGGCTTGCCCGCTCCCACAGGTCGTGCACAGCTTTCAAAAGCGGTGGTGTTCCTGTgg                                                                                                                           <  1:3503973‑M1/60‑1 (MQ=255)
                                                                                            ggtttGAAATCGGTGGTGGATTATTCGCGGGCTTGCCCGCTCCCACAGGTCGTGCACAGCTTTCAAAAGCGGTGGTGTTCCTGTGGGAGCGGGCAAGCCCGCAAATGGGCCGGCACAGGCTGTTACTGCTGGTGAGCCGGCGCGCTTTGc                                                           <  2:4420286/147‑1 (MQ=255)
                                                                                              gTTGAAATCGGTGGTGGATTATTCGCGGGCTTGCCCGCTCCCACAGGTCGTGcaca                                                                                                                                                       >  1:245474/1‑56 (MQ=255)
                                                                                              gTTGAAATCGGTGGTGGATTATTCGCGGGCTTGCCCGCTCCCACAGGTCGTGCACCGCTTTCAAAAGCGGTGGTGTTCCTGTGGGAGCGGGCAAGCCCGCAAATGGGCCGGCACAGGCTGTTACTGCTGGTGAGCCGGCGCGCTTTGccc                                                         <  2:3314856/150‑1 (MQ=255)
                                                                                              gTTGAAATCGGTGGTGGATTATTCGCGGGCTTGCCCGCTCCCACAGGTCGTGCACAGCTTTCAAAAGCGGTGGTGTTCCTGTGGGAGCGGGCAAGCCCGCAAATGGGCCGGCACAGGCGGTTAATGCTGGTGAGCCGGCGCGCTTTTccc                                                         >  1:1921267/1‑150 (MQ=255)
                                                                                               ttGAAATCGGTGGTGGATTATTCGCGGGCTTGCCCGCTCCCACAGGTCGTGCACAGCTTTCAAAAGCGGTGGTGTTCCTGTGGGAGCGGGCAAGCCCGCAAATGGGCCGGCACAGGCTGTTACTGCTGGTGAGCCGGCGCGCTTTGCCCg                                                        >  2:3725747/1‑150 (MQ=255)
                                                                                                     tCGGTGGTGGATTATTCGCGGGCTTGCCCGCTCCCACCGTTCTGGCACAGCTTTCAACCGCGGTGGTGTTCCTGTGGGCGCGGGCAAGCCCGCAAATGGGCCGGCACAGGCTGTTACTGCTGGTGAGCCGGCGCGCTTTGCCCGTGCTTg                                                  <  1:4840251/150‑1 (MQ=38)
                                                                                                        gtggtgGATTATTCGCGGGCTTGCCCGCTCCCACCGGTCGTGCACAGCTTTCAAAAGCGGTGGTGTTCCTGTGGGAGCGGGCAAGCCCGCAAATGGGCCGGCACAGGCTGTTACTGCTGGTGAGCCGGCGCGCTTTGCCCGTGCTTGGcc                                               <  2:7593590/150‑1 (MQ=255)
                                                                                                           gtggATTATTCGCGGGCTTGCCCGCTCCCACAGGTCGTGCCCAGCTTTCAAAAGCGGTGGTGTTCCGGTGGGAGCGGGCAAGCCCGCAAATGGGCCGGCACAGGCTGTTACTGCTGGTGAGCCGGCGCGCTTTGCCCGTGCTTGGCCACt                                            >  2:3910718/1‑150 (MQ=255)
                                                                                                           gtggATTATTCGCGGGCTTGCCCGCTCCCACAGGTCGTGCACAGCTTTCAAAAGCGGTGGTGTTCCGGTGGGAGCGGGCAAGCCCGCAAATGGGCCGGCACAGGCTGTTACTGCTGGTGAGCCGGCGCGCTTTGCCCGTGCTTGGCCACt                                            >  2:3403768/1‑150 (MQ=255)
                                                                                                             ggATTATTCGCGGGCTTGCCCGCTCCCACAGGTCGTGCACAGCTTTCAAAAGCGGTGGTGTTCCTGTGGGAGCGGGCAAGCCCGCAAATGGGCCGGCACAGGCTGTTACTGCTGGTGAGCCGGCGCGCTTTGCCCGTGCTTGGCCACTTc                                          >  2:2540141/1‑150 (MQ=255)
                                                                                                              gATTATTCGCGGGCTTGCCCGCTCCCACAGGTCGTGCACAGCTTTCAAAAGCGGTGGTGTTCCTGTGGGAGCGGGCAAGCCCGCAAATGGGCCGGCACAGGCTGTTACTGCTGGTGAGCCGGCGCGCTTTGCCCGTGCTTGGCCACTTcc                                         >  1:7440221/1‑150 (MQ=255)
                                                                                                               tttttccgcgGGCTTGCCCGCTCCAACTGGTCGTGCACCGCTTCCAAACGGGGTGGTGTTCGTGTGGGCGGGGGCAAGCCCGCAAATGGGCCGGCACAGGCTGTTACTGCTGGTGAGCCGGCGCGCTTTGCCCGTGCTTGGCCACTTccc                                        <  1:2321660/144‑1 (MQ=255)
                                                                                                               attattCGCGGGCTTGCCCGCTCCCACAGGTCGTGCACAGCTTTCAAAAGCGGTGGTGTTCCTGTGGGAGCGGGCAAGCCCGCAAATGGGCCGGCACAGGCTGTTACTGCTGGTGAGCCGGCGCGCTTTGCCCGTGCTTGGCCACTTccc                                        <  1:5986388/150‑1 (MQ=255)
                                                                                                                ttattCGCGGGCTTGCCCGCTCCCACAGGTCGTGCACAGCTTTCAAAAGCGGTGGTGTTCCTGTGGGAGCGGGCAAGCCCGCAAATGGGCCGGCACAGGCTGTTACTGCTGGTGAGCCGGCGCGCTTTGCCCGTGCTTGGCCACTTCCCg                                       >  2:1775393/1‑150 (MQ=255)
                                                                                                                 tattCGCGGGCTTGCCCGCTCCCACAGGTCGTGCACAGCTTTCAAAAGCGGTGGTGTTCCTGTGGGAGCGGGCAAGCCCGCAAATGGGCCGGCACAGGCTGTTACTGCTGGTGAGCCGGCGCGCTTTGCCCGTGCTTGGCCACTTCCCgg                                      <  2:3522161/150‑1 (MQ=255)
                                                                                                                      GCGGGCTTGCCCGCTCCCACAGGTCGTGCACAGCTTTCAAAAGCGGTGGTGTTCCTGTGGGAGCGGGCAAGCCCGCAAATGGGCCGGCACAGGCTGTTACTGCTGGTGAGCCGGCGCGCTTTGCCCGTGCTTGGCCACTTCCCGGCTAtt                                 <  1:3754955‑M1/150‑1 (MQ=255)
                                                                                                                      GCGGGCTTGCCCGCACCCACAGGTCCTCCACACCTTACAAAACCGGTGGTGTACCTGTGGGAGCGGGCAAGCCCGCAAATGGGCCGGCACAGGCTGTTACTGCTGGTGAGCCGGCGCGCATTGCCCGTGCTTGGCCACTTCCCGGCTAtt                                 >  2:1594351‑M1/1‑150 (MQ=255)
                                                                                                                        GGGCTTGCCCGCTCCCACAGGTCGTGCACAGCTTTCAAAAGCGGTGGTGTTCCTGTGGGAGCGGGCAAGCCCGCAAATGGGCCGGCACAGGCTGTTACTGCTGGTGAGCCGGCGCGCTTTGCCCGTGCTTGGCCACTTCCCGGCTATTGc                               <  2:7682924‑M1/150‑1 (MQ=255)
                                                                                                                        GGGCTTGCCCGCTCCCACAGGAACACCACCGCTTTTGAAAGCTGTGCACGACCTGTGGGAGCGGGCAAGCCCGCAAATGGGCCGGCACAGGCTGTTACTGCTGGTGAGCCGGCGCGCTTTGCCCGTGCTTGGCCACTTCCCGGCTATTGc                               <  2:4306601‑M1/150‑1 (MQ=255)
                                                                                                                          GCTTGCCCGCTCCCACAGGTCGTGCACAGCTTTCAAAAGCGGTGGTGTTCCTGTGGGAGCGGGCAAGCCCGCAAATGGGCCGGCACAGGCTGTTACTGCTGGTGAGCCGGCGCGCTTTGCCCGTGCTTGGCCACTTCCCGGCTATTGCCa                             <  2:5831120‑M1/150‑1 (MQ=255)
                                                                                                                          GCTTGCCCGCTCCCACAGGAACACCACCGCTTTTGAAAGCTGTGCACGACCTGTGGGAGCGGGCAAGCCCGCAAATGGGCCGGCACAGGCTGTTACTGCTGGTGAGCCGGCGCGCTTTGCCCGTGCTTGGCCACTTCCCGGCTATTGCCa                             <  1:4730171‑M1/150‑1 (MQ=255)
                                                                                                                          GCTTGCCCGCTCCCACAGGAACACCACCGCTTTTGAAAGCTGTGCACGACCTGTGGGAGCGGGCAAGCCCGCAAATGGGCCGGCACAGGCTGTTACTGCTGGTGAGCCGGCGCGCTTTGCCCGTGCTTGGCCACTTCCCGGCTATTGCCa                             <  2:8551620‑M1/150‑1 (MQ=255)
                                                                                                                           CTTGCCCGCTCCCACAGGTCGTGCACAGCTTTCAAAAGCGGTGGTGTTCCTGTGGGAGCGGGCAAGCCCGCAAATGGGCCGGCACAGGCTGTTACTGCTGGTGAGCCGGCGCGCTTTGCCCGTGCTTGGCCACTTCCCGGCTATTGCCAt                            <  2:7706923‑M1/150‑1 (MQ=255)
                                                                                                                            TTGCCCGCTCCCACAGGTCGTGCACAGCTTTCAAAAGCGGTGGTGTTCCTGTGGGAGCGGGCAAGCCCGCAAATGGGCCGGCACAGGCTGTTACTGCTGGTGAGCCGGCGCGCTTTGCCCGTGCTTGGCCACTTCCCGGCTATTGCCATa                           <  2:1907253‑M1/150‑1 (MQ=255)
                                                                                                                                ccGCACCCACAGGTCGTGCACAGCTTTCAAAAGCGGTGGTGTTcctgtgggagcgggcaagcccgcgaataatccaccaccgatttcaacccgaaccaccataaggcagccgctatgcgtcccgcccgcgccctgatcgacctccaggcc                       <  2:5453387‑M2/150‑108 (MQ=255)
                                                                                                                                CCGCTCCCCCCGGTCGTGCACAGCTTTCAAAAGCGGTGGTGTTCCTGTGGGAGCGGGCAAGCCCGCAAATGGGCCGGCACAGGCTGTTACTGCTGGTGAGCCGGCGCGCTTTGCCCGTGCTTGGCCACTTCCCGGCTATTGCCATAACgg                       >  2:6663288‑M1/1‑150 (MQ=255)
                                                                                                                                 CGCTCCCACAGGTCGTGCACAGCTTTCAAAAGCGGTGGTGTTCCTGTGGGAGCGGGCAAGCCCGCAAATGGGCCGGCACAGGCTGTTACTGCTGGTGAGCCGGCGCGCTTTGCCCGTGCTTGGCCACTTCCCGGCTATTGCCATAACggg                      >  2:4725959‑M1/1‑150 (MQ=255)
                                                                                                                                 CGCTCCCACAGGTCGTGCACAGCTTTCAAAAGCGGTGGTGTTCCTGTGGGAGCGGGCAAGCCCGCAAATGGGCCGGCACAGGCTGTTACTGCTGGTGAGCCGGCGCGCTTTGCCCGTGCTTGGCCACTTCCCGGCTATTGCCATAACggg                      >  2:7295836‑M1/1‑150 (MQ=255)
                                                                                                                                  GCTCCCACTGGGGGTGCACCGCTTTCAAACGCGGTGGTGTGCCTGTGGGAGCGGGCAAGCCCGCAAATGGGCCGGCACAGGCTGTTACTGCTGGTGAGCCGGCGCGCTTTGCCCGTGCTTGGCCACTTCCCGGCTATTGCCATAACGGga                     >  1:3437859‑M1/1‑150 (MQ=255)
                                                                                                                                  GCTCCCACAGGTCGTGCACAGCTTTCAAAAGCGGTGGTGTTCCTGTGGGAGCGGGCAAGCCCGCAAATGGGCCGGCACAGGCTGTTACTGCTGGTGAGCCGGCGCGCTTTGCCCGTGCTTGGCCACTTCCCGGCTATTGCCATAACGGga                     >  2:6090876‑M1/1‑150 (MQ=255)
                                                                                                                                  GCTCCCACAGGTCGTGCACAGCTTTCAAAAGCGGTGGTGTTCCTGTGGGAGCGGGCAAGCCCGCAAATGGGCCGGCACAGGCTGTTACTGCTGGTGAGCCGGCGCGCTTTGCCCGTGCTTGGCCACTTCCCGGCTATTGCCATAACGGga                     >  2:492896‑M1/1‑150 (MQ=255)
                                                                                                                                  GCTCCCACAGGTCGGGCACCCCTTTCAAACGCGGTTGTGTTCGTGTGGGCGGGGGCAAGCCCGCAAATGGGCCGGCACAGGCTGGTACTGCTGGTGAGCCGGCGCGCTTTGCCCGTGCTTGGCCACTTTCCGGCTATTGCCATAACGGga                     >  1:8071941‑M1/1‑150 (MQ=255)
                                                                                                                                    TCCCACTGGTCGGGCACCGCTTTCAAACGGGGTGGTGTTCCTGTGGGAGCGGGCAAGCCCGCAAATGGGCCGGCACAGGCTGTTACTGCTGGTGAGCCGGCGCGCTTTGCCCGTGCTTGGCCACTTCCCGGCTATTGCCATAACGGgaga                   >  1:7274751‑M1/1‑150 (MQ=255)
                                                                                                                                    TCCCACAGGTCGTGCACAGCTTTCAAAAGCGGTGGTGTTCCTGTGGGAGCGGGCAAGCCCGCAAATGGGCCGGCACAGGCTGTTACTGCTGGTGAGCCGGCGCGCTTTGCCCGTGCTTGGCCACTTCCCGGCTATTGCCATAACGGgaga                   <  1:6135101‑M1/150‑1 (MQ=255)
                                                                                                                                    TCCCACAGGTCGTGCACAGCTTTCAAAAGCGGTGGTGTTCCTGTGGGAGCGGGCAAGCCCGCAAATGGGCCGGCACAGGCTGTTACTGCTGGTGAGCCGGCGCGCTTTGCCCGTGCTTGGCCACTTCCCGGCTATTGCCATAACGGgaga                   <  1:6322702‑M1/150‑1 (MQ=255)
                                                                                                                                      ccACAGGTCGTGCACAGCTTTCAAAAGCGGTGGTGTTcctgtgggagcgggcaagcccgcgaataatccaccaccgatttcaacccgaaccaccataaggcagccgctatgcgtcccgcccgcgccctgatcgacctccaggccctgcgc                 <  1:1331071‑M2/150‑114 (MQ=255)
                                                                                                                                      ccACAGGTCGTGCACAGCTTTCAAAAGCGGTGGTGTTcctgtgggagcgggcaagcccgcaaataatccaccaccgatttcaacccgaaccaccataaggcagccgctatgcgtcccgcccgcgccctgatcgacctccaggccctgcgc                 <  1:3161236‑M2/150‑114 (MQ=255)
                                                                                                                                      CCACAGGTCGTGCACAGCTTTCAAAAGCGGTGGTGTTCCTGTGGGAGCGGGCAAGCCCGCAAATGGGCCGGCACAGGCTGTTACTGCTGGTGAGCCGGCGCGCTTTGCCCGTGCTTGGCCACTTCCCGGCTATTGCCATAACGGGAGata                 <  2:6473191‑M1/150‑1 (MQ=255)
                                                                                                                                       cacaGGTCGTGCACAGCTTTCAAAAGCGGTGGTGTTcctgtgggagcgggcaagcccgcaaataatccaccaccgatttcaacccgaaccaccataaggcagccgctatgcgtcccgcccgcgccctgatcgacctccaggccctgcgcc                <  2:589528‑M2/150‑115 (MQ=255)
                                                                                                                                        acaGGTCGTGCACAGCTTTCAAAAGCGGTGGTGTTcctgtgggagcgggcaagcccgcgaataatccaccaccgatttcaacccgaaccaccataaggcagccgctatgcgtcccgcccgcgccctgatcgacctccaggccctgcgcca               >  2:5078958‑M2/1‑35 (MQ=255)
                                                                                                                                        ACAGGTCGTGCACAGCTTTCAAAAGCGGTGGTGTTCCTGTGGGAGCGGGCAAGCCCGCAAATGGGCCGGCACAGGCTGTTACTGCTGGTGAGCCGGCGCGCTTTGCCCGTGCTTGGCCACTTCCCGGCTATTGCCATAAc                         <  2:3771486‑M1/140‑1 (MQ=255)
                                                                                                                                           GGTCGTGCACAGCTTTCAAAAGCGGTGGTGTTCCTGTGGGAGCGGGCAAGCCCGCAAATGGGCCGGCACAGGCTGTTACTGCTGGTGAGCCGGCGCGCTTTGCCCGTGCTTGGCCACTTCCCGGCTATTGCCATAACGGGAGATATCCAg            >  2:6322702‑M1/1‑150 (MQ=255)
                                                                                                                                            GTCGTGCACAGCTTTCAAAAGCGGTGGTGTTCCTGTGGGAGCGGGCAAGCCCGCAAATGGGCCGGCACAGGCTGTTACTGCTGGTGAGCCGGCGCGCTTTGCCCGTGCTTGGCCACTTCCCGGCTATTGCCATAACGGGAGATATCCAgg           <  2:5266707‑M1/150‑1 (MQ=255)
                                                                                                                                                 gCACAGCTTTCAAAAGCGGTGGTGTTCCTGTGGGAGCGGGCAAGCCCGCAAATGGGCCGGCACAGGCTGTTACTGCTGGTGAGCCGGCGCGCTTTGCCCGTGCTTGGCCACTTCCCGGCTATTGCCATAACGGGAGATATCCAGGCCTTc      <  2:1546541/150‑1 (MQ=255)
                                                                                                                                                   acaGCTTTCAAAAGCGGTGGTGTTCCTGTGGGAGCGGGCAAGCCCGCAAATGGGCCGGCACAGGCTGTTACTGCTGGTGAGCCGGCGCGCTTTGCCCGTGCTTGGCCACTTCCCGGCTATTGCCATAACGGGAGATATCCAGGCCTTCgg    >  1:6425472/1‑150 (MQ=255)
                                                                                                                                                   acaGCTTTCAAAAGCGGTGGTGTTCCTGTGGGAGCGGGCAAGCCCGCAAATGGGCCGGCACAGGCTGTTACTGCTGGTGAGCCGGCGCGCTTTGCCCGTGCTTGGCCACTTCCCGGCTATTGCCATAACGGGAGATATCCAGGCCTTCgg    >  2:8771359/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCAGTTCACGGGCCAGACGGTAGTTGTGGCGCAGGGCCTGGAGGTCGATCAGGGCGCGGGCGGGACGCATAGCGGCTGCCTTATGGTGGTTCGGGTTGAAATCGGTGGTGGATTATTCGCGGGCTTGCCCGCTCCCACAGGTCGTGCACAGCTTTCAAAAGCGGTGGTGTTCCTGTGGGAGCGGGCAAGCCCGCAAATGGGCCGGCACAGGCTGTTACTGCTGGTGAGCCGGCGCGCTTTGCCCGTGCTTGGCCACTTCCCGGCTATTGCCATAACGGGAGATATCCAGGCCTTCGGCA  >  NC_002947/6021108‑6021406

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: