Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,491,276 A→T 5.5% L53Q (CTG→CAG)  clpV ← protein ClpV1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,491,2760AT 5.5% 261.7 / 3.8 84L53Q (CTG→CAG) clpVprotein ClpV1
Reads supporting (aligned to +/- strand):  ref base A (40/39);  new base T (2/3);  total (42/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATGTGCGGCGACAGCGCCGGGGTGCGGGTATTACCTTTCTTGAAGGTGCCCAGGGCGGTGACCAGCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAGGCACAGGGCCGCCGCGCTTTCGGTGGCCGTGCGGCTGGTTTCGTTGAGCTTGGCGAACAGAGACTTCAGGTTCACAAGGCGACCCCATCGTAAGGAATGTGGAAGA  >  NC_002947/3491128‑3491425
                                                                                                                                                    |                                                                                                                                                     
aaTGTGCGGCGACAGCGCCGGGGTGCGGGTATTACCTTTCTTGAAGGTGCCCAGGGCGGTGACCAGCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTgcag                                                                                                                                                      <  2:951467/150‑1 (MQ=255)
     gcggcgACAGCGCCGGGGTGCGGGTATTACCTTTCTTGAAGGTGCCCAGGGCGGTGACCAGCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTgcagcagca                                                                                                                                                 >  1:1774350/1‑150 (MQ=255)
     gcggcgACAGCGCCGGGGTGCGGGTATTACCTTTCTTGAAGGTGCCCAGGGCGGTGACCAGCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTgcagcagca                                                                                                                                                 <  2:2033219/150‑1 (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)
                     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)
                                    tttCTTGAAGGTGCCCAGGGCGGTGACCAGCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTcc                                                                                                                                         <  2:2175485/127‑1 (MQ=255)
                                    tttCTTGAAGGTGCCCAGGGCGGTGACCAGCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTcc                                                                                                                                         >  1:2175485/1‑127 (MQ=255)
                                        ttGAAGGTGCCCAGGGCGGTGACCAGCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTc                                                                                                                                          <  1:1260447/122‑1 (MQ=255)
                                        ttGAAGGTGCCCAGGGCGGTGACCAGCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTc                                                                                                                                          >  2:1260447/1‑122 (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                                                                                                                                                       <  2:870902/97‑1 (MQ=255)
                                                    aGGGCGGTGACCAGCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTgca                                                                                                                                                       >  1:870902/1‑97 (MQ=255)
                                                    aGGGCGGTGACCAGCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTCCACTTc                                                                                                                                    <  1:790672/116‑1 (MQ=255)
                                                    aGGGCGGTGACCAGCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTCCACTTc                                                                                                                                    >  2:790672/1‑116 (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                                                                                        <  1:2002100/150‑1 (MQ=255)
                                                              ccAGCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAGGCACAGGGCCGCCGCGCttt                                                                                        >  2:1936456/1‑150 (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                                                                                               >  1:988852/1‑129 (MQ=255)
                                                                            aGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAGGCACAGGgccgcc                                                                                               <  2:988852/129‑1 (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                                                               >  2:1886011/1‑144 (MQ=255)
                                                                                             cTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAGGCACAGGGCCGCCGCGCTTTCGGTGGCCGTGCGGCTGGTTTCGtt                                                               <  1:1886011/144‑1 (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                                                                                                                    <  1:1876828/80‑1 (MQ=255)
                                                                                                        gcAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTAATGGTGc                                                                                                                    >  2:1876828/1‑80 (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                                                                                                       <  1:2115736/77‑1 (MQ=255)
                                                                                                                        gTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAGGcaca                                                                                                       >  2:2115736/1‑77 (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                                                                                                           >  2:974412/1‑60 (MQ=255)
                                                                                                                                     ttgtCGAGCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAgg                                                                                                           <  1:974412/60‑1 (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               <  1:367806/147‑1 (MQ=255)
                                                                                                                                          gAGCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAGGCACAGGGCCGCCGCGCTTTCGGTGGCCGTGCGGCTGGTTTCGTTGAGCTTGGCGAACAGAGACTTCAGGTTCACAAGGCGACCCCATCGTaa               >  2:367806/1‑147 (MQ=255)
                                                                                                                                            gCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAGGCACAGGGCCGCCGCGCTTTCGGTGGCCTTGCGGCTGGTTTCGTTGAGCTTGGCGAACAGAGACTTCAGGTTCACAAGGCGACCCCATCGTAAGGAAt          >  1:1333667/1‑150 (MQ=255)
                                                                                                                                                  gcagcagcaGGTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAGGCACAGGGCCGCCGCGCTTTCGGTGGCCGTGCGGCTGGTTTCGTTGAGCTTgg                                                       <  2:1319097/99‑1 (MQ=255)
                                                                                                                                                  gcagcagcaGGTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAGGCACAGGGCCGCCGCGCTTTCGGTGGCCGTGCGGCTGGTTTCGTTGAGCTTgg                                                       >  1:1319097/1‑99 (MQ=255)
                                                                                                                                                  gcagcagcaGGTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAGGCACAGGGCCGCCGCGCTTTCGGTGGCCGTGCGGCTGGTTTCGTTGAGCTTGGCGAACAGAGACTTCAGGTTCACAAGGCGACCCCATCGTAAGGAATGTGGaa    >  1:2048101/1‑150 (MQ=255)
                                                                                                                                                  gcagcagcaGGTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAGGCACAGGGCCGCCGCGCTTTCGGTGGCCGTGCGGCTGGTTTCGTTGAGCTTGGCGAACAGAGACTTCAGGTTCACAAGGCGACCCCATCGTAAGGAATGTGGaa    <  1:207536/150‑1 (MQ=255)
                                                                                                                                                  gcagcagcaGGTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAGGCACAGGGCCGCCGCGCTTTCGGTGGCCGTGCGGCTGGTTTCGTTGAGCTTGGCGAACAGAGACTTCAGGTTCACAAGGCGACCCCATCGTAAGGAATGTGGaa    <  1:2146427/150‑1 (MQ=255)
                                                                                                                                                   cagcagcagGTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAGGCTCAGGGCCGCCGCGCTTTCGGTGGCCGTGCGGCTGGTTTCGTTGAGCTTGGCGGACAGAGACTTCAGGTTCACAAGGCGACCCCATCGTAAGGAAGGTGGAAg   >  2:1228310/1‑150 (MQ=255)
                                                                                                                                                    agcagcagGTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAGGCACAGGGCCGCCGCGCTTTCGGTGGCCGTGCGGCTGGTTTCGTTGAGCTTGGCGAACAGAGACTTCAGGTTCACAAGGCGACCCCATCGTAAGGAATGTGGAAGa  <  1:150568/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AATGTGCGGCGACAGCGCCGGGGTGCGGGTATTACCTTTCTTGAAGGTGCCCAGGGCGGTGACCAGCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAACTGCAGCAGCAGGTGCTCCACTTCCACCTCGTAATGGTGCTCCGACAGGCACAGGGCCGCCGCGCTTTCGGTGGCCGTGCGGCTGGTTTCGTTGAGCTTGGCGAACAGAGACTTCAGGTTCACAAGGCGACCCCATCGTAAGGAATGTGGAAGA  >  NC_002947/3491128‑3491425

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: