Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,491,271 A→C 6.3% L55V (TTG→GTG)  clpV ← protein ClpV1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,491,2710AC 6.3% 275.7 / 6.1 95L55V (TTG→GTG) clpVprotein ClpV1
Reads supporting (aligned to +/- strand):  ref base A (44/45);  new base C (3/3);  total (47/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.50e-01

GGTAATGTGCGGCGACAGCGCCGGGGTGCGGGTATTACCTTTCTTGAAGGTGCCCAGGGCGGTGACCAGCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAGGCACAGGGCCGCCGCGCTTTCGGTGGCCGTGCGGCTGGTTTCGTTGAGCTTGGCGAACAGAGACTTCAGGTTCACAAGGCGACCCCATCGTAAGGAAT  >  NC_002947/3491125‑3491417
                                                                                                                                                  |                                                                                                                                                  
ggTAATGTGCGGCGACAGCGCCGGGGTGCGGGTATTACCTTTCTTGAAGGTGCCCAGGGCGGTGACCAGCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTg                                                                                                                                                 >  2:642479/1‑150 (MQ=255)
   aaTGTGCGGCGACAGCGCCGGGGTGCGGGTATTACCTTTCTTGAAGGTGCCCAGGGCGGTGACCAGCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTgcag                                                                                                                                              <  2:951467/150‑1 (MQ=255)
        gcggcgACAGCGCCGGGGTGCGGGTATTACCTTTCTTGAAGGTGCCCAGGGCGGTGACCAGCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTgcagcagca                                                                                                                                         <  2:2033219/150‑1 (MQ=255)
        gcggcgACAGCGCCGGGGTGCGGGTATTACCTTTCTTGAAGGTGCCCAGGGCGGTGACCAGCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTgcagcagca                                                                                                                                         >  1:1774350/1‑150 (MQ=255)
         cggcgACAGCGCCGGGGTGCGGGTATTACCTTTCTTGAAGGTGCCCAGGGCGGTGACCAGCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTGcagcagcag                                                                                                                                        <  1:457510/150‑1 (MQ=255)
           gcgACAGCGCCGGGGTGCGGGTATTACCTTTCTTGAAGGTGCCCAGGGCGGTGACCAGCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGt                                                                                                                                      >  2:1475224/1‑150 (MQ=255)
            cgACAGCGCCGGGGTGCGGGTATTACCTTTCTTGCCGGTGGCCAGGGCGGTGACCAGCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTg                                                                                                                                     <  2:1774350/150‑1 (MQ=255)
             gACAGCGCCGGGGTGCGGGTATTACCTTTCTTGAAGGTGCCCAGGGCGGTGACCAGCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCACCTGCTGCTGCAGTTGc                                                                                                                                    >  2:907784/1‑150 (MQ=255)
              aCAGCGCCGGGGTGCGGGTATTACCTTTCTTGAAGGTGCCCAGGGCGGTGACCAGCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCACCTGCTGCTGCAGTTGCt                                                                                                                                   <  1:910465/150‑1 (MQ=255)
               cAGCGCCGGGGTGCGGGTATTACCTTTCTTGAAGGTGCCCAGGGCGGTGACCAGCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCACCTGCTGCTGCAGTTGCTc                                                                                                                                  <  1:907784/150‑1 (MQ=255)
                  cgcCGGGGTGCGGGTATTACCTTTCTTGAAGGTGCCCAGGGCGGTGACCAGCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTCCAc                                                                                                                               >  2:889740/1‑150 (MQ=255)
                       gggTGCGGGTATTACCTTTCTTGAAGGTGCCCAGGGCGGTGACCAGCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACt                                                                                                                                                  <  1:2029785/126‑1 (MQ=255)
                       gggTGCGGGTATTACCTTTCTTGAAGGTGCCCAGGGCGGTGACCAGCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACt                                                                                                                                                  >  2:2029785/1‑126 (MQ=255)
                        ggTGCGGGTATTACCTTTCTTGAAGGTGCCCAGGGCGGTGACCAGCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTCCACTTCCAc                                                                                                                         <  1:729544/150‑1 (MQ=255)
                         gTGCGGGTATTACCTTTCTTGAAGGTGCCCAGGGCGGTGACCAGCTGCGCTTGCAGGCGTTCGGCGCCGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTCCACTTCCAcc                                                                                                                        <  1:1335203/150‑1 (MQ=255)
                            cGGGTATTACCTTTCTTGAAGGTGCCCAGGGCGGTGACCAGCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCg                                                                                                                     <  1:1316993/150‑1 (MQ=255)
                                  ttACCTTTCTTGAAGGTGCCCAGGGCGGTGACCAGCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTAATgg                                                                                                               >  2:2001986/1‑150 (MQ=255)
                                      cTTTCTTGAAGGTGCCCAGGGCGGTGACCAGCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTAATGGTGCt                                                                                                           >  2:1168373/1‑150 (MQ=255)
                                       tttCTTGAAGGTGCCCAGGGCGGTGACCAGCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCaa                                                                                                                                                    >  1:1977358/1‑108 (MQ=255)
                                       tttCTTGAAGGTGCCCAGGGCGGTGACCAGCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCaa                                                                                                                                                    <  2:1977358/108‑1 (MQ=255)
                                       tttCTTGAAGGTGCCCAGGGCGGTGACCAGCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTcc                                                                                                                                 <  2:2175485/127‑1 (MQ=255)
                                       tttCTTGAAGGTGCCCAGGGCGGTGACCAGCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTcc                                                                                                                                 >  1:2175485/1‑127 (MQ=255)
                                           ttGAAGGTGCCCAGGGCGGTGACCAGCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTc                                                                                                                                  >  2:1260447/1‑122 (MQ=255)
                                           ttGAAGGTGCCCAGGGCGGTGACCAGCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTc                                                                                                                                  <  1:1260447/122‑1 (MQ=255)
                                             gAAGGTGCCCAGGGCGGTGACCAGCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAg                                                                                                    <  2:53288/150‑1 (MQ=255)
                                               aGGTGCCCAGGGCGGTGACCAGCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAGGc                                                                                                  <  1:1168373/150‑1 (MQ=255)
                                                  tGCCCAGGGCGGTGACCAGCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAGGcaca                                                                                               <  1:2145433/150‑1 (MQ=255)
                                                   gCCCAGGGCGGTGACCAGCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAGGCACAg                                                                                              >  1:1417000/1‑150 (MQ=255)
                                                     ccAGGGCGGTGACCAGCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTAATg                                                                                                                <  2:1076450/130‑1 (MQ=255)
                                                     ccAGGGCGGTGACCAGCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTAATg                                                                                                                >  1:1076450/1‑130 (MQ=255)
                                                     ccAGGGCGGTGACCAGCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAGGCACAggg                                                                                            >  1:1695289/1‑150 (MQ=255)
                                                      cAGGGCGGTGACCAGCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAGGCACAGGgc                                                                                           >  2:1295334/1‑150 (MQ=255)
                                                       aGGGCGGTGACCAGCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTgca                                                                                                                                               >  1:870902/1‑97 (MQ=255)
                                                       aGGGCGGTGACCAGCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTgca                                                                                                                                               <  2:870902/97‑1 (MQ=255)
                                                       aGGGCGGTGACCAGCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTCCACTTc                                                                                                                            >  2:790672/1‑116 (MQ=255)
                                                       aGGGCGGTGACCAGCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTCCACTTc                                                                                                                            <  1:790672/116‑1 (MQ=255)
                                                         ggCGGTGACCAGCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAGGCACAGGgccgc                                                                                        <  2:1652845/150‑1 (MQ=255)
                                                            ggTGACCAGCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAGGca                                                                                                 >  1:447246/1‑138 (MQ=255)
                                                            ggTGACCAGCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAGGca                                                                                                 <  2:447246/138‑1 (MQ=255)
                                                              tGACCAGCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAGGCACAGGGCCGCCgcgc                                                                                   >  2:955827/1‑150 (MQ=255)
                                                                 ccAGCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTgcagca                                                                                                                                            >  2:1779463/1‑90 (MQ=255)
                                                                 ccAGCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTgcagca                                                                                                                                            <  1:1779463/90‑1 (MQ=255)
                                                                 ccAGCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAGGCACAGGGCCGCCGCGCttt                                                                                >  2:1936456/1‑150 (MQ=255)
                                                                 ccAGCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAGGCACAGGGCCGCCGCGCttt                                                                                <  1:2002100/150‑1 (MQ=255)
                                                                  cAGCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAGGCACAGGGCCGCCGCGCTTTc                                                                               >  2:1606204/1‑150 (MQ=255)
                                                                   aGCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAGGCACAGGGCCGCCGCGCTTTCg                                                                              <  1:838836/150‑1 (MQ=255)
                                                                    gCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTaa                                                                                                                  <  2:305899/113‑1 (MQ=255)
                                                                    gCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTaa                                                                                                                  >  1:305899/1‑113 (MQ=255)
                                                                    gCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAGGCACAGGGCCGCCGCGCTTTCgg                                                                             >  2:310733/1‑150 (MQ=255)
                                                                      tGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAGGCACAGGGCCGCCGCGCTTTCGGTg                                                                           <  1:642479/150‑1 (MQ=255)
                                                                        cgcTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAGGCACAGGGCCGCCGCGCTTTCGGTGGc                                                                         <  2:1295251/150‑1 (MQ=255)
                                                                         gcTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAGGCACAGGGCCGCCGCGCTTTCGGTGGcc                                                                        >  1:1529510/1‑150 (MQ=255)
                                                                           ttGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAGGCACAGGGCCGCCGCGCTTTCGGTGGCCGt                                                                      >  2:150568/1‑150 (MQ=255)
                                                                               aGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAGGCACAGGgccgcc                                                                                       <  2:988852/129‑1 (MQ=255)
                                                                               aGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAGGCACAGGgccgcc                                                                                       >  1:988852/1‑129 (MQ=255)
                                                                               aGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAGGCACAGGGCCGCCGCGCTTTCGGTGGCCGTGCgg                                                                  >  2:709202/1‑150 (MQ=255)
                                                                                ggCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAGGCACAGGGCCGCCGCGCTTTCGGTGGCCGTGCGGc                                                                 >  1:1315222/1‑150 (MQ=255)
                                                                                   gTTCGACGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAGGCACAGGGCCGCCGCGCTTTCGGTGGCCGTGCGGCTgg                                                              >  1:1118145/1‑150 (MQ=255)
                                                                                     tCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCCCAGCTTCCACCTCGTATTTGTCTGCCTGGAGGCACAGGGCCGGCGCGCTTTCGGTGGCCGTGCGGCTGGtt                                                            <  1:2255752/150‑1 (MQ=255)
                                                                                         cgacgaCCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAGGCACAGGGCCGCCGCGCTTTCGGTGGCCGTGCGGCTGGTTTCGt                                                        <  1:1936456/150‑1 (MQ=255)
                                                                                          gacgacCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAGGCACAGGGCCGCCGCGCTTTCGGTGGCCGTGCGGCTGGTTTCGtt                                                       <  1:310733/150‑1 (MQ=255)
                                                                                             gacCTGGTAGTGGCGCAGGAACGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAGGCACAGGGCCGCCGCGCTTTCGGTGGCCGTGCGGCTGGTTTCGTTGAg                                                    >  2:217199/1‑150 (MQ=255)
                                                                                                cTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAGGCACAGGGCCGCCGCGCTTTCGGTGGCCGTGCGGCTGGTTTCGtt                                                       <  1:1886011/144‑1 (MQ=255)
                                                                                                cTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAGGCACAGGGCCGCCGCGCTTTCGGTGGCCGTGCGGCTGGTTTCGtt                                                       >  2:1886011/1‑144 (MQ=255)
                                                                                                cTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAGGCACAGGGCCGCCGCGCTTTCGGTGGCCGTGCGGCTGGTTTCGTTGAGCtt                                                 >  1:802935/1‑150 (MQ=255)
                                                                                                 tGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAGGCACAGGGCCGCCGCGCTTTCGGTGGCCATGCGGCTGGTTTCGTTGAGCTTg                                                >  2:2338172/1‑150 (MQ=255)
                                                                                                   gTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAGGCACAGGGCCGCCGCGCTTTCGGTGGCCGTGCGGCTGGTTTCGTTGAGCTTGGc                                              <  2:1315222/150‑1 (MQ=255)
                                                                                                   gTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAGGCACAGGGCCGCCGCGCTTTCGGTGGCCGTGCGGCTGGTTTCGTTGAGCTTGGc                                              >  2:1346839/1‑150 (MQ=255)
                                                                                                         gcgcAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAGGCACAGGGCCGCCGCGCTTTCGGTGGCCGTGCGGCTGGTTTCGTTGAGCTTGGCGAACag                                        <  1:955827/150‑1 (MQ=255)
                                                                                                         gcgcAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAGGCACAGGGCCGCCGCGCTTTCGGTGGCCGTGCGGCTGGTTTCGTTGAGCTTGGCGAACag                                        >  2:1113020/1‑150 (MQ=255)
                                                                                                          cgcAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAGGCACAGGGCCGCCGCGCTTTCGGTGGCCGTGCGGCTGGTTTCGTTGAGCTTGGCGAACaga                                       <  1:255000/150‑1 (MQ=255)
                                                                                                           gcAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTAATGGTGc                                                                                                            >  2:1876828/1‑80 (MQ=255)
                                                                                                           gcAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTAATGGTGc                                                                                                            <  1:1876828/80‑1 (MQ=255)
                                                                                                           gcAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAGGCACAGGGCCGCCGCGCTTTCGGTGGCCGTGCGGCTGGTTTCGTTGAGCTTGGCCAACagag                                      >  1:109568/1‑150 (MQ=255)
                                                                                                            cAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAGGCACAGGGCCGCCGCGCTTTCGGTGGCCGTGCGGCTGGTTTCGTTGAGCTTGGCGAACAgaga                                     <  2:1428602/150‑1 (MQ=255)
                                                                                                                  cGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAGGCACAGGGCCGCCGCGCTTTCGGTGGCCGTGCGGCTGGTTTCGTTGAGCTTGGCGAACAGAGACTTCAg                               <  1:1272497/150‑1 (MQ=255)
                                                                                                                       ccAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAGGCACAGGGCCGCCGCGCTTTCGGTGGCCGTGCGGCTGGTTTCGTTGAGCTTGGCGAACAGAGACTTCAGGTTca                          <  1:1346839/150‑1 (MQ=255)
                                                                                                                       ccAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAGGCACAGGGCCGCCGCGCTTTCGGTGGCCGTGCGGCTGGTTTCGTTGAGCTTGGCGAACAGAGACTTCAGGTTca                          <  1:889740/150‑1 (MQ=255)
                                                                                                                           gTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAGGcaca                                                                                               >  2:2115736/1‑77 (MQ=255)
                                                                                                                           gTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAGGcaca                                                                                               <  1:2115736/77‑1 (MQ=255)
                                                                                                                            tCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAGGCACAGGGCCGCCGCGCTTTCGGTGGCCGTGCGGCTGGTTTCGTTGAGCTTGGCGAACAGAGACTTCAGGTTCACAAgg                     >  2:2146427/1‑150 (MQ=255)
                                                                                                                               cTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAGGCACAGGGCCGCCGCGCTTTCGGTGGCCGTGCGGCTGGt                                                             >  1:2335041/1‑107 (MQ=255)
                                                                                                                               cTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAGGCACAGGGCCGCCGCGCTTTCGGTGGCCGTGCGGCTGGt                                                             <  2:2335041/107‑1 (MQ=255)
                                                                                                                                  tCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAGGCACAGGGCAGCCGCGCTTTCGGTGGCCGTGCGGATGGTATCGTTTGTGTTGGCTAACAGAGACTTCAGGTTCACAAGGCGAccc               >  1:1343207/1‑150 (MQ=255)
                                                                                                                                   cGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAGGCACAGGGCCGCCGCGCTTTCGGTGGCCGTGCGGCTGGTTTCGTTGAGCTTGGCGAACAGAGACTTCAg                               <  2:2004242/133‑1 (MQ=255)
                                                                                                                                   cGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAGGCACAGGGCCGCCGCGCTTTCGGTGGCCGTGCGGCTGGTTTCGTTGAGCTTGGCGAACAGAGACTTCAg                               >  1:2004242/1‑133 (MQ=255)
                                                                                                                                    gatgttgtCCAGCACATCAAGCAGCAGGTGATACAGTACCGCATCGTGATGGTTCTCCGACATGCACAGGGCCGCCGCGCGTTCGGTGGCCGTGCAGCTGGTTTCGTTGAGCTTGGCGAACAGAGACTTCAGGTTCACAAGGCGACCCCa             >  2:392663/3‑150 (MQ=255)
                                                                                                                                      tgttgtCGAGCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAGGCACAGGGCCGCCGCGCTTTCGGTGGCCGTGCGGCTGGTTTCGTTGAGCTTGGCGAACAGAGACTTCAGGTTCACAAGGCGACCCCATc           <  2:1118145/150‑1 (MQ=255)
                                                                                                                                        ttgtCGAGCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAgg                                                                                                   <  1:974412/60‑1 (MQ=255)
                                                                                                                                        ttgtCGAGCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAgg                                                                                                   >  2:974412/1‑60 (MQ=255)
                                                                                                                                          gtggAGCACCTGCTGCTGCAGTTGCTCGACTTCCACCTCGTAATGGTGCTCCGACAGGCACAGGGCCGCCGCGCTTTCGGTGGCCGTGCGGCTGGTTTCGTTGAGCTTGGCGAACAGAGACTTCAGGTTCACAAGGCGACCCCATCGTaa       <  1:320259/147‑1 (MQ=255)
                                                                                                                                          gtggAGCACCTGCTGCTGCAGTTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAGGCACAGGGCCGCCGCGCTTTCGGTGGCCGTGCGGCTGGTTTCGTTGAGCTTGGCGAACAGAGACTTCAGGTTCACAAGGCGACCCCATCGTaa       >  1:827925/4‑150 (MQ=255)
                                                                                                                                          gtCGAGTAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAGGCACAGGGCCGCCGCGCTTTCGGTGGCCGTGCGGCTGGTTTCGTTGAGCTTGGCGAACAGAGACTTCAGGTTCACAAGGCGACCCCATCGTaa       <  2:19357/150‑1 (MQ=255)
                                                                                                                                          gtCGAGCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAGGCACAGGGCCGCCGCGCTTTCGGTGGCCGTGCGGCTGGTTTCGTTGAGCTTGGCGAACAGAGACTTCAGGTTCACAAGGCGACCCCATCGTaa       <  2:19243/150‑1 (MQ=255)
                                                                                                                                             gAGCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAGGCACAGGGCCGCCGCGCTTTCGGTGGCCGTGCGGCTGGTTTCGTTGAGCTTGGCGAACAGAGACTTCAGGTTCACAAGGCGACCCCATCGTaa       >  2:367806/1‑147 (MQ=255)
                                                                                                                                             gAGCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAGGCACAGGGCCGCCGCGCTTTCGGTGGCCGTGCGGCTGGTTTCGTTGAGCTTGGCGAACAGAGACTTCAGGTTCACAAGGCGACCCCATCGTaa       <  1:367806/147‑1 (MQ=255)
                                                                                                                                               gCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAGGCACAGGGCCGCCGCGCTTTCGGTGGCCTTGCGGCTGGTTTCGTTGAGCTTGGCGAACAGAGACTTCAGGTTCACAAGGCGACCCCATCGTAAGGAAt  >  1:1333667/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GGTAATGTGCGGCGACAGCGCCGGGGTGCGGGTATTACCTTTCTTGAAGGTGCCCAGGGCGGTGACCAGCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAGGCACAGGGCCGCCGCGCTTTCGGTGGCCGTGCGGCTGGTTTCGTTGAGCTTGGCGAACAGAGACTTCAGGTTCACAAGGCGACCCCATCGTAAGGAAT  >  NC_002947/3491125‑3491417

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