Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,551,327 A→C 5.5% intergenic (+44/+82) groL → / ← pykA chaperonin GroEL/pyruvate kinase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,551,3270AC 5.5% 321.6 / 9.0 116intergenic (+44/+82)groL/pykAchaperonin GroEL/pyruvate kinase II
Reads supporting (aligned to +/- strand):  ref base A (58/51);  new base C (3/4);  total (61/55)
Fisher's exact test for biased strand distribution p-value = 7.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

Only 100 of 116 total aligned reads displayed.

GATGATCACCACCGAAGCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATGGTGT  >  NC_002947/1551181‑1551470
                                                                                                                                                  |                                                                                                                                               
gatgatCACCACCGAAGACATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCTCGCGaa                                                                                                                                              >  2:948969/1‑150 (MQ=255)
 atgatCACCACCGAAGCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAAc                                                                                                                                             >  1:1254319/1‑150 (MQ=255)
 atgatCACCACCGAAGCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGATGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGACCCGGTTCGAGAGaaa                                                                                                                                             >  2:1170693/1‑149 (MQ=255)
 atgatCACCACCGAAGCCATGGCTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGATGGCATGGGTGGCATGGGCGAAATGATGTAAGCCAGCCTTACCCCCCTAAACCATCAAAAAGCCCCGGTTCGAGAGaaa                                                                                                                                             >  2:368659/1‑149 (MQ=255)
     tCACCACCGAAGCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGAATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCggg                                                                                                                                         <  1:691288/150‑1 (MQ=255)
       cccaccGAAGCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCTCGCGAACCGGGGc                                                                                                                                       <  2:1563205/149‑1 (MQ=255)
        ccaccGAAGCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCt                                                                                                                                      >  1:900705/1‑150 (MQ=255)
        ccaccGAAGCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTACATCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGca                                                                                                                                      >  2:187328/1‑149 (MQ=255)
         caccGAAGCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGTCGCGGGGCCCAGTTTGTGGTGCTGGCCAAGTCGCCGGCACCATCAAAAAGCCCCGGTCCGCGAGAACCGGGGttt                                                                                                                                     >  1:860793/1‑150 (MQ=255)
           ccGAAGCCATGGTTGCCGACCTGCCTAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCTCGCGAACCGGGGCttttt                                                                                                                                  <  1:2112619/150‑1 (MQ=255)
            cGAAGCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCttttt                                                                                                                                  <  1:2234807/150‑1 (MQ=255)
              aaGCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCgg                                                                                                                                          >  1:375466/1‑140 (MQ=255)
                gCCATGGTTGCCTACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCTTGATATGGTTGGCATGGGTAGCATGGTCGGCATGATGTAAGCCAGCCTTACCTCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTAtt                                                                                                                              <  2:778425/150‑1 (MQ=255)
                 ccATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTg                                                                                                                             >  2:2207887/1‑150 (MQ=255)
                  cATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTAtg                                                                                                                              <  1:881927/148‑2 (MQ=255)
                  cATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTAtg                                                                                                                              >  2:881927/1‑147 (MQ=255)
                   aTGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCt                                                                                                                           >  1:1169167/1‑150 (MQ=255)
                      gTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTc                                                                                                                        >  1:736085/1‑150 (MQ=255)
                       ttGCCGACCTGCCGGAAGATATGCCAGCTGCTGGCATGCCTGATATGCGTGGCATGCGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCg                                                                                                                       <  1:1870829/150‑1 (MQ=255)
                           cGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGc                                                                                                                   <  2:2177469/150‑1 (MQ=255)
                              ccTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGCGTGGCATGGCTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCctc                                                                                                                >  1:1214634/1‑150 (MQ=255)
                               cTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGTGGCTTTTTTATTGCTGTCGGGGCctct                                                                                                               >  2:2166638/1‑150 (MQ=255)
                                  ccGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCg                                                                                                            <  2:736085/150‑1 (MQ=255)
                                    ggAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCCTGATGTAAGCAAGAGAACCTGCCCTGCACCATCAAAACACACCGGTTCGCGAGAGCCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCgca                                                                                                          >  1:522621/1‑150 (MQ=255)
                                    ggAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCgca                                                                                                          <  1:2207887/150‑1 (MQ=255)
                                    ggAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCgca                                                                                                          >  1:2208376/1‑150 (MQ=255)
                                     gAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCgcag                                                                                                         >  1:1454351/1‑150 (MQ=255)
                                     gAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCCGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCgcag                                                                                                         <  2:1254319/150‑1 (MQ=255)
                                        gACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGcac                                                                                                      <  2:1454351/150‑1 (MQ=255)
                                            aGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAgg                                                                                                  <  2:2299991/150‑1 (MQ=255)
                                             gCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGTCCGTGTCGGAGCACAAGGc                                                                                                 <  2:1154123/150‑1 (MQ=255)
                                              ccAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCt                                                                                                >  1:2177743/1‑150 (MQ=255)
                                                     ctgGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGTCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTg                                                                                         <  1:368659/150‑1 (MQ=255)
                                                      tgGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGa                                                                                        <  1:1920563/150‑1 (MQ=255)
                                                      tgGCATGCCTGATATGGCTGGCATGTGTGGCATCGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGCTCGCGAGAACCGGGGCTTTTTTATTGCGGTCGGCGCCACCTCGCAGCACAAGCGTGCGCttta                                                                                        <  2:272664/150‑3 (MQ=255)
                                                           tGCCTGATATGGGTGGCATGGGTGGCCTTGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCtt                                                                                   >  1:1230381/1‑150 (MQ=255)
                                                           tGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCCGTGCGGTGGCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCtt                                                                                   <  2:1230381/150‑1 (MQ=255)
                                                           tGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCtt                                                                                   <  2:900705/150‑1 (MQ=255)
                                                             ccTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCgca                                                                                                          >  2:2166071/1‑125 (MQ=255)
                                                               tGATATGGGTGGCATGGGTGGCATGGGCGGCATGCTGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGc                                                                               >  2:335077/1‑150 (MQ=255)
                                                               tGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTAtt                                                                                                                              >  2:569546/1‑103 (MQ=255)
                                                               tGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGt                                                                                                                         >  1:1747628/1‑108 (MQ=255)
                                                               tGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGc                                                                               >  1:1218654/1‑150 (MQ=255)
                                                                        tggcatgggtggcatgggCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCt                                                                                                                           <  1:1760656/97‑1 (MQ=255)
                                                                        tggcatgggtggcatgggCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCt                                                                                                                           >  2:1760656/1‑97 (MQ=255)
                                                                          gcatgggtggcatgggCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCtg                                                                             <  2:1838615/141‑1 (MQ=255)
                                                                          gcatgggtggcatgggCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCtg                                                                             >  1:1838615/1‑141 (MQ=255)
                                                                           catgggtggcatgggCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGt                                                                   <  2:1218654/150‑1 (MQ=255)
                                                                           catgggtggcatgggCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCCACATCAGGGGCTTTTTTATTGCTGGCGTGGCCTCTTCTCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGt                                                                   <  1:870884/150‑1 (MQ=255)
                                                                             tgggtggcatgggCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTgg                                                                 <  2:1151001/150‑1 (MQ=255)
                                                                                     atgggCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCg                                                         >  1:686059/1‑150 (MQ=255)
                                                                                                    tAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCg                                          <  2:2177743/150‑1 (MQ=255)
                                                                                                       gccagccTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACg                                       >  2:1891032/1‑150 (MQ=255)
                                                                                                         cagccTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTg                                     >  2:1509195/1‑150 (MQ=255)
                                                                                                          agccTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCtg                                                                             >  1:107949/1‑109 (MQ=255)
                                                                                                                  cccccTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGt                                                                   >  2:1396222/1‑111 (MQ=255)
                                                                                                                  cccccTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGt                                                                   <  1:1396222/111‑1 (MQ=255)
                                                                                                                  cccccTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCtt                            >  1:2188090/1‑150 (MQ=255)
                                                                                                                     ccTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCa                                   >  1:640057/1‑140 (MQ=255)
                                                                                                                     ccTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCa                                   <  2:640057/140‑1 (MQ=255)
                                                                                                                      cTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATg                        <  2:1386427/150‑1 (MQ=255)
                                                                                                                       tGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGc                       <  1:335077/150‑1 (MQ=255)
                                                                                                                        gCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGCTAAACGTACCAACGGATCACCGACGTGCAGGATCTTCATGcc                      >  2:1133109/1‑150 (MQ=255)
                                                                                                                         cACCATCAAAAAGCCCCGGTTCGCGACAACCGGGGCTTTTTAATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCa                     <  1:1904628/150‑1 (MQ=255)
                                                                                                                            cATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGaaa                                                                        <  2:2270892/96‑1 (MQ=255)
                                                                                                                            cATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTg                  <  1:1891032/150‑1 (MQ=255)
                                                                                                                            cATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCCTGAGCCTTGTGCTGTGaaa                                                                        >  1:2270892/1‑96 (MQ=255)
                                                                                                                             aTCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTGTATTGCTGTCGGGACCTCGTAGCAGCACAAGGCTGGGCTTGAGCCTTGTGCTGTGGAAAGATTGGGTCAACCGACGAGCGAATAACAGACATGCAGGAGCTTCATGCCATTgg                 >  2:648878/1‑150 (MQ=255)
                                                                                                                              tCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGGTCTTCATGCGATTgtt                >  1:965971/1‑148 (MQ=255)
                                                                                                                              tCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGt                >  1:737402/1‑150 (MQ=255)
                                                                                                                              tCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGt                <  1:1860969/150‑1 (MQ=255)
                                                                                                                                aaaaaGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGt                                                                   <  1:831854/97‑1 (MQ=255)
                                                                                                                                aaaaaGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGt                                                                   >  2:831854/1‑97 (MQ=255)
                                                                                                                                aaaaaGCCCCGGTTCGCGAGAACCCGCTACGTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGc              <  1:1170693/150‑1 (MQ=255)
                                                                                                                                 aaaaGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGc                       >  2:1569629/1‑140 (MQ=255)
                                                                                                                                 aaaaGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGc                       <  1:1569629/140‑1 (MQ=255)
                                                                                                                                 aaaaGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGcc             >  2:1765028/1‑150 (MQ=255)
                                                                                                                                  aaaGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTgg                                                                 >  2:339717/1‑97 (MQ=255)
                                                                                                                                  aaaGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTgg                                                                 <  1:339717/97‑1 (MQ=255)
                                                                                                                                   aaGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCAc           <  2:1711976/150‑1 (MQ=255)
                                                                                                                                     gCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCgg                                                                                                                      >  1:1669046/1‑41 (MQ=255)
                                                                                                                                     gCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCgg                                                                                                                      <  2:1669046/41‑1 (MQ=255)
                                                                                                                                     gCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTTAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCACGCCATTGGTGCCACCg         <  2:1169167/150‑1 (MQ=255)
                                                                                                                                     gCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCg         >  2:1760535/1‑150 (MQ=255)
                                                                                                                                      ccccGGTTCTCGCGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCtgt                                                                            <  2:62525/82‑1 (MQ=255)
                                                                                                                                      ccccGGTTCTCGCGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCtgt                                                                            >  1:62525/1‑82 (MQ=255)
                                                                                                                                      ccccGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCTGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGa        >  1:1650041/1‑150 (MQ=255)
                                                                                                                                      ccccGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCAt                         <  2:2291596/133‑1 (MQ=255)
                                                                                                                                      ccccGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCAt                         >  1:2291596/1‑133 (MQ=255)
                                                                                                                                       cccGGTTCTCGCGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGc                                                                                      <  2:1953924/71‑1 (MQ=38)
                                                                                                                                       cccGGTTCTCGCGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGc                                                                                      >  1:1953924/1‑71 (MQ=255)
                                                                                                                                        ccGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCCGCCACCGTACGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATg      >  2:1374825/1‑150 (MQ=255)
                                                                                                                                        ccGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATg      >  1:2094644/1‑150 (MQ=255)
                                                                                                                                        ccGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATg      >  2:1328730/1‑150 (MQ=255)
                                                                                                                                         cGGTTCGCGAGAACCGGGGCTTTTTTCTGTGAGGCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATgg     >  1:1814101/1‑150 (MQ=255)
                                                                                                                                         cGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATgg     <  2:1650041/150‑1 (MQ=255)
                                                                                                                                          ggTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATGgt    <  1:1509195/150‑1 (MQ=255)
                                                                                                                                           gTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGaaa                                                                        >  2:518205/1‑81 (MQ=255)
                                                                                                                                           gTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGaaa                                                                        <  1:518205/81‑1 (MQ=255)
                                                                                                                                           gTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATGgtg   <  2:522621/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
GATGATCACCACCGAAGCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATGGTGT  >  NC_002947/1551181‑1551470

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