Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,831,819 G→A 5.6% R219Q (CGG→CAG)  yfiP → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,831,8190GA 5.6% 270.4 / 4.8 89R219Q (CGG→CAG) yfiPhypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (41/43);  new base A (2/3);  total (43/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.79e-01

CCGCCTGCGGCGCTCCAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGCTTGATCTCCCTTGAC  >  NC_002947/1831670‑1831968
                                                                                                                                                     |                                                                                                                                                     
ccGCCTGCGGCGCTCCAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAgcggcg                                                                                                                                                       <  1:1136068/150‑1 (MQ=255)
  gCCTGCGGCGCTCCGAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCggcggc                                                                                                                                                     >  1:1101857/1‑150 (MQ=255)
  gCCTGCGGCGCTCCAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCggcggc                                                                                                                                                     >  2:306450/1‑150 (MQ=255)
            cTCCAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCAGCCGCCGCTTgc                                                                                                                                           <  1:1855328/150‑2 (MQ=255)
              acAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGAGGTGCCTGGACCTTGCTGGCGATGGGCAGGCTTCGCAGGCGCTGGATGCGTATCTGGATGTGTTCAGCCTGCATTACCTTAGCGGCAAGCGGCGGCTGCCGCTTGACg                                                                                                                                         >  2:190640/2‑150 (MQ=255)
                  gcgcGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCa                                                                                                                                     <  2:479148/150‑1 (MQ=255)
                   cgcgATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAg                                                                                                                                    <  2:1399252/150‑1 (MQ=255)
                      gatgatCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGgac                                                                                                                                 >  2:2185217/1‑150 (MQ=255)
                             aTTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGAcaccc                                                                                                                          >  2:701123/1‑150 (MQ=255)
                             aTTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGAcaccc                                                                                                                          >  2:1780814/1‑150 (MQ=255)
                             aTTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATACCTATCTAGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGAcaccc                                                                                                                          >  2:188630/1‑150 (MQ=255)
                              tttCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCAGCCGCCGCTTgcc                                                                                                                                          >  2:1497742/1‑131 (MQ=255)
                              tttCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCAGCCGCCGCTTgcc                                                                                                                                          <  1:1497742/133‑3 (MQ=255)
                                     aCGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGt                                                                                                                  <  2:1101857/150‑1 (MQ=255)
                                      cGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGtt                                                                                                                 <  1:1093728/150‑1 (MQ=255)
                                       ggCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGTCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGttt                                                                                                                <  1:1594788/150‑1 (MQ=255)
                                        gCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTg                                                                                                               <  1:1780814/150‑1 (MQ=255)
                                        gCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTg                                                                                                               <  2:1150024/150‑1 (MQ=255)
                                         cAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGga                                                                                                                                  >  1:1472178/1‑130 (MQ=255)
                                         cAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGga                                                                                                                                  <  2:1472178/130‑1 (MQ=255)
                                             aGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATAc                                                                                                          <  1:934921/150‑1 (MQ=255)
                                              gTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATAcc                                                                                                         <  1:2185217/150‑1 (MQ=255)
                                               ttGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCt                                                                                                        >  1:2042247/1‑150 (MQ=255)
                                                 gcggcgATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCATCAGCGTTTGCATACCTTc                                                                                                      <  1:1355604/150‑1 (MQ=255)
                                                 gcggcgATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTc                                                                                                      <  1:306450/150‑1 (MQ=255)
                                                    gcgATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTg                                                                                                               >  2:2336996/1‑138 (MQ=255)
                                                    gcgATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTg                                                                                                               <  1:2336996/138‑1 (MQ=255)
                                                    gcgATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCta                                                                                                   >  1:1862886/1‑150 (MQ=255)
                                                      gATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCtata                                                                                                 >  2:5592/1‑150 (MQ=255)
                                                      gATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTCAGCGGCCAGCGGCGGCTGCCGCTTGACGATCAGGATATCCCCCTTCACCGTTTGCATACCTTCCgata                                                                                                 <  2:1780324/150‑1 (MQ=255)
                                                       aTGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGTATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAg                                                                                                >  1:766181/1‑150 (MQ=255)
                                                           gCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGc                                                                                                                                               >  1:1535629/1‑99 (MQ=255)
                                                           gCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGc                                                                                                                                               <  2:1535629/99‑1 (MQ=255)
                                                           gCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACccacca                                                                                                                      <  2:684471/124‑1 (MQ=255)
                                                           gCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACccacca                                                                                                                      >  1:684471/1‑124 (MQ=255)
                                                           gCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCa                                                                                            >  2:1868374/1‑150 (MQ=255)
                                                            ccTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGtt                                                                                                                 <  1:1615196/128‑1 (MQ=255)
                                                            ccTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGtt                                                                                                                 >  2:1615196/1‑128 (MQ=255)
                                                            ccTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCaa                                                                                           <  1:1868374/150‑1 (MQ=255)
                                                             cTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAg                                                                                          >  1:447667/1‑150 (MQ=255)
                                                                 aCCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGc                                                                                      <  2:2042247/150‑1 (MQ=255)
                                                                  ccTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGca                                                                                     >  1:2106769/1‑150 (MQ=255)
                                                                     tGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACaa                                                                                  <  2:763683/150‑1 (MQ=255)
                                                                       cTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAAcc                                                                                >  1:202461/1‑150 (MQ=255)
                                                                        tGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAAccc                                                                               <  1:140907/150‑1 (MQ=255)
                                                                        tGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAAccc                                                                               <  2:1862886/150‑1 (MQ=255)
                                                                             aTGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCt                                                                          >  1:1965998/1‑150 (MQ=255)
                                                                              tGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCtt                                                                         <  2:2127008/150‑1 (MQ=255)
                                                                                gCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCg                                                                                                                   <  2:2338556/106‑1 (MQ=255)
                                                                                gCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCg                                                                                                                   >  1:2338556/1‑106 (MQ=255)
                                                                                 ccAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCt                                                                                                    >  1:382631/1‑120 (MQ=255)
                                                                                 ccAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCt                                                                                                    <  2:382631/120‑1 (MQ=255)
                                                                                 ccAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTgg                                                                      >  2:420522/1‑150 (MQ=255)
                                                                                  cAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGaca                                                                                                                             >  1:1507426/1‑94 (MQ=255)
                                                                                  cAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGaca                                                                                                                             <  2:1507426/94‑1 (MQ=255)
                                                                                  cAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCt                                                                                                        <  2:963467/115‑1 (MQ=255)
                                                                                  cAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCt                                                                                                        >  1:963467/1‑115 (MQ=255)
                                                                                  cAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGt                                                                     >  1:2070163/1‑150 (MQ=255)
                                                                                   aGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGtt                                                                    <  2:2106769/150‑1 (MQ=255)
                                                                                     gCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCa                                                                  <  2:202461/150‑1 (MQ=255)
                                                                                       ttCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATa                                                                <  2:2070163/150‑1 (MQ=255)
                                                                                                gcTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCg                                                                                                                   <  2:733193/90‑1 (MQ=255)
                                                                                                gcTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCg                                                                                                                   >  1:733193/1‑90 (MQ=255)
                                                                                                    gATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCa                                                                                                                                     >  2:2093183/1‑68 (MQ=255)
                                                                                                    gATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCa                                                                                                                                     <  1:2093183/68‑1 (MQ=255)
                                                                                                          tATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCa                                             >  1:1436024/1‑150 (MQ=255)
                                                                                                           aTCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGt                                                                     <  2:881828/125‑1 (MQ=255)
                                                                                                           aTCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGt                                                                     >  1:881828/1‑125 (MQ=255)
                                                                                                            tCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAgg                                           >  1:507336/1‑150 (MQ=255)
                                                                                                                 aTGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTggg                                      >  1:1342077/1‑150 (MQ=255)
                                                                                                                  tgtgTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGCCCCGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGc                                     <  1:1675368/150‑1 (MQ=255)
                                                                                                                    tgtTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCCGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAg                                            >  1:2100680/1‑141 (MQ=255)
                                                                                                                    tgtTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCCGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAg                                            <  2:2100680/141‑1 (MQ=255)
                                                                                                                      ttCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGAcacc                                                                                                                           <  2:298036/60‑1 (MQ=255)
                                                                                                                      ttCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGAcacc                                                                                                                           >  1:298036/1‑60 (MQ=255)
                                                                                                                      ttCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCgaga                                 <  2:1965998/150‑1 (MQ=255)
                                                                                                                      ttCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCgaga                                 >  1:2324721/1‑150 (MQ=255)
                                                                                                                              gCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAggg                         >  1:366736/1‑150 (MQ=255)
                                                                                                                                 ttACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGa                      <  1:190640/150‑1 (MQ=255)
                                                                                                                                    ccTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCa                            <  1:1469994/141‑1 (MQ=255)
                                                                                                                                    ccTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCa                            >  2:1469994/1‑141 (MQ=255)
                                                                                                                                      tGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGc                 >  1:1041949/1‑150 (MQ=255)
                                                                                                                                      tGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGAAGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGc                 >  1:1118385/1‑150 (MQ=255)
                                                                                                                                       gAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGCt                <  1:701123/150‑1 (MQ=255)
                                                                                                                                        aGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGCtt               <  1:188630/150‑1 (MQ=255)
                                                                                                                                          cGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGc                                                                                                                    <  1:1256881/47‑1 (MQ=255)
                                                                                                                                          cGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGc                                                                                                                    >  2:1256881/1‑47 (MQ=255)
                                                                                                                                          cGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGc                 <  2:1118715/146‑1 (MQ=255)
                                                                                                                                          cGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGc                 >  1:1118715/1‑146 (MQ=255)
                                                                                                                                           gtcAAGCGGCAGCCGCCGCTTGCCGAGCAGGACGACACCCACCAGCGTTTGCATACCTTc                                                                                                      >  2:2010948/3‑60 (MQ=25)
                                                                                                                                           gtcAAGCGGCAGCCGCCGCTTGCCGAGCAGGACGACACCCACCAGCGTTTGCATACCTTc                                                                                                      <  1:2010948/58‑1 (MQ=25)
                                                                                                                                            gCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGCTTGAtc           >  2:1406765/1‑150 (MQ=255)
                                                                                                                                                  ggcggcTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGAtt                    <  2:291604/135‑1 (MQ=255)
                                                                                                                                                  ggcggcTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGAtt                    >  1:291604/1‑135 (MQ=255)
                                                                                                                                                    cggcTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCCGCCTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGCTTGATCTCCCTTGa   <  2:1247607/150‑1 (MQ=255)
                                                                                                                                                     ggcTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGCTTGATCTCCCTTGAc  <  1:5592/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCGCCTGCGGCGCTCCAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGCTTGATCTCCCTTGAC  >  NC_002947/1831670‑1831968

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: