Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 408,432 A→G 10.4% G745G (GGT→GGC aceE ← pyruvate dehydrogenase E1 component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947408,4320AG10.4% 255.1 / 14.3 94G745G (GGT→GGCaceEpyruvate dehydrogenase E1 component
Reads supporting (aligned to +/- strand):  ref base A (41/43);  new base G (5/5);  total (46/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

GAGGCGATGACCGGGCCTTTGCGACCGTTCAGGCACTGCTCGACGTAGGTCTGCTGTGGCTTCTGGCCAGGCTTCAGGCGGTTGGCGCGTTCCACGGCCAGGCCGTCGCGACGCAGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGCACCGACGTTGAACTCTTCACGCAGGATCTTCGCGGCTTCGCGGACTTCGCGCAGGATGGTGCCGGAGCCCATCAGCTGTACGTGGTGTGCGGCTTCGCGGGTGTCCTCTTCCAGCAGGTACATGCCCTTGATGATGCCTTCCTCGACACC  >  NC_002947/408283‑408581
                                                                                                                                                     |                                                                                                                                                     
gagGCGATGACCGGGCCTTTGCGACCGTTCAGGCACTGCTCGACGTAGGTCTGCTGTGGCTTCTGGCCAGGCTTCAGGCGGTTGGCGCGTTCCACGGCCAGGCCGTCGCGACGCAGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGCa                                                                                                                                                       <  2:168481/150‑1 (MQ=255)
gagGCGATGACCGGGCCTTTGCGACCGTTCAGGCACTGCTCGACGTAGGTCTGCTGTGGCTTCTGGCCAGGCTTCAGGCGGTTGGCGCGTTCCACGGCCAGGCCGTCGCGACGCAGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGCa                                                                                                                                                       >  2:1849918/1‑150 (MQ=255)
  ggCGCTGACAGGGCCTTTGCGACCGTTCAGGCACTGCTCGACGTAGGTCTGCTGTGGCTTCTGGCCAGGCTTCAGGCGGTTGGCGCGTTCCACGGCCAGGCCGTCGCGACGCAGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGCAcc                                                                                                                                                     >  2:98958/1‑150 (MQ=255)
  ggCGATGACCGGGCCTTTGCGACCGTTCAGGCACTGATCGTCGTAGGTCTGCTGTGGCTTCTGGCCAGGCTTCAGGCGGTTGGCGCGTTCCACGGCCAGGCCGTCGCGACGCAGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGCAcc                                                                                                                                                     <  1:2343167/150‑1 (MQ=255)
   gCGATGACCGGGCCTTTGCGACCGTTCAGGCACTGCTCGACGTAGGTCTGCTGTGGCTTCTGGCCAGGCTTCAGGCGGTTGGCGCGTTCCACGGCCAGGCCGTCGCGACGCAGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGCACCg                                                                                                                                                    >  2:2376717/1‑150 (MQ=255)
   gCGATGACCGGGCCTTTGCGACCGTTCAGGCACTGCTCGACGTAGGTCTGCTGTGGCTTCTGGCCAGGCTTCAGGCGGTTGGCGCGTTCCACGGCCAGGCCGTCGCGACGCAGTTCGTTGAAGCTGGTGACGATCCAAATGTCGGTAcca                                                                                                                                                    >  2:1957118/1‑149 (MQ=255)
     gATGACCGGGCCTTTGCGACCGTTCAGGCACTGCTCGACGTAGGTCTGCTGTGGCTTCTGGCCAGGCTTCAGGCGGTTGGCGCGTTCCACGGCCAGGCCGTCGCGACGCAGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGCACCGAc                                                                                                                                                  <  2:769347/150‑1 (MQ=255)
        gACCGGGCCTTTGCGACCGTTCAGGCACTGCTCGACGTAGGTCTGCTGTGGCTTCTGGCCAGGCTTCAGGCGGTTGGCGCGTTCCACGGCCAGGCCGTCGCGACGCAGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGCACCGACGtt                                                                                                                                               <  2:1644423/150‑1 (MQ=255)
           cGGGCCTTTGCGACCGTTCAGGCACTGCTCGACGTAGGTCTGCTGTGGCTTCTGGCCAGGCTTCAGGCGGTTGGCGCGTTCCACGGCCAGGCCGTCGCGACGCAGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGTGCCGACGtgt                                                                                                                                              >  2:881562/1‑146 (MQ=255)
           cGGGCCTTTGCGACCGTTCAGGCACTGCTCGACGTAGGTCTGCTGTGGCTTCTGGCCAGGCTTCAGGCGGTTGGCGCGTTCCACGGCCAGGCCGTCGCGACGCAGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGTGCCGACGtgt                                                                                                                                              <  1:881562/148‑3 (MQ=255)
           cGGGCCTTTGCGACCGTTCAGGCACTGCTCGACGTAGGTCTGCTGTGGCTTCTGGCCAGGCTTCAGGCGGTTGGCGCGTTCCACGGCCAGGCCGTCGCGACGCAGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGCACCGACGTTGaa                                                                                                                                            >  2:384280/1‑150 (MQ=255)
           cGGGCCTTTGCGACCGTTCAGGCACTGCTCGACGTAGGTCTGCTGTGGCTTCTGGCCAGGCTTCAGGCGGTTGGCGCGTTCCACGGCCAGGCCGTCGCGACGCAGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGCACCGACGTTGaa                                                                                                                                            <  2:717763/150‑1 (MQ=255)
           cGGGCCTTTGCGACCGTTCAGGCACTGCTCGACGTAGGTCTGCTGTGGCTTCTGGCCAGGCTTCAGGCGGTTGGAGCGTGCTACGGCCGCGCCGGCGCTACGCAGTTCGTTGAAGCAGGTGACGCTTCACACGTCGGCACCGACGTTGaa                                                                                                                                            >  1:644708/1‑150 (MQ=255)
            gggCCTTTGCGACCGTTCAGGCACTGCTCGACGTAGGTCTGCTGTGGCTTCTGGCCAGGCTTCAGGCGGTTGGCGCGTTCCACGGCCAGGCCGTCGCGACGCAGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGCACCGACGTTGAAc                                                                                                                                           >  2:1552962/1‑150 (MQ=255)
                cTTTGCGACCGTTCAGGCACTGCTCGACGTAGGTCTGCTGTGGCTTCTGGCCAGGCTTCAGGCGGTTGGCGCGTTCCACGGCCAGGCCGTCGCGACGCAGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGCACCGACGTTGAACTCtt                                                                                                                                       <  2:2423493/150‑1 (MQ=255)
                 tttGCGACCGTTCAGGCACTGCTCGACGTAGGTCTGCTGTGGCTTCTGGCCAGGCTTCAGGCGGTTGGCGCGTTCCACGTCCAGGCCGTCGCGACGCAGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGCACCGACGTTGAACTCTTc                                                                                                                                      >  1:612511/1‑150 (MQ=255)
                 tttGCGACCGTTCAGGCACTGCTCGACGTAGGTCTGCTGTGGCTTCTGGCCAGGCTTCAGGCGGTTGGCGCGTTCCACGGCCAGGCCGTCGCGACGCAGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGCACCGACGTTGAACTCTTc                                                                                                                                      >  2:1481999/1‑150 (MQ=255)
                      gACCGTTCAGGCACTGCTCGACGTAGGTCTGCTGTGGCTTCTGGCCAGGCTTCAGGCGGTTGGCGCGTTCCACGGCCAGGCCGTCGCGACGCAGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGCACCGACGTTGAACTCTTCACGCa                                                                                                                                 <  2:1952808/150‑1 (MQ=255)
                      gACCGTTCAGGCACTGCTCGACGTAGGTCTGCTGTGGCTTCTGGCCAGGCTTCAGGCGGTTGGCGCGTTCCACGGCCAGGCCGTCGCGACGCAGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGCACCGACGTTGAACTCTTCACGCa                                                                                                                                 >  1:1423155/1‑150 (MQ=255)
                      gACCGTTCAGGCACTGCTCGACGTAGGTCTGCTGTGGCTTCTGGCCAGGCTTCAGGCGGTTGGCGCGTTCCACGGCCAGGCCGTCGCGACGCAGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGCACCGACGTTGAACTCTTCACGCa                                                                                                                                 <  1:2046986/150‑1 (MQ=255)
                      gACCGTTCAGGCACTGCTCGACGTAGGTCTGCTGTGGCTTCTGGCCAGGCTTCAGGCGGTTGGCGCGTTCCACGGCCAGGCCGTCGCGACGCAGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGCACCGACGTTGAACTCTTCACGCa                                                                                                                                 <  1:2047208/150‑1 (MQ=255)
                       aCCGTTCAGGCACTGCTCGACGTAGGTCTGCTGTGGCTTCTGGCCAGGCTTCAGGCGGTTGGCGCGTTCCACGGCCAGGCCGTCGCGACGCAGTTCGTTGAAGCTGGTGACGCCCCACACGTCGGCACCGACGTTGAACTCTTCACGCAg                                                                                                                                >  2:1125267/1‑150 (MQ=255)
                        ccGTTCAGGCACTGCTCGACGTAGGTCTGCTGTGGCTTCTGGCCAGGCTTCAGGCGGTTGGCGCGTTCCACGGCCAGGCCGTCGCGACGCAGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGCACCGACGTTGAACTCTTCACGCAgg                                                                                                                               <  2:1334262/150‑1 (MQ=255)
                        ccGTTCAGGCACTGCTCGACGTAGGTCTGCTGTGGCTTCTGGCCAGGCTGCAGGCGGTTGGCGCGTTCCACGGCCAGGCCGTCGCGACGCAGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGCACCGACGTTGAACTCTTCACGCAgg                                                                                                                               <  2:1832262/150‑1 (MQ=255)
                              aGGCACTGCTCGACGTAGGTCTGCTGTGGCTTCTGGCCAGGCTTCAGGCGGTTGGCGCGTTCCACGGCCAGGCCGTCGCGACGCAGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGCACCGACGTTGAACTCTTCACGCAGGATcttc                                                                                                                         <  2:1671449/150‑1 (MQ=255)
                                  aCTGCTCGACGTAGGTCTGCTGTGGCTTCTGGCCAGGCTTCAGGCGGTTGGCGCGTTCCACGGCCAGGCCGTCGCGACGCAGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGCACCGACGTTGAACTCTTCACGCAGGATcttcgcgg                                                                                                                     >  1:133149/1‑150 (MQ=255)
                                  aCTGCTCGACGTAGGTCTGCTGTGGCTTCTGGCCAGGCTTCAGGCGGTTGGCGCGTTCCACGGCCAGGCCGTCGCGACGCAGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGCACCGACGTTGAACTCTTCACGCAGGATcttc                                                                                                                         >  2:916819/1‑146 (MQ=255)
                                  aCTGCTCGACGTAGGTCTGCTGTGGCTTCTGGCCAGGCTTCAGGCGGTTGGCGCGTTCCACGGCCAGGCCGTCGCGACGCAGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGCACCGACGTTGAACTCTTCACGCAGGATcttc                                                                                                                         <  1:916819/146‑1 (MQ=255)
                                    tGCTCGACGTAGGTCTGCTGTGGCTTCTGGCCAGGCTTCAGGCGGTTGGCGCGTTCCACGGCCAGGCCGTCGCGACGCAGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGCACCGACGTTGAACTCTTCACGCAGGATcttcgcggct                                                                                                                   >  2:438993/1‑150 (MQ=255)
                                         ggcGTAGGTCTGCTGTGGCTTCTGGCCATGCTTTAGGCGGTTGGCGCGTTCCACGGCCAGGCCGTCGCGACGCCGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGCACCGACGTTGAACTCTTCACGCAGGATcttcgcggcttcgcg                                                                                                              <  2:1234171/148‑1 (MQ=255)
                                         gACGTAGGTCTGCTGTGGCTTCTGGCCAGGCTTCAGGCGGTTGGCGCGTTCCACGGCCAGGTCGTCGCGACGCAGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGCACCGACGTTGAACTCTTCACGCAGGATcttcgcggcttcgcg                                                                                                              <  2:912242/150‑1 (MQ=255)
                                          aCGTAGGTCTGCTGTGGCTTCTGGCCAGGCTTCAGGCGGTTGGCGCGTTCCACGGCCAGGCCGTCGCGACGCAGTTCGTGGAAGCTGGGGACGCTCCACACGTCAGCACCGACGTTGAACTCTTCACGCAGGATcttcgcggcttcgcgg                                                                                                             >  2:458933/1‑150 (MQ=255)
                                                 tCTGCTGTGGCTTCTGGCCAGGCTTCAGGCGGTTGGCGCGTTCCACGGCCAGGCCGTCGCGACGCAGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGCACCGACGTTGAACTCTTCACGCAGGTTCTTCGCGGCTTCGCGGACTTcgc                                                                                                      >  1:2296181/1‑150 (MQ=255)
                                                 tCTGCTGTGGCTTCTGGCCAGGCTTCAGGCGGTTGGCGCGTTCCACGGCCAGGCCGTCGCGACGCAGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGCACCGACGTTGAACTCTTCACGCAGGATCTTCGCGGCTTCGCGGACTTcgc                                                                                                      >  1:1398763/1‑150 (MQ=255)
                                                       gtgGCTTCTGGCCAGGCTTCAGGCGGTTGGCGCGTTCCACGGCCAGGCCGTCGCGACGCAGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGCACCGACGTTGAACTCTTCACGCAGGATCTTCGCGGCTTCGCGGACTTCGCGCAGGa                                                                                                >  2:2287949/1‑150 (MQ=255)
                                                         ggCTTCTGGCCAGGCTTCAGGCGGTTGGCGCGTTCCACGGCCAGGCCGTCGCGACGCAGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGCACCGACGTTGAACTCTTCACGCAGGATCTTCGCGGCTTCGCGGACTTCGCGCAGGATg                                                                                              <  2:133149/150‑1 (MQ=255)
                                                              cTGGCCAGGCTTCAGGCGGTTGGCGCGTTCCACGGCCAGGCCGTCGCGACGCAGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGCACCGACGTTGAACTCTTCACGCAGGATCTTCGCGGCTTCGCGGACTTCGCGCAGGATGGTGcc                                                                                         <  2:1423155/150‑1 (MQ=255)
                                                                  ccAGGCTTCAGGCGGTTGGCGCGTTCCACGGCCAGGCCGTCGCGACGCAGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGCACCGACGTTGAACTCTTCACGCAGGATCTTCGCGGCTTCGCGGACTTCGCGCAGGATGGTGCCGGAg                                                                                     <  2:2148193/150‑1 (MQ=255)
                                                                  ccAGGCTTCAGGCGGTTGGCGCGTTCCACGGCCAGGCCGTCGCGACGCAGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGCACCGACGTTGAACTCTTCACGCAGGATCTTCGCGGCTTCGCGGACTTCGCGCAGGATGGTGCCGGAg                                                                                     <  1:384280/150‑1 (MQ=255)
                                                                   cAGGCTTCAGGCGGTTGGCGCGTTCCACGGCCAGGCCGTCGCGACGCAGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGCACCGACGTTGAACTCTTCACGCAGGATCTTCGCGGCTTCGCGGACTTCGCGCAGGATGGTGCCGGAGc                                                                                    >  2:1491490/1‑150 (MQ=255)
                                                                   cAGGCTTCAGGCGGTTGGCGCGTTCCACGGCCAGGCCGTCGCGACGCAGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGCACCGACGTTGAACTCTTCACGCAGGATCTTCGCGGCTTCGCGGACTTCGCGCAGGATGGTGCCGGAGc                                                                                    <  2:644708/150‑1 (MQ=255)
                                                                   cAGGCTTCAGGCGGTTGGCGCGTTCCACGGCCAGGCCGTCGCGACGCAGTTCGTTGAAGCTGGTGACGCTCCACAAGTGGGCACCGACGTGGAACTCTTCACGCAGGATCTTCGCGGCTTCGCGGACTTCGCGCAGGATGGTGCCGGAGc                                                                                    >  1:883533/1‑150 (MQ=255)
                                                                      gCTTCAGGCGGTTGGCGCGTTCCACGGCCAGGCCGTCGCGACGCAGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGCACCGACGTTGAACTCTTCACGCAGGATCTTCGCGGCTTCGCGGACTTCGCGCAGGATGGTGCCGGAGCCCa                                                                                 >  2:1570426/1‑150 (MQ=255)
                                                                       cTTCAGGCGGTTGGCGCGTTCCACGGCCAGGCCGTCGCGACGCAGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGCACCGACGTTGAACTCTTCACGCAGGATCTTCGCGGCTTCGCGGACTTCGCGCAGGATGGTGCCGGAGCCCAt                                                                                <  1:2316594/150‑1 (MQ=255)
                                                                           aGGCGGTTGGCGCGTTCCACGGCCAGGCCGTCGCGACGCAGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGCACCGACGTTGAACTCTTCACGCAGGATCTTCGCGGCTTCGCGGACTTCGCGCCGGATGGTGCCGGAGccttacaaa                                                                            >  2:734282/1‑143 (MQ=255)
                                                                           aGGCGGTTGGCGCGTTCCACGGCCAGGCCGTCGCGACGCAGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGCACCGACGTTGAACTCTTCACGCAGGATCTTCGCGGCTTCGCGGACTTCGCGCAGGATGGTGCCGGAGCCCATCAGc                                                                            >  2:1343376/1‑150 (MQ=255)
                                                                                    gcgcGTTCCACGGCCAGGCCGTCGCGACGCAGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGCACCGACGTTGAACTCTTCACGCAGGATCTTCGCGGCTTCGCGGACTTCGCGCAGGATGGTGCCGGAGCCCATCAGCTGTACgtgg                                                                   >  1:2027514/1‑150 (MQ=255)
                                                                                     cgcgTTCCACGTCCAGGCCGTCGCGACGCAGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGCACCGACGTTGAACTCTTCACGCAGGATCTTCGCGGCTTCGCGGACTTCGCGCAGGATGGTGCCGGAGCCCATCAGCTGTACgtggt                                                                  <  2:612511/150‑1 (MQ=255)
                                                                                       cgTTCCACGGCCAGGCCGTCGCGACGCAGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGCACCGACGTTGAACTCTTCACGCAGGATCTTCGCGGCTTCGCGGACTTCGCGCAGGATGGTGCCGGAGCCCATCAGCTGTACGTGgtgt                                                                <  1:1957118/150‑1 (MQ=255)
                                                                                        gTTCCACGGCCAGGCCGTCGCGACGCAGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGCACCGACGTTGAActct                                                                                                                                        <  2:1668929/77‑1 (MQ=255)
                                                                                        gTTCCACGGCCAGGCCGTCGCGACGCAGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGCACCGACGTTGAActct                                                                                                                                        >  1:1668929/1‑77 (MQ=255)
                                                                                          gcgaCGGCCTGGCCGTCGCGACGCAGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGCACCGACGTTGAACTCTTCACGCAGGATCTTCGCGGCTTCGCGGACTTCGCGCAGGATGGTGCCGGAGCCCATCAGCTGTACGTGGTGTGCg                                                             >  2:477932/4‑150 (MQ=255)
                                                                                           ccACGGCCAGGCCGTCGCGACGCAGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGCACCGACGTTGAACTCTTCACGCAGGATcttcgcggcttcgcg                                                                                                              >  2:2098444/1‑100 (MQ=255)
                                                                                           ccACGGCCAGGCCGTCGCGACGCAGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGCACCGACGTTGAACTCTTCACGCAGGATcttcgcggcttcgcg                                                                                                              <  1:2098444/100‑1 (MQ=255)
                                                                                            cACGGCCAGGCCGTCGCGACGCAGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGCACCGACGTTGAACTCTTCACGCAGGATCTTCGCGGCTTCGCGGACTTCGCGCAGGATGGTGCCGGAGCCCATCAGCTGTACGTGGTGTGCGGc                                                           >  1:1313479/1‑150 (MQ=255)
                                                                                             aCGGCCAGGCCGTCGCGACGCAGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGCACCGACGTTGAACTCTTCACGCAGGATCTTCGCGGCTTCGCGGACTTCGCGCAGGATGGTGCCGGAGCCCATCAGCTGTACGTGGTGTGc                                                              <  2:1369410/146‑1 (MQ=255)
                                                                                             aCGGCCAGGCCGTCGCGACGCAGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGCACCGACGTTGAACTCTTCACGCAGGATCTTCGCGGCTTCGCGGACTTCGCGCAGGATGGTGCCGGAGCCCATCAGCTGTACGTGGTGTGc                                                              >  1:1369410/1‑146 (MQ=255)
                                                                                                 ccAGGCCGTCGCGACGCAGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGCACCGACGTTGAACTCTTCACGCAGGATCTTCGCGGCTTCGCGGACTTCGCGCAGGATGGTGCCGGAGCCCATCAGCTGTACGTGGTGTGCGGCTTcgc                                                      >  2:644066/1‑150 (MQ=255)
                                                                                                   aGGCCGTCGCGACGCAGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGCACCGACGTTGAACTCTTCACGCATTTTCTTCGCGGCTTCGCGGACTTCGCGCAGGATGGTGCCGGAGCCCATCAGCTGTACGTGGTGTGCGGCTTCGCgg                                                    <  1:2376717/150‑1 (MQ=255)
                                                                                                    ggCCGTCGCGACGCAGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGCACCGACGTTGAACTCTTCACGCAGGATCTTCGCGGCTTCGCGGACTTCGCGCAGGATGGTGCCGGAGCCCATCAGCTGTACGTGGTGTGCGGCTTcgcg                                                     >  2:2126282/1‑148 (MQ=255)
                                                                                                    ggCCGTCGCGACGCAGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGCACCGACGTTGAACTCTTCACGCAGGATCTTCGCGGCTTCGCGGACTTCGCGCAGGATGGTGCCGGAGCCCATCAGCTGTACGTGGTGTGCGGCTTcgcg                                                     <  1:2126282/148‑1 (MQ=255)
                                                                                                    ggCCGTCGCGACGCAGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGCACCGACGTTGAACTCTTCACGCAGGATCTTCGCGGCTTCGCGGACTTCGCGCAGGATGGTGCCGGAGCCCATCAGCTGTACGTGGTGTGCGGCTTCGCggg                                                   <  2:1116935/150‑1 (MQ=255)
                                                                                                     ggcgccgcgACGCAGTTCGTGGAAGCTGGTGACGATCCACACGTCGGCACGGACGTTGAACTCTTCACGCAGGAGCTTCGCGGCTTCGCGGACTTCGCGCAGGATGGTGCCGGAGCCCATCAGCTGTACGTGGTGTGCGGCTTCGCGGgt                                                  <  2:1161770/145‑1 (MQ=255)
                                                                                                       cGTCGCGACGCAGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGCACCGACGTTGAACTCTTCACGCAGGATCTTCGCGGCTTCGCGGACTTCGCGCAGGATGGTGCCGGAGCCCATCAGCTGTACGTGGTGTGCGGCTTCGCGGgtgt                                                >  2:2387404/1‑150 (MQ=255)
                                                                                                           gcgACGCAGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGCACCGACGTTGAACTCTTCACGCAGGATCTTCGCGGCTTCGCGGACTTCGCGCAGGATGGTGCCGGAGCCCATCAGCTGTACGTGGTGTGc                                                              >  1:1301543/1‑132 (MQ=255)
                                                                                                           gcgACGCAGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGCACCGACGCTGAACTCTTCACGCAGGATCTTCGCGGCTTCGCGGACTTCGCGCAGGATGGTGCCGGAGCCCATCAGCTGTACGTGGTGTGc                                                              <  2:1301543/132‑1 (MQ=255)
                                                                                                             gACGCAGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGCACCGACGTTGAACTCTTCACGCAGGATCTTCGCGGCTTCGCGGACTTCGCGCAGGATGGTGCCGGAGCCCATCAGCTGTACGTGGTGTGCGGCTTCGCGGGTGTCCTCtt                                          >  2:200325/1‑150 (MQ=255)
                                                                                                                    ttCGTTGAAGCTGGTGACGCTCCACACGTCGGCACCGACGTTGAACTCTTCACGCAGGATCTTCGCGGCTTCGCGGACTTGCGGCAGGATGGTGCCGGAGCCCATCAGCTGTACGTGGTGTGCGGCTTCGCGGGTGTCCTCTTCcagcag                                   >  1:1333998/1‑150 (MQ=255)
                                                                                                                     tCGTTGAAGCTGGTGACGCTCCACACGTCGGCACCGACGTTGAACTCTTCACGCAGGATCTTCGCGGCTTCGCGGACTTCGCGCAGGATGGTGCCGGAGCCCATCAGCTGTACGTGGTGTGCGGCTTCGCGGGTGTCCTCTTCCAGCAgg                                  <  1:1481999/150‑1 (MQ=255)
                                                                                                                       gTTGAAGCTGGTGACGCTCCACACGTCGGCACCGACGTTGAACTCTTCACGCAGGATCTTCGCGGCTTCGCGGACTTCGCGCAGGATGGTGCCGGAGCCCATCAGCTGTACGTGGTGTGCGGCTTCGCGGGTGTCCTCTTCCAGCAGGTa                                >  1:2155193/1‑150 (MQ=255)
                                                                                                                                  tGACGCTCCACGCGTCGGCGCCGACGTTGATCTCTTCACGCAGGATGTTCGCGGCTTCGCGGACGTCGCGCAGGATGGTGCCGGAGCCCATCAGCTGTACGTGGTGTGCGGCTTCGCGGGTGTCCTCTTCCAGCAGGTACATGCCCTtga                     <  2:749585/150‑1 (MQ=255)
                                                                                                                                  tGACGCTCCACACGTCGGCACCGACGTTGAACTCTTCACGCAGGATCTTCGCGGCTTCGCGGACTTCGCGCAGGATGGTGCCGGAGCCCATCAGCTGTACGTGGTGTGCGGCTTCGCGGGTGTCCTCTTGAAGCAGGTACATggcttgta                     <  2:524773/150‑8 (MQ=255)
                                                                                                                                  tGACGCTCCACACGTCGGCACCGACGTTGAACTCTTCACGCAGGATCTTCGCGGCTTCGCGGACTTCGCGCAGGATGGTGCCGGAGCCCATCAGCTGTACGTGGTGTGCGGCTTCGCGGGTGTCCTCTTCCAGCAGGTACATGCCCTtga                     <  1:502554/150‑1 (MQ=255)
                                                                                                                                   gACGCTCCACACGTCGGCACCGACGTTGAACTCTTCACGCAGGATCTTCGCGGCTTCGCGGACTTCGCGCAGGATGGTGCCGGAGCCCATCAGCTGTACGTGGTGTGCGGCTTCGCGGGTGTCCTCTTCCAGCAGGTACATGCCCTtgat                    >  2:2077799/1‑150 (MQ=255)
                                                                                                                                   gACGCTCCACACGTCGGCACCGACGTTGAACTCTTCACGCAGGATCTTCGCGGCTTCGCGGACTTCGCGCAGGATGGTGCCGGAGCCCATCAGCTGTACGTGGTGTGCGGCTTCGCGGGTGTCCTCTTCCAGCAGGTACATGCCCTtgat                    <  1:1491490/150‑1 (MQ=255)
                                                                                                                                      gCTCCACACGTCGGCACCGACGTTGAACTCTTCACGCAGGATCTTCGCGGCTTCGCGGACTTCGCGCAGGATGGTGCCGGAGCCCATCAGCTGTACGTGGTGTGCGGCTTCGCGGGTGTCCTCTTCCAGCAGGTACATGCCCTTgatgat                 <  1:1858670/150‑1 (MQ=255)
                                                                                                                                      gCTCCACACGTCGGCACCGACGTTGAACTCTTCACGCAGGATCTTCGCGGCTTCGCGGACTTCGCGCAGGATGGTGCCGGAGCCCATCAGCTGTACGTGGTGTGCGGCTTCGCGGGTGTCCTCTTCCAGCAGGTACATGCCCTTgatgat                 >  1:1512669/1‑150 (MQ=255)
                                                                                                                                       cTCCACACGTCGGCACCGACGTTGAACTCTTCACGCAGGATCTTCGCGGCTTCGCGGACTTCGCGCAGGATGGTGCCGGAGCCCATCAGCTGTACGTGGTGTGCGGCTTCGCGGGTGTCCTCTTCCAGCAGGTACATGCCCTTGatgatg                <  1:70569/150‑1 (MQ=255)
                                                                                                                                        tCCACACGTCGGCACCGACGTTGAACTCTTCACGCAGGATCTTCGCGGCTTCGAGGCTCACGCGCGGTATGGCGCCGAAGACCATCAGCTGTACGTGGTGTGCGGCTTCGCGGGTGTCCTCTTCCAGCAGGTACATGCCCTTGATGATGc               <  1:2387404/150‑1 (MQ=255)
                                                                                                                                         ttcaacGTCGGTGCCGACGTTGAACTCTTCACGCAGGATCTTCGCGGCTTCGCGGACTTCGCGCAGGATGGTGCCGGAGCCCATCAGCTGTACGTGGTGTGCGGCTTCGCGGgt                                                  >  1:1185442/5‑114 (MQ=255)
                                                                                                                                         ttcaacGTCGGTGCCGACGTTGAACTCTTCACGCAGGATCTTCGCGGCTTCGCGGACTTCGCGCAGGATGGTGCCGGAGCCCATCAGCTGTACGTGGTGTGCGGCTTCGCGGgt                                                  <  2:1185442/110‑1 (MQ=255)
                                                                                                                                         ccACACGTCGGCACCGACGTTGAACTCTTCACGCAGGATCTTCGCGGCTTCGAATATTTCGCTGAGGATGGTGCCGGAGCCCATCAGCTGTACGTGGTGTGCGGCTTCGCGGGTGTCCTCTTCCAGCAGGTACATGCCCTTGATGATGcc              <  2:323296/150‑1 (MQ=255)
                                                                                                                                           caacGTCGGTGCCGACGTTGAACTCTTCACGCAGGATCTTCGCGGCTTCGCGGACTTCGCGCAGGATGGTGCCGGAGCCCATCAGCTGTACGTGGTGTGCGGCTTCGCGGgtgt                                                >  1:1319436/3‑114 (MQ=255)
                                                                                                                                           caacGTCGGTGCCGACGTTGAACTCTTCACGCAGGATCTTCGCGGCTTCGCGGACTTCGCGCAGGATGGTGCCGGAGCCCATCAGCTGTACGTGGTGTGCGGCTTCGCGGgtgt                                                <  2:1319436/112‑1 (MQ=255)
                                                                                                                                           caacGTCGGTGCCGACGTTGAACTCTTCACGCAGGATCTTCGCGGCTTCGCGGACTTCGCGCAGGATGGTGCCGGAGCCCATCAGCTGTACGTGGTGTGCGGCTTCGCGGGTGTCCTCTTCCAGCAGGTACATGCCCTTGATGATGCCtt            >  1:2418442/3‑150 (MQ=255)
                                                                                                                                            cacGTCGGCACCGACGTTGAACTCTTCACGCAGGATCTTCGCGGCTTCGCGGACTTCGCGCAGCTCGTCTGCGGAGCCCATCAGCTGTACGTGGTGTGCGGCTTCGCGGGTGTCCTCTTCCAGCAGGTACATGCCCTTGATGATGCCTTc           <  1:644066/150‑1 (MQ=255)
                                                                                                                                            aacGTCGGTGCCGACGTTGAACTCTTCACGCAGGATCTTCGCGGCTTCGCGGACTTCGCGCAGGATGGTGCCGGAGCCCAt                                                                                >  2:2340545/2‑81 (MQ=255)
                                                                                                                                            aacGTCGGTGCCGACGTTGAACTCTTCACGCAGGATCTTCGCGGCTTCGCGGACTTCGCGCAGGATGGTGCCGGAGCCCAt                                                                                <  1:2340545/80‑1 (MQ=255)
                                                                                                                                             acGTCGGCACCGACGTTGAACTCTTCACGCAGGATCTTCGCGGCTTCGCGGACTTCGCGCAGGATGGTGCCGGAGCCCATCAGCTGTACGTGGTGTGCGGCTTCGCGGGTGTCCTCTTCCAGCAGGTACATGCCCTTGATGATGCCTTcc          <  2:1398763/150‑1 (MQ=255)
                                                                                                                                              cGTCGGCACCGACGTTGAACTCTTCACGCAGGATCTTCGCGGCTTCGCGGACTTCGCGCAGGATGGTGCCGGAGCCCATCAGCTGTACGTGgtgt                                                                >  2:1904503/1‑95 (MQ=255)
                                                                                                                                              cGTCGGCACCGACGTTGAACTCTTCACGCAGGATCTTCGCGGCTTCGCGGACTTCGCGCAGGATGGTGCCGGAGCCCATCAGCTGTACGTGgtgt                                                                <  1:1904503/95‑1 (MQ=255)
                                                                                                                                               gTCGGCACCGACGTTGAACTCTTCACGCAGGATCTTCGCGGCTTCGCGGACTTCGCGCAGGATGGTGCCGGAGCCCATCAGCTGTACGTGGTGTGCGGCTTCGCGGGTGTCCTCTTCCAGCAGGTACATGCCCTTGATGATGCCTTCCTc        >  2:681918/1‑150 (MQ=255)
                                                                                                                                               gTCGGCACCGACGTTGAACTCTTCACGCAGGATCTTCGCGGCTTCGCGGACTTCGCGCAGGATGGTGCCGGAGCCCATCAGCTGTACGTGGTGTGCGGCTTCGCGGGTGTCCTCTTCCAGCAGGTACATGCCCTTGATGATGCCTTCCTc        <  1:2077799/150‑1 (MQ=255)
                                                                                                                                               gTCGGCACCGACGTTGAACTCTTCACGCAGGATCTTCGCGGCTTCGCGGACTTCGCGCAGGATGGTGCCGGAGCCCATCAGCTGTACGTGGTGTGCGGCTTCGCGGGTGTCCTCTTCCAGCAGGTACATGCCCTTGATGATGCCTTCCTc        >  1:802825/1‑150 (MQ=255)
                                                                                                                                                   gCACCGACGTTGAACTCTTCACGCAGGATCTTCGCGGCTTCGCGGACTTCGCGCAGGATGGTGCCGGAGCCCATCAGCTGTACGTGGTGTGCGGCTTcgc                                                      <  1:2381258/100‑1 (MQ=255)
                                                                                                                                                   gCACCGACGTTGAACTCTTCACGCAGGATCTTCGCGGCTTCGCGGACTTCGCGCAGGATGGTGCCGGAGCCCATCAGCTGTACGTGGTGTGCGGCTTcgc                                                      >  2:2381258/1‑100 (MQ=255)
                                                                                                                                                     aCCGACGTTGAACTCTTCACGCAGGATCTTCGCGGCTTCGCGGACTTCGCGCAGGATGGTGCGGGAGCCCATCAGCTGTACGTGGTGTGCGGCTTCGCGGGTGTCCTCTTCCAGCAGGTACATGCCCTTGATGATGCCTTCCTCGACAcc  >  1:2416263/1‑150 (MQ=255)
                                                                                                                                                     aCCGACGTTGAACTCTTCACGCAGGATCTTCGCGGCTTCGCGGACTTCGCGCAGGATGGTGCCGGAGCCCATCAGCTGTACGTGGTGTGCGGCTTCGCGGGTGTCCTCTTCCAGCAGGTACATGCCCTTGATGATGCCTTCCTCGACAcc  <  1:681918/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GAGGCGATGACCGGGCCTTTGCGACCGTTCAGGCACTGCTCGACGTAGGTCTGCTGTGGCTTCTGGCCAGGCTTCAGGCGGTTGGCGCGTTCCACGGCCAGGCCGTCGCGACGCAGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGCACCGACGTTGAACTCTTCACGCAGGATCTTCGCGGCTTCGCGGACTTCGCGCAGGATGGTGCCGGAGCCCATCAGCTGTACGTGGTGTGCGGCTTCGCGGGTGTCCTCTTCCAGCAGGTACATGCCCTTGATGATGCCTTCCTCGACACC  >  NC_002947/408283‑408581

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