Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 197,402 C→T T286I (ACA→ATA)  rseP → intramembrane zinc metalloprotease RseP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913197,4020CT90.6% 161.2 / 4.8 64T286I (ACA→ATA) rsePintramembrane zinc metalloprotease RseP
Reads supporting (aligned to +/- strand):  ref base C (2/4);  new base T (30/28);  total (32/32)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.40e-01

CCTGATAAAGAAGATCCGGTATCTTCGCTGGGGATTCGTCCTCGTGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGCAGTGGGTGACCTTTGTGATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGTCGGTATTGAGCCGAAAGTCATTCCTTTGCCAGATGAGTATAAAGTTGTACGCCAGTATGGGCCGTTCAACGCCATCGTCGAAGCCACGGACAAAACGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATCGCCA  >  NC_000913/197155‑197641
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                               
ccTGATAAAGATCATCCGGAATCATCGCTGGGGATTCGTCCTCTTGGGCCACAAATTGAACCTGTAATGGAAAATGTGCAGCCACACTCGGCGGCATGCAAGGCCGGTTTGCAAGCAGGCGACAGGCTCGTTAAAGTCGAGGGTCAGCCCTTAACGCAGTGGGTGACCTTTGTGATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGCAGTGCCGGGGAGACCCTTGTCTTTGATAtt                                                                                                                                                                                                                                              >  2:539532/1‑251 (MQ=255)
 cTGATAAAGAAGATCCGGTATCTTCGCTGGGGATTCGTCCTCGTGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGCAGTGGGTGACCTTTGTGATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGATAtt                                                                                                                                                                                                                                              >  2:349075/1‑250 (MQ=255)
     aaaaGAAGATCCGGTATCTTCGCTGGGGATTCGTCCTCGTGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGCAGTGGGTGACCTTTGTGATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGATATTAATcc                                                                                                                                                                                                                                         <  2:148674/250‑1 (MQ=255)
       aagaagATCCGGTATCTTCGCTGGGGATTCGTCCTCGTGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGCAGTGGGTGACCTTTGTGATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGATATTAATcc                                                                                                                                                                                                                                         >  2:550935/1‑249 (MQ=255)
       aagaagATCCGGTATCTTCGCTGGGGATTCGTCCTCGTGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGCAGTGGGTGACCTTTGTGATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGATATTAATcc                                                                                                                                                                                                                                         <  1:550935/249‑1 (MQ=255)
          aataTCCGGTATCTCCGCGTGGGATTCGTCCTCGTGGGTCGCTAAGTGAATCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCATGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGCAGTGGGTGACCTTTGTGTTGCTTGGCCTGGATATCCCTGGTAAACCCTTAGCGTTAGAAATCGAAAGGCAGTGTAGTCCCTTGTCTTGGATCTTAGTCCCGgag                                                                                                                                                                                                                                    <  2:74079/248‑1 (MQ=255)
              tCCGGTATCTTCGCTGGGGATTCGTCCTCGTGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGCAGTGGGTGACCTTTGTGATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGATATTAATCCCGGAGAGt                                                                                                                                                                                                                                 >  1:391227/1‑250 (MQ=255)
                 ggtaacgtCGCTGGGGATTCGTCCTCGTGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGCAGTGGGTGACCTTTGTGATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTATAAAGCGAAAGGCAGGGGAGGCCCTTGTCTTTGATATTAATTCCGGAGAGTAAAc                                                                                                                                                                                                                             <  2:599766/244‑1 (MQ=255)
                 ggTATCTTCGCTGGGGATTCGTCCTCGTGGGCCGCAACTTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGCAGTGGGTGACCTTTGTGTTGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGATATTAATCCCGGAGAGAAAAc                                                                                                                                                                                                                             >  2:4597/1‑251 (MQ=255)
                   tATCTTCGCTGGGGATTCGTCCTCGTGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGCAGTGGGTGACCTTTGTGATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGATATTAATCCCGGAGAGTAAACCg                                                                                                                                                                                                                           >  2:63712/1‑251 (MQ=255)
                           cTGGGGATTCGTCCTCGTGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGCAGTGGGTGACCTTTGTGATGCTTGTCCGGGCTAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGATATTAATCCCGGGGAGTAAACCGGGTAATgg                                                                                                                                                                                                                   >  2:346586/1‑251 (MQ=255)
                                    cGTCCTCGTGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGCAGTGGGTGACCTTTGTGATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGATATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGa                                                                                                                                                                                                           <  1:150486/250‑1 (MQ=255)
                                      cccTCGTGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGCAGTGGGTGACCTTTGTGATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGATATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTg                                                                                                                                                                                                        <  2:228594/250‑1 (MQ=255)
                                            ggggCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGCAGTGGGTGACCTTTGTGATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGCGGAGTCCCTTGTCTTTGATATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTg                                                                                                                                                                                                  <  1:305812/250‑1 (MQ=255)
                                                ccGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGCAGTGGGTGACCTTTGTGATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGATATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGTCgg                                                                                                                                                                                              >  2:545964/1‑251 (MQ=255)
                                                  gCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGCAGTGGGTGACCTTTGTGATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGATATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGTCGGTa                                                                                                                                                                                            <  1:348958/251‑1 (MQ=255)
                                                   cAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGCAGTGGGTGACCTTTGTGATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGATATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGTCGGTAt                                                                                                                                                                                           <  1:365312/251‑1 (MQ=255)
                                                   cAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGCAGTGGGTGACCTTTGTGATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGATATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGTCGGTAt                                                                                                                                                                                           <  1:235696/251‑1 (MQ=255)
                                                     aaTTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGACAGCCCTTAACGCAGTGGGTGACCTTTGTGATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGATATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGTCGGTATTg                                                                                                                                                                                         >  1:257797/1‑251 (MQ=255)
                                                       ttGAACCTGTACTGGAAAATGTGCAGCCGAACTCGGCGGCAAGCACGGCAGGTTTGCAAGCAGGCGACAGGAACGTTAAAGTCGAGGGTCAACCCTTAAAAAAGTGGGTGAGCTTTGAGATGCTTATCCGGGATAAACCGGGTAAATACTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCACTAGACTTTGATATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGTCGGTATTg                                                                                                                                                                                         >  1:628773/1‑249 (MQ=255)
                                                        tGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGCAGTGGGTGACCTTTGTGATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGATATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGTCGGTATTGAg                                                                                                                                                                                       <  1:534543/250‑1 (MQ=255)
                                                               gTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGCAGTGGGTGACCTTTGTGATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGATATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGTCGGTATTGAGCCGAAAGt                                                                                                                                                                               <  1:299225/251‑1 (MQ=255)
                                                                              cAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGAGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGCAGTGGGTGACCTTTGTGATGATTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGATATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGTCGGTATTGAGCCGAAAGTCATTCCTTTGcc                                                                                                                                                                   >  2:264618/1‑248 (MQ=255)
                                                                                         ggcggcAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGCAGTGGGTGCCCTTTGTGATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGATATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGTCGGTATTGAGCCGAAAGTCATTCCTTTGCCAGATGAGTATAAAg                                                                                                                                                     <  1:392431/251‑1 (MQ=255)
                                                                                         ggcggcAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGCAGTGGGTGACCTTTGTGATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGATATTAAGCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGGCGGTATTGAGCCTAAAGTCATTCCTTTGCCAGATGAGTATAAAg                                                                                                                                                     <  1:290390/251‑1 (MQ=255)
                                                                                           cggcAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGCAGTGGGTGACCTTTGTGATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGATATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGTCGGTATTGAGCCGAAAGTCATTCCTTTGCCAGATGAGTATAAAGtt                                                                                                                                                   >  1:561563/1‑251 (MQ=255)
                                                                                            ggcAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGCAGTGGGTGACCTTTGTGATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGATATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGTCGGTATTGAGCCGAAAGTCATTCCGTTGCCAGATGAGTATAAAGTTg                                                                                                                                                  <  2:73586/251‑1 (MQ=255)
                                                                                                  caagGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGCAGTGGGTGACCTTTGTGATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGTCGGTATTGAGCCGAAAGTCATTCCTTTGCCAGATGAGTATAAAGTTGTACGcc                                                                                                                                            >  1:507607/1‑251 (MQ=255)
                                                                                                                 aGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGCAGTGGGTGACCTTTGTGATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGATATTAATCCCGGAGAGTAACCCGGGTAATGGTAAAGCGATTGGTTTTGTCGGTATTGAGCCGAAAGTCATTCCTTTGCCAGATGGGTATAAAGTTGTACGCCAGTATGGGCCGTTc                                                                                                                              >  1:269239/1‑250 (MQ=255)
                                                                                                                            ggATCGTTAAAGGCGATGGTCAGCCCTTAACGCAGTGGGTGACCTTTGTGATGCTTGGCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGATATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGTCGGTATTGAGCCGAAAGTCATTCCTTTGCCAGATGAGTATAAAGTTGTACGCCAGTATGGGCCGTTCAACGCCAtcgtc                                                                                                                  <  1:546204/251‑1 (MQ=255)
                                                                                                                                cGTTAAAGAAGATAGTCAGCCCTTCACGCAGTGGGTGACCTTTGCGATGCTTGTCCGGGATAACCCGGATAAATCCTTTGAGTTAGAACTCGAAAGGCAGGGGAGTCCCTTGTCTTTGATATTAATCCCGGAGAGTAAACCGGGTACTGGTAAAGCGATTGGTTTTGTCGGTATTGAGCCGAAAGTCATTCCTTTGACAGATGAGTATAAAGTTGTACGCCAGTATGGGCCGTTCAGCGCCATCGTCGaac                                                                                                              >  2:583685/1‑250 (MQ=255)
                                                                                                                                         cGATGGTCAGCCCTTAACGCAGTGGGTGACCTTTGTGATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGTCGGTATTGAGCCGAAAGTCATTCCTTTGCCAGATGAGTATAAAGTTGTACGCCAGTATGGGCCGTTCAACGCCATCGTCGAAGCCACGGACa                                                                                                     <  1:413675/251‑1 (MQ=255)
                                                                                                                                          gATGGTCAGCCCTTAACGCAGTGGGTGACCTTTGTGATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGATATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGTCGGTATTGAGCCGAAAGTCATTCCTTTGCCAGATGAGTATAAAGTTGTACGCCAGTATGGGCCGTTCAACGCCATCGTCGAAGCCACGGACaa                                                                                                    <  1:561038/251‑1 (MQ=255)
                                                                                                                                           aTGGTCAGCCCTTAACGCAGTGGGTGACCTTTGTGATGCTTGTCCGGGCTAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGATATTAATCCCGGAGAGTAACCCGGGTAATGGTAAAGCGATTGGTTTTGTCGGTATTGAGCCGAAAGTCCTTCCTTTGCCAGATGAGTATAAAGTTGTACGCCAGTATGGGCCGTTCAACGCCATCGTCGAAGCCACGGACaa                                                                                                    >  1:51484/1‑250 (MQ=255)
                                                                                                                                             ggTCAGCCCTTAACGCAGTGGGTGACCGTTGTGATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGATATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGTCGGTATTGAGCCGAAAGTCATTCCTTTGCCAGATGAGTATAAAGTTGTACGCCAGTATGGGCCGTTCAACGCCATCGTCGAAGCCACGGACAAAAc                                                                                                 <  1:517764/251‑1 (MQ=255)
                                                                                                                                                cAGCCCTTAACGCAGTGGGTGACCTTTGTGATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGATATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGTCGGTATTGAGCCGAAAGTCATTCCTTTGCCAGATGAGTATAAAGTTGTACGCCAGTATGGGCCGTTCAACGCCATCGTCGAAGCCACGGACAAAACGTg                                                                                              <  1:467293/251‑1 (MQ=255)
                                                                                                                                                cAGCCCTTAACGCAGTGGGTGACCTTTGTGATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGATATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGTCGGTATTGAGCCGAAAGTCATTCCTTTGCCAGATGAGTATAAAGTTGTACGCCAGTATGGGCCGTTCAACGCCATCGTCGAAGCCACGGACAAAACGTg                                                                                              <  2:387862/251‑1 (MQ=255)
                                                                                                                                                    ccTTAACGCAGTGGGTGACCTTTGTGATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGATATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGTCGGTATTGAGCCGAAAGTCATTCCTTTGCCAGATGAGTATAAAGTTGTACGCCAGTCTGGGCCGTTCAACGCCATCGTCGAAGCCACGGACAAAACGTGGc                                                                                            >  2:139889/1‑249 (MQ=255)
                                                                                                                                                     cTTAACGCAGTGGGTGACCTTTGTGATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGATATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGTCGGTATTGAGCCGAAAGTCATTCCTTTGCCAGATGAGTATAAAGTTGTACGCCAGTATGGGCCGTTCAACGCCATCGTCGAAGCCACGGACAAAACGTGGCAGc                                                                                         >  1:292066/1‑251 (MQ=255)
                                                                                                                                                         aCGCAGTGGGTGACCTTTGTGATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGATATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGTCGGTATTGAGCCGAAAGTCATTCCTTTGCCAGATGAGTATAAAGTTGTACGCCAGTATGGGCCGTTCAACGCCATCGTCGAAGCCACGGACAAAACGTGGCAGCtgat                                                                                     >  1:332848/1‑251 (MQ=255)
                                                                                                                                                           gCAGTGGGTGACCTTTGTGATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGATATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGTCGGTATTGAGCCGAAAGTCATTCCTTTGCCAGATGAGTATAAAGTTGTACGCCAGTATGGGCCGTTCAACGCCATCGTCGAAGCCACGGACAAAACGTGGCAGCtgatga                                                                                   <  1:85851/251‑1 (MQ=255)
                                                                                                                                                            cAGTGGGTGACCTTTGTGATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGATATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGTCGGTATTGAGCCGAAAGTCATTCCTTTGCCAGATGAGTATAAAGTTGTACGCCAGTATGGGCCGTTCAACGCCATCGTCGAAGCCACGGACAACACGTGGCAGCtgatga                                                                                   >  2:612227/1‑250 (MQ=255)
                                                                                                                                                             aGTGGGTGACCTTTGTGTTGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGATATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGCTTGGTTTTGTCGGTATTGAGCCGAAAGTCATTCCTTTGCCAGATGAGTATAAAGATGTACGCCAGTATGGGCCGTTCAACGCCATCGTCGAAGCAACGGACAAAACGTGGCAGCtgatga                                                                                   >  2:160025/1‑249 (MQ=255)
                                                                                                                                                                gggTGGCCTTTGTGATGCGTGTCCGTGATAACCCGGGTAAATCCTTAGCGGTAGAAATCGAAAGGCAGGGGAGGCCCTTGTCTTTGATAGTAATCCCGGAGAGTAGACCGGGTAATGGTAAAGCGATTGGTTTTGTCGGTATTGAGCTGAAAGACATTCCGTTGCCAGAGGAGTATAAAGTTGTACGCCAGTATGGGGCGTTCAGCTCCATCGTCGAAGCCACGGACAAAACGTGGCAGCTGATGAAGCt                                                                               <  2:205844/250‑1 (MQ=255)
                                                                                                                                                                gggTGACCTTTGTGATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGTCGGTATTGAGCCGAAAGTCATTCCTTTGCCAGATGAGTATAAAGTTGTACGCCAGTATGGGCCGTTCAACGCCATCGTCGAAGCCACGGACAAAACGTGGCAGCTGATGAAGCTg                                                                              <  1:480089/251‑1 (MQ=255)
                                                                                                                                                                             gATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGATATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGTCGGTATTGAGCCGAAAGTCATTCCTTTGCCAGATGAGTATAAAGTTGTACGCCAGTATGGGCCGTTCAACGCCATCGTCGAAGCCACGGACAAAACGTGGCAGCTGATGAAGCTGACGGTCAGTATg                                                                  <  1:545964/250‑1 (MQ=255)
                                                                                                                                                                                    gTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGATATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGTCGGTATTGAGCCGAAAGTCATTCCTTTGCCAGATGAGTATAAAGTTGTACGCCAGTATGGGCCGTTCAACGCCATCGTCGAAGCCACGGACAAAACGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGaaa                                                          <  2:154262/251‑1 (MQ=255)
                                                                                                                                                                                    gTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGATATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGTCGGTATTGAGCCGAAAGTCATTCCTTTGCCAGATGAGTATAAAGTTGTACGCCAGTATGGGCCGTTCAACGCCATCGTCGAAGCCACGGACAAAACGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGaa                                                           <  1:117256/250‑1 (MQ=255)
                                                                                                                                                                                        gggATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGATATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGTCGGTATTGAGCCGAAAGTCATTCCTTTGCCAGATGAGTATAAAGTTGTACGCCAGTATGGGCCGTTCAACGCCATCGTCGAAGCCACGGACAAAACGTGGCAGCTGATGAAGCTGAGGTTCAGTATGCTGGGAAAATTg                                                      <  1:512524/251‑1 (MQ=255)
                                                                                                                                                                                        gggATAACCCGGGTAAATCCTTAGCGTTAGAAAGCGAAAGGCAGGGGAGTCCCTTGTCTTTGATATTAATCCCGGAGAGTAAACCGGTTAATGGTAAAGCGATTGGTTTTGTCGGTATTGAGCCGAAAGTCATTCCTTTGCCAGATGAGTATAAAGTTGTACGCCAGTATGGGCCGTTCAACGCCATCGTCGAAGCCACGGACAAAACGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTg                                                      <  1:393988/251‑1 (MQ=255)
                                                                                                                                                                                           aTAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGTCGGTATTGAGCCGAAAGTCATTCCTTTGCCAGATGAGTATAAAGTTGTACGCCAGTATGGGCCGTTCAACGCCATCGTCGAAGCCACGGACAAAACGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATc                                                   <  1:478963/251‑1 (MQ=255)
                                                                                                                                                                                                 cGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGATATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGTCGGTATTGAGCCGAAAGTCATTCCTTTGCCAGATGAGTATAAAGTTGTACGCCAGTATGGGCCGTTCAACGCCATCGTCGAAGCCACGGACAAAACGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCggg                                             >  2:95675/1‑250 (MQ=255)
                                                                                                                                                                                                   ggTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGATATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGTCGGTATTGAGCCGAAAGTCATTCCTTTGCCAGATGAGTATAAAGTTGTACGCCAGTATGGGCCGTTCAACGCCATCGTCGAAGCCACGGACAAAACGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTg                                            >  2:439945/1‑250 (MQ=255)
                                                                                                                                                                                                    gTAAATCCTCAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGATATTAATCCCGGTGAGTAACCCGGGTACTGGTAAAGCGATTGGTTTTGTCGGTATTGAGCCGAAAGCCATTCCTTTGCCAGATGCGTATAAAGTTGTACGCCAGTATGGGCCGTTCAACGCCATCGTAGAAGACACGGAAAAAACGTGGCAGCTGATGACGCTGACGGTCAGAATGCTGGGAAAATTGATCACCGGTGAt                                          >  2:557030/1‑251 (MQ=255)
                                                                                                                                                                                                            ttAGCGTTAGAACCCGAAAGGCAGGGGGGTCCCTTGTCTTTGATATTAATCCCGGCGAGAAAACCGGGTAATGGTAAAGCGATTGGTTTTGTCGGTATTGAGCCGAACGTCATTCCTTTGCCAGATGCGTATACAGTTGTACGCCAGTATGGGCCGTTCAACGCCATCGTCGCAGCCACGGACAAAACGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGGGAAACt                                  >  2:308076/1‑251 (MQ=255)
                                                                                                                                                                                                                     gAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGATATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGTCGGTATTGAGCCGAAAGTCATTCCTTTGCCAGATGAGTATAAAGTTGTACGCCAGTATGGGCCGTTCAACGCCATCGTCGAAGACACGGACAAAACGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCt                         >  1:243895/1‑251 (MQ=255)
                                                                                                                                                                                                                      aaagcGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGTCGGTATTGAGCCGAAAGTCATTCCTTTGCCAGATGAGTATAAAGTTGTACGCCAGTATGGGCCGTTCAACGCCATCGTCGAAGCCACGGACAAAACGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTc                        <  2:209686/247‑1 (MQ=255)
                                                                                                                                                                                                                           gAAAGGCAGGGGAGTCCCTTGTCTTTGATATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGTCGGTATTGAGCCGAAAGTCATTCCTTTGCCAGATGAGTATAAAGTTGTACGCCAGTATGGGCCGTTCAACGCCATCGACGAAGCCACGGACAAAACGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGAGATGTGAAACTGAACAACCTCAGTgg                   >  1:236275/1‑251 (MQ=255)
                                                                                                                                                                                                                            aaaGGCAGGGGAGTACCTTGTCTTTGATATTAATCCCGGAGAGTAAACAGGGTAATGGTAAAGCGATTGGTTTTGTCGGTATTGAGCCGAAAGTCATTCCTTTGCCAGATGAGTATAAAGTTGTACGCCAGTATGGGCCGTTCAACGCCATCGTCGAAGCCACGGACAAAACGTGGCAGCTGATGAAGCTGACGGTAAGTATGCTGGGAAAATTGATCACAGGTGATGTGAAACTGAACAACCTCAGTggg                  >  1:480923/1‑251 (MQ=255)
                                                                                                                                                                                                                                   ggggAGTCCCTTGTCTTTGATATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGTCGGTATTGAGCCGAAAGTCATTCCTTTGCCAGATGAGTATAAAGTTGTACGCCAGTATGGGCCGTTCAACGCCATCGTCGAAGCCACGGACAAAACGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGAtct           >  1:627029/1‑251 (MQ=255)
                                                                                                                                                                                                                                    gggAGTCCCTTGTCTTTGATATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGTCGGTATTGAGCCGAAAGTCATTCCTTTGCCAGATGAGTATAAAGTTGTACGCCAGTATGGGCCGTTCAACGCCATCGTCGAAGCCACGGACAAAACGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGAtctc          >  2:117962/1‑251 (MQ=255)
                                                                                                                                                                                                                                       aGTCCCTTGTCTTTGACATTAATAACGGAGAGTAAAACGGGTAATGGTAAAGCGATAGGTTTTGTCGGTATTGAGCCGAAAGTCATTCCTTTGCCAGATGAGTTTAAAGTTGTACGCCAGTATCGGCCGTTCAACGCCATAGTCGAAGCCACGGACAAAACGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAACATGATCACCGGCGATGTGAAACTCAACAACCTCAGTGGGCCGATCTCTAt       >  2:392354/1‑251 (MQ=255)
                                                                                                                                                                                                                                        gTCCCTTGTCTTTGATATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGTCGGTATTGAGCCGAAAGTCATTCCTTTGCCAGATGAGGATAAAGTTGTACGCCAGTATGGGCCGTTCAACGCCATCGTCGAAGCCACGGACAAAACGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCGCTATc      <  1:529506/251‑1 (MQ=255)
                                                                                                                                                                                                                                            cTTGTCTTTGATATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGTCGGTATTGAGCCGAAAGTCATTCCTTTGCCAGATGAGTATAAAGTTGTACGCCAGTATGGGCCGTTCAACGCCATCGTCGAAGCCACGGACAAAACGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCCGTGGGCCGATCTCTATCGccc  >  2:129830/1‑250 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                               
CCTGATAAAGAAGATCCGGTATCTTCGCTGGGGATTCGTCCTCGTGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGCAGTGGGTGACCTTTGTGATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGTCGGTATTGAGCCGAAAGTCATTCCTTTGCCAGATGAGTATAAAGTTGTACGCCAGTATGGGCCGTTCAACGCCATCGTCGAAGCCACGGACAAAACGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATCGCCA  >  NC_000913/197155‑197641

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: