Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,885,008 C→A 5.0% E139* (GAA→TAA) ‡ yedI ← inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,885,0080CA 5.0% 317.9 / 2.4 98E139* (GAA→TAA) ‡yedIinner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (49/44);  new base A (3/2);  total (52/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 = 9.61e-01

ACCAATGGTCATCACCACCGCAATGCCGGAAAGCACGATGATCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCGAAACCCTCGAAGCACAGGAACGCGCCACCCAGCATCAACAAC  >  NC_002947/4884862‑4885156
                                                                                                                                                  |                                                                                                                                                    
acCAATGGTCATCACCACCGCAATGCCGGAAAGCACGATGATCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCaaa                                                                                                                                                   >  1:598151/1‑150 (MQ=255)
 ccAATGGTCATCACCACCGCAATGCCGGAAAGCACGATGATCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTCCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCaaaa                                                                                                                                                  <  1:88000/150‑1 (MQ=255)
  cAATGGTCATCACCACCGCAATGCCGGAAAGCACGATGATCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCaaaa                                                                                                                                                  >  2:1815607/1‑149 (MQ=255)
  cAATGGTCATCACCACCGCAATGCCGGAAAGCACGATGATCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCaaaa                                                                                                                                                  <  1:1815607/149‑1 (MQ=255)
     tGGTCATCACCACCGCAATGCCGGAAAGCACGATGATCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGcca                                                                                                                                              <  1:2264803/150‑1 (MQ=255)
        tcatcaCCACCGCAATGCCGGAAAGCACGATGATCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGccacca                                                                                                                                           >  1:947649/1‑150 (MQ=255)
          atcaCCACCGCAATGCCGGAAAGCACGATGATCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGccaagatc                                                                                                                                         >  2:2358372/1‑145 (MQ=255)
           ccaccacCGCAATGTCGGAAAGCACGATGATCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGt                                                                                                                                        <  2:1525445/149‑1 (MQ=255)
            caccacCGCAATGCCGGAAAGCACGATGATCTGCGGGGTCAGCGGCGAGTCAGCGAAGATACACAAGGTAATGGCGATGATTTCTGCCGACAGAATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTc                                                                                                                                       >  2:482879/1‑150 (MQ=255)
               caccGCAATGCCGGAAAGCACGATGATCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGGCCTCTTGAGCTGGGCGGGTGCAATATCCACCAGGTCaaa                                                                                                                                    >  2:872506/1‑149 (MQ=255)
                 ccGCAATGCCGGAAAGCACGATGATCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAacac                                                                                                                                  >  1:1098741/1‑150 (MQ=255)
                 ccGCAATGCCGGAAAGCACGATGATCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAacac                                                                                                                                  >  1:771374/1‑150 (MQ=255)
                   gCAATGCCGGAAAGCACGATGATCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTg                                                                                                                                <  1:1038724/150‑1 (MQ=255)
                          cGGAAAGCACGATGCTCTGCTGGGTCAGCGTCGATGCATCCATGATTCCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCgg                                                                                                                         >  1:1863512/1‑150 (MQ=255)
                          cGGAAAGCACGATGATCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGGTGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCgg                                                                                                                         >  1:763932/1‑150 (MQ=255)
                          cGGAAAGCACGATGATCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCgg                                                                                                                         <  2:2210560/150‑1 (MQ=255)
                           ggAAAGCACGATGATCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGc                                                                                                                        >  2:2240588/1‑150 (MQ=255)
                            gAAAGCACGATGATCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCa                                                                                                                       >  2:1881534/1‑150 (MQ=255)
                            gAAAGCACGATGATCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCa                                                                                                                       >  1:626560/1‑150 (MQ=255)
                               aGCACGATGATCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACg                                                                                                                    <  2:1465147/150‑1 (MQ=255)
                                  aCGATGATCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGcc                                                                                                                 >  1:117895/1‑150 (MQ=255)
                                    gatgatCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGcctc                                                                                                               <  2:1863512/150‑1 (MQ=255)
                                    gatgatCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGcctc                                                                                                               >  2:1218307/1‑150 (MQ=255)
                                    gatgatCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTGAAAAAGccaa                                                                                                                                             <  2:184882/120‑2 (MQ=255)
                                    gatgatCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTGAAAAAGccaa                                                                                                                                             >  1:184882/1‑119 (MQ=255)
                                       gatCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCtt                                                                                                            <  2:2003684/150‑1 (MQ=255)
                                         tCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTCTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTgc                                                                                                          >  2:2004434/1‑150 (MQ=255)
                                            gctgGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTgcgcg                                                                                                       >  1:590047/1‑150 (MQ=255)
                                                  tCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACt                                                                                                                                 <  2:355358/118‑1 (MQ=255)
                                                  tCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACt                                                                                                                                 >  1:355358/1‑118 (MQ=255)
                                                        gcgAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCtgtg                                                                                           >  2:276928/1‑150 (MQ=255)
                                                         cgAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGc                                                                                          <  2:626560/150‑1 (MQ=255)
                                                          gAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGcc                                                                                         <  1:206469/150‑1 (MQ=255)
                                                          gAGTCAGCCACGATAACCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGcc                                                                                         <  2:750760/150‑1 (MQ=255)
                                                               acacACGATGCCCAATGGAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATc                                                                                    >  2:225244/4‑150 (MQ=255)
                                                                gCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCACTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGcctc                                                                                                               >  2:2115114/1‑122 (MQ=255)
                                                                gCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCACTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGcctc                                                                                                               <  1:2115114/122‑1 (MQ=255)
                                                                  cACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCctc                                                                                 <  2:590047/150‑1 (MQ=255)
                                                                          ccAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGGGCGCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGc                                                                         >  1:1833567/1‑150 (MQ=255)
                                                                          ccAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCa                                                                                      <  1:1022804/137‑1 (MQ=255)
                                                                          ccAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCa                                                                                      >  2:1022804/1‑137 (MQ=255)
                                                                          ccAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGc                                                                         >  2:1878444/1‑150 (MQ=255)
                                                                          ccAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCCTGCGCGCCTCGTGCTGTGCCTCATCCTCTTcc                                                                             <  2:568075/146‑1 (MQ=255)
                                                                          ccAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCCTGCGCGCCTCGTGCTGTGCCTCATCCTCTTcc                                                                             >  1:568075/1‑146 (MQ=255)
                                                                           cAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCt                                                                        >  1:619516/1‑150 (MQ=255)
                                                                           cAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCGGTTGCTCATCCTCTTCCTTGCt                                                                        >  2:561148/1‑150 (MQ=255)
                                                                              ggTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCtgtg                                                                     >  1:2293058/1‑150 (MQ=255)
                                                                                tAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCATTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCa                                                                   >  1:1166689/1‑150 (MQ=255)
                                                                                 aaTGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCaa                                                                  <  1:1796188/150‑1 (MQ=255)
                                                                                   tGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGa                                                                >  2:1665685/1‑150 (MQ=255)
                                                                                     gCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCTTTGCTGTGCAAGAAc                                                              <  1:225244/150‑1 (MQ=255)
                                                                                     gCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAAc                                                              >  2:2263453/1‑150 (MQ=255)
                                                                                      cGATGGTTTCTGCCGACAGGATGAAGTCGGTACGCAGCGCCCCCTTGATCTTGGCTTTTTCAAAGGCCAACAGGTCAACACTGGCATCGGGGACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACt                                                             >  2:461021/1‑150 (MQ=255)
                                                                                       gatgatTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACtt                                                            <  1:2207262/150‑1 (MQ=255)
                                                                                         tgatTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTtgt                                                          <  1:1881534/150‑1 (MQ=255)
                                                                                         ggatTTCTGCCCACAGGATGATCTAGGTGAGCCCCGCCCTGCTTATCGGGTCTTTTGCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGTGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTtgt                                                          <  2:1833567/149‑1 (MQ=255)
                                                                                           atTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTgc                                                        <  2:1098741/150‑1 (MQ=255)
                                                                                                   cGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCATATTTCTTCTTGTGCAAGAACTTGTGCGCCAGCtt                                                >  2:1508407/1‑150 (MQ=255)
                                                                                                    gACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCCAGACCTTGTGAGCCCGCTtc                                               >  2:2315979/1‑150 (MQ=255)
                                                                                                      cAGGGTGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTtctc                                             <  2:763932/150‑1 (MQ=255)
                                                                                                       aGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCg                                            <  1:1218307/150‑1 (MQ=255)
                                                                                                         gATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAAc                                                              <  2:1699595/130‑1 (MQ=255)
                                                                                                         gATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAAc                                                              >  1:1699595/1‑130 (MQ=255)
                                                                                                         gATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTtctc                                             <  2:742225/147‑1 (MQ=255)
                                                                                                         gATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTtctc                                             >  1:742225/1‑147 (MQ=255)
                                                                                                          aTGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCGaaa                                         <  2:2119646/150‑1 (MQ=255)
                                                                                                          aTGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCGaaa                                         <  1:276928/150‑1 (MQ=255)
                                                                                                          aTGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCGCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCGaaa                                         >  2:1646332/1‑150 (MQ=255)
                                                                                                              aGTCGGTACGCACCGCCCCCTTGATCTTGGCTTCTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCGAAACCCt                                     <  1:2004434/150‑1 (MQ=255)
                                                                                                               gTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCGGCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCGAAACCCTc                                    >  1:263001/1‑150 (MQ=255)
                                                                                                                  ggTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCGAAACCCTCGaa                                 <  1:2263453/150‑1 (MQ=255)
                                                                                                                  ggTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCGAAACCCTCGaa                                 <  2:619516/150‑1 (MQ=255)
                                                                                                                    tACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCGAAACCCTCGAAGc                               >  1:930544/1‑150 (MQ=255)
                                                                                                                       gCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCGAAACCCTCGAAGcaca                            >  1:964254/1‑150 (MQ=255)
                                                                                                                       gCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCGAAACCCTCGAAGcaca                            <  2:947649/150‑1 (MQ=255)
                                                                                                                        cACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCTACCAGGTCAACACTGGCATGCGCAACGGCCTCCTTGCGCGCCTCGTTCTGTGCGTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGATTTTCGAACCCCTCAACGcacta                           >  2:337584/1‑148 (MQ=255)
                                                                                                                         aCCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTcc                                                                             >  2:750287/1‑99 (MQ=255)
                                                                                                                         aCCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTcc                                                                             <  1:750287/99‑1 (MQ=255)
                                                                                                                         aCCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCtgt                                                                      >  2:1769059/1‑106 (MQ=255)
                                                                                                                         aCCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCtgt                                                                      <  1:1769059/106‑1 (MQ=255)
                                                                                                                             cccccTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCGAAACCCTCGAAGCACAGGAAcg                      >  2:1105356/1‑150 (MQ=255)
                                                                                                                                ccTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCtt                                                                                                            >  2:720399/1‑61 (MQ=255)
                                                                                                                                ccTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCtt                                                                                                            <  1:720399/61‑1 (MQ=255)
                                                                                                                                    gATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCa                                                                                      >  2:53456/1‑79 (MQ=255)
                                                                                                                                    gATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCa                                                                                      <  1:53456/79‑1 (MQ=255)
                                                                                                                                      tggtGGCTTTTGAAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCtg                                                                                             <  1:712444/67‑1 (MQ=38)
                                                                                                                                      tggtGGCTTTTGAAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCtg                                                                                             >  2:712444/4‑70 (MQ=38)
                                                                                                                                      tCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCACCGTCTCGAAACCCTCGAAGCACAGGAACGCGCCACCCa             <  1:1508407/150‑1 (MQ=255)
                                                                                                                                        ttGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGc                                                                                          >  2:665656/1‑71 (MQ=255)
                                                                                                                                        ttGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGc                                                                                          <  1:665656/71‑1 (MQ=255)
                                                                                                                                        ttGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCGAAACCCTCGAAGCACAGGAACGCGCCACCCAGc           >  2:2206301/1‑150 (MQ=255)
                                                                                                                                         tGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCGAAACCCTCGAAGCACAGGAACGCGCCACCCAGCa          <  2:117895/150‑1 (MQ=255)
                                                                                                                                         tGGCTTTTGAAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCGAAACCCTCGAAGCACAGGAACGCGCCACCCAGCa          >  1:448892/1‑150 (MQ=255)
                                                                                                                                             tttttCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGGGGCTCGTGCTTTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCGAAACCCTCGAAGCACAGGAACGCGCCACCCAGCCTcaa      >  2:586008/1‑150 (MQ=255)
                                                                                                                                              ttttCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCGAAACCCTCGAAGCACAGGAACGCGCCACCCAGCATcaac     <  1:1878444/150‑1 (MQ=255)
                                                                                                                                              tcttCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCGAAACCCTCGAAGCACAGGAACGCGCCACCCAGCATcaac     <  1:1611775/148‑1 (MQ=255)
                                                                                                                                                ttCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCGAAACCCTCGAAGCACAGGAACGCGCCACCCAGCATcaacaa   <  1:2206301/150‑1 (MQ=255)
                                                                                                                                                 tCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCGAAACCCTCGAAGCACAGGAACGCGCCACCCAGCATCaacaac  <  1:2240588/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
ACCAATGGTCATCACCACCGCAATGCCGGAAAGCACGATGATCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCGAAACCCTCGAAGCACAGGAACGCGCCACCCAGCATCAACAAC  >  NC_002947/4884862‑4885156

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