Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 408,431 C→T 10.3% A746T (GCC→ACC)  aceE ← pyruvate dehydrogenase E1 component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947408,4310CT10.3% 249.8 / 15.2 96A746T (GCC→ACC) aceEpyruvate dehydrogenase E1 component
Reads supporting (aligned to +/- strand):  ref base C (41/45);  new base T (6/4);  total (47/49)
Fisher's exact test for biased strand distribution p-value = 5.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01

GAGGCGATGACCGGGCCTTTGCGACCGTTCAGGCACTGCTCGACGTAGGTCTGCTGTGGCTTCTGGCCAGGCTTCAGGCGGTTGGCGCGTTCCACGGCCAGGCCGTCGCGACGCAGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGCACCGACGTTGAACTCTTCACGCAGGATCTTCGCGGCTTCGCGGACTTCGCGCAGGATGGTGCCGGAGCCCATCAGCTGTACGTGGTGTGCGGCTTCGCGGGTGTCCTCTTCCAGCAGGTACATGCCCTTGATGATGCCTTCCTC  >  NC_002947/408283‑408575
                                                                                                                                                    |                                                                                                                                                
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                                                                                                                                        <  1:881562/148‑3 (MQ=255)
           cGGGCCTTTGCGACCGTTCAGGCACTGCTCGACGTAGGTCTGCTGTGGCTTCTGGCCAGGCTTCAGGCGGTTGGCGCGTTCCACGGCCAGGCCGTCGCGACGCAGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGTGCCGACGtgt                                                                                                                                        >  2:881562/1‑146 (MQ=255)
           cGGGCCTTTGCGACCGTTCAGGCACTGCTCGACGTAGGTCTGCTGTGGCTTCTGGCCAGGCTTCAGGCGGTTGGCGCGTTCCACGGCCAGGCCGTCGCGACGCAGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGCACCGACGTTGaa                                                                                                                                      <  2:717763/150‑1 (MQ=255)
           cGGGCCTTTGCGACCGTTCAGGCACTGCTCGACGTAGGTCTGCTGTGGCTTCTGGCCAGGCTTCAGGCGGTTGGCGCGTTCCACGGCCAGGCCGTCGCGACGCAGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGCACCGACGTTGaa                                                                                                                                      >  2:384280/1‑150 (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                                                                                                                           <  1:2046986/150‑1 (MQ=255)
                      gACCGTTCAGGCACTGCTCGACGTAGGTCTGCTGTGGCTTCTGGCCAGGCTTCAGGCGGTTGGCGCGTTCCACGGCCAGGCCGTCGCGACGCAGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGCACCGACGTTGAACTCTTCACGCa                                                                                                                           <  1:2047208/150‑1 (MQ=255)
                      gACCGTTCAGGCACTGCTCGACGTAGGTCTGCTGTGGCTTCTGGCCAGGCTTCAGGCGGTTGGCGCGTTCCACGGCCAGGCCGTCGCGACGCAGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGCACCGACGTTGAACTCTTCACGCa                                                                                                                           <  2:1952808/150‑1 (MQ=255)
                      gACCGTTCAGGCACTGCTCGACGTAGGTCTGCTGTGGCTTCTGGCCAGGCTTCAGGCGGTTGGCGCGTTCCACGGCCAGGCCGTCGCGACGCAGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGCACCGACGTTGAACTCTTCACGCa                                                                                                                           >  1:1423155/1‑150 (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                                                                               <  1:384280/150‑1 (MQ=255)
                                                                  ccAGGCTTCAGGCGGTTGGCGCGTTCCACGGCCAGGCCGTCGCGACGCAGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGCACCGACGTTGAACTCTTCACGCAGGATCTTCGCGGCTTCGCGGACTTCGCGCAGGATGGTGCCGGAg                                                                               <  2:2148193/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                                                                                                                                  >  1:1668929/1‑77 (MQ=255)
                                                                                        gTTCCACGGCCAGGCCGTCGCGACGCAGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGCACCGACGTTGAActct                                                                                                                                  <  2:1668929/77‑1 (MQ=255)
                                                                                          gcgaCGGCCTGGCCGTCGCGACGCAGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGCACCGACGTTGAACTCTTCACGCAGGATCTTCGCGGCTTCGCGGACTTCGCGCAGGATGGTGCCGGAGCCCATCAGCTGTACGTGGTGTGCg                                                       >  2:477932/4‑150 (MQ=255)
                                                                                           ccACGGCCAGGCCGTCGCGACGCAGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGCACCGACGTTGAACTCTTCACGCAGGATcttcgcggcttcgcg                                                                                                        <  1:2098444/100‑1 (MQ=255)
                                                                                           ccACGGCCAGGCCGTCGCGACGCAGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGCACCGACGTTGAACTCTTCACGCAGGATcttcgcggcttcgcg                                                                                                        >  2:2098444/1‑100 (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:1512669/1‑150 (MQ=255)
                                                                                                                                      gCTCCACACGTCGGCACCGACGTTGAACTCTTCACGCAGGATCTTCGCGGCTTCGCGGACTTCGCGCAGGATGGTGCCGGAGCCCATCAGCTGTACGTGGTGTGCGGCTTCGCGGGTGTCCTCTTCCAGCAGGTACATGCCCTTgatgat           <  1:1858670/150‑1 (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                                          <  2:1319436/112‑1 (MQ=255)
                                                                                                                                           caacGTCGGTGCCGACGTTGAACTCTTCACGCAGGATCTTCGCGGCTTCGCGGACTTCGCGCAGGATGGTGCCGGAGCCCATCAGCTGTACGTGGTGTGCGGCTTCGCGGgtgt                                          >  1:1319436/3‑114 (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                                                          <  1:1904503/95‑1 (MQ=255)
                                                                                                                                              cGTCGGCACCGACGTTGAACTCTTCACGCAGGATCTTCGCGGCTTCGCGGACTTCGCGCAGGATGGTGCCGGAGCCCATCAGCTGTACGTGgtgt                                                          >  2:1904503/1‑95 (MQ=255)
                                                                                                                                               gTCGGCACCGACGTTGAACTCTTCACGCAGGATCTTCGCGGCTTCGCGGACTTCGCGCAGGATGGTGCCGGAGCCCATCAGCTGTACGTGGTGTGCGGCTTCGCGGGTGTCCTCTTCCAGCAGGTACATGCCCTTGATGATGCCTTCCTc  <  1:2077799/150‑1 (MQ=255)
                                                                                                                                               gTCGGCACCGACGTTGAACTCTTCACGCAGGATCTTCGCGGCTTCGCGGACTTCGCGCAGGATGGTGCCGGAGCCCATCAGCTGTACGTGGTGTGCGGCTTCGCGGGTGTCCTCTTCCAGCAGGTACATGCCCTTGATGATGCCTTCCTc  >  1:802825/1‑150 (MQ=255)
                                                                                                                                               gTCGGCACCGACGTTGAACTCTTCACGCAGGATCTTCGCGGCTTCGCGGACTTCGCGCAGGATGGTGCCGGAGCCCATCAGCTGTACGTGGTGTGCGGCTTCGCGGGTGTCCTCTTCCAGCAGGTACATGCCCTTGATGATGCCTTCCTc  >  2:681918/1‑150 (MQ=255)
                                                                                                                                                   gCACCGACGTTGAACTCTTCACGCAGGATCTTCGCGGCTTCGCGGACTTCGCGCAGGATGGTGCCGGAGCCCATCAGCTGTACGTGGTGTGCGGCTTcgc                                                <  1:2381258/100‑1 (MQ=255)
                                                                                                                                                   gCACCGACGTTGAACTCTTCACGCAGGATCTTCGCGGCTTCGCGGACTTCGCGCAGGATGGTGCCGGAGCCCATCAGCTGTACGTGGTGTGCGGCTTcgc                                                >  2:2381258/1‑100 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
GAGGCGATGACCGGGCCTTTGCGACCGTTCAGGCACTGCTCGACGTAGGTCTGCTGTGGCTTCTGGCCAGGCTTCAGGCGGTTGGCGCGTTCCACGGCCAGGCCGTCGCGACGCAGTTCGTTGAAGCTGGTGACGCTCCACACGTCGGCACCGACGTTGAACTCTTCACGCAGGATCTTCGCGGCTTCGCGGACTTCGCGCAGGATGGTGCCGGAGCCCATCAGCTGTACGTGGTGTGCGGCTTCGCGGGTGTCCTCTTCCAGCAGGTACATGCCCTTGATGATGCCTTCCTC  >  NC_002947/408283‑408575

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: