Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,057,150 C→T 6.1% intergenic (+162/+28) sibC → / ← serA sRNA antisense regulator of toxic IbsC protein/D‑3‑phosphoglycerate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,057,1500CT 6.1% 423.4 / 14.2 133intergenic (+162/+28)sibC/serAsRNA antisense regulator of toxic IbsC protein/D‑3‑phosphoglycerate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (60/62);  new base T (4/4);  total (66/67)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.25e-02

Only 100 of 136 total aligned reads displayed.

GAGACTGACTGTTAATAAGCGCTGAAACTTATGAGTAACAGTACAATCAGTATGATGACAAGTCGCATCATAACCCTTCTCCTTCAAGCCCTCGCTTCGGTGAGGGCTTTACCGTTACAGCCCCATGCTGCCCTGCCATCGTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCAGCGCAGTTAGCACGCCCGGACGGTTTTCGTGGATGTGCATCAGACGACGCCCACCGT  >  NC_000913/3056879‑3057407
                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                 
gagaCTGACTGTTAATAAGCGCTGAAACTTATGAGTAACAGTACAATCAGTATGATGACAAGTCGCATCATAACCCTTCTCCTTCAAGCCCTCGCTTCGGTGAGGGCTTTACCGTTACAGCCCCATGCTGCCCTGCCATCGTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGc                                                                                                                                                                                                                                                                >  4:303954/1‑275 (MQ=255)
      gactgTTAATAAGCGCTGAAACTTATGAGTAACAGTACAATCAGTATGATGACAAGTCGCATCATAACCCTTCTCCTTCAAGCCCTCGCTTCGGTGAGGGCTTTACCGTTACAGCCCCATGCTGCCCTGCCATCGTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTg                                                                                                                                                                                                                                                          <  2:303955/275‑1 (MQ=255)
             aataaGCGCTGAAACTTATGAGTAACAGTACAATCAGTATGATGACAAGTCGCATCATAACCCTTCTCCTTCAAGCCCTCGCTTCGGTGAGGGCTTTACCGTTACAGCCCCATGCTGCCCTGCCATCGTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCa                                                                                                                                                                                                                                                   >  1:68328/1‑275 (MQ=255)
                aaaagagaaaaCTTATGAGTAACATTACAATCAGGAGGATGACAAGTCGCATCATAACCCTTCTCCTTCAAGCCCTCGCTTCGGTGAGGGCTTTACCGTTACAGCCCCATGCTGCCCTGCCATCGTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCa                                                                                                                                                                                                                                                   <  4:127627/264‑1 (MQ=255)
                    atcaaaaCTTAATAGTAAAATAACAATCANNTTGATAACAAGTGCAACCACACCCCTACTCNTTCAAGCCCCCGCTTCGGTGAGGCCTTTACCGTTACAGCCCCATGCTGCCCTGCCATCGTAAATCCCCATTAAANAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAgg                                                                                                                                                                                                                                            <  3:68328/271‑1 (MQ=255)
                          aCTTATGAGTAACAGTAAAATCAGTATGATGACAAGTCGCCTCACAACCCTTCTCCTTCAAGCCCTCGATTCGGTGAGGGCTTTACCGTTACAGCCCCATGCTGCCCTGCCATCGTAAAGCCCCATTAAATAAACACACCGCATTGTTCTGACTTTGATTGATTTTCTGGAGCAGACTCGCAAAGTAGCATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAAt                                                                                                                                                                                                                                       <  1:128988/274‑1 (MQ=255)
                                      ccgTACAGACAGTATGATGACAAGTCGCATCAAAACCCTTCTCCTTCAAGCCCTCGCTTCGGTGAGGGCTTTACGGTTACAGCCCCATGCCGCCCCGCCATCGTAAATCCCCATTAAATAAACACAACGCGTTGACCCGACTGTGATTTATTTTACGGTGCAGACTCTCAAAGTAGCAGACGCAACTGGGCGACGGTGTGGAGAAGGGAAAAAAAAAAGGGCCAGCCAGCGACCTACCCGTCGATTTTCCCAGACGGGGAAGTCGCAcagaaaa                                                                                                                                                                                                                           >  4:145753/3‑270 (MQ=255)
                                      cAGTACAATCAGTATGATGACAAGTCGCATCATAACCCTTCTCCTTCAAGCCCTCGCTTCGGTGAGGGCTTTACCGTTACAGCCCCATGCTGCCCTGCCATCGTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGAtt                                                                                                                                                                                                                                                       >  2:127628/1‑246 (MQ=255)
                                            aaTCAGTATGATGACAAGTCGCATCATAACCCTTCTCCTTCAAGCCCTCGCTTCGGTGAGGGCTTTACCGTTACAGCCCCATGCTGCCCTGCCATCGTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGgcgc                                                                                                                                                                                                                    <  2:73873/275‑1 (MQ=255)
                                                 gTATGATGACAAGTCGCATCATAACCCTTCTCCTTCAAGCCCTCGCTTCGGTGAGGGCTTTACCGTTACAGCCCCATGCTGCCCTGCCATCGTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAGGGGATAAAAAAACGGGCAAGTCAGTGACATGCCCGTTGATTTTCAGAGAAGGGGAATTTATACCGCCGCCGGGCCCGAATg                                                                                                                                                                                                               >  4:233257/1‑275 (MQ=255)
                                                 gTATGATGACAAGTCGCATCATAACCCTTCTCCTTCAAGCCCTCGCTTCGGTGAGGGCTTTACCGTTACAGCCCCATGCTGCCCTGCCATCGTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAAAAGGGGAATTAGTACAGCAGACGGGCGCGact                                                                                                                                                                                                                >  2:439259/1‑272 (MQ=255)
                                                    ttaatccAAGTCGCATCCTAACCCTTCTCCTTCAAGCCCGCGCTTCGGTGAGGGCTTTACCGTTACAGCCCCATGCTGCCCTGCCATCGTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTa                                                                                                                                                                                                            <  4:311661/269‑1 (MQ=255)
                                                            aGTCGCATCATAACCCTTCTCCTTCAAGCCCTCGCTTCGGTGAGGGCTTTACCGTTACAGCCCCATGCTGCCCTGCCATCGTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCCAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCagaa                                                                                                                                                                                                                                                >  4:264462/1‑230 (MQ=255)
                                                                gCATCATAACCCTTCTCCTTCAAGCCCTCGCTTCGGTGAGGGCTTTACCGTTACAGCCCCATGCTGCCCTGCCATCGTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATg                                                                                                                                                                                                               <  4:233482/260‑1 (MQ=255)
                                                                              ctcCTTCAAGCCCTCGCTTCGGTGAGGGCTTTACCGTTACAGCCCCATGCTGCCCTGCCATCGTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGGATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCa                                                                                                                                                                                   >  2:461254/1‑274 (MQ=255)
                                                                                  ttCAAGCCCTCGCTTCGGTGAGGGCTTTACCGTTACAGCCCCATGCTGCCCTGCCATCGTAAATCCCCATTAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGcc                                                                                                                                                                                                                                                               <  2:121219/193‑1 (MQ=255)
                                                                                                             tACCGTGACAGCCCCATGCAGCCCTGACATCGTAAATCCCCATGAAATAAACACAACGCACTTATCTTACTTTAATTTATTTTCTAGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAGCGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTtca                                                                                                                                                   <  1:118366/275‑1 (MQ=255)
                                                                                                             accgtaccaccccccATCTGCCCTGCCATCTTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTTATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTtcaa                                                                                                                                                  <  4:461253/264‑1 (MQ=255)
                                                                                                              aCCGTTACAGCCCCATGCTGCCCTGCCATCGTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTtcaa                                                                                                                                                  <  2:233258/275‑1 (MQ=255)
                                                                                                                  ttACAGCCCCATGCTGCCCTGCCATCGTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCagag                                                                                                                                                                                                                                                >  2:207677/1‑177 (MQ=255)
                                                                                                                      aGCCCCATGCTGCCCTGCCATCGTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAgg                                                                                                                                                                                                                                            <  4:207676/177‑1 (MQ=255)
                                                                                                                       gCCCCATGCTGCCCTGCCATCGTAAATCCCCAATAAATAAACACAACGCATTGATCTGATTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTtcaatatcaata                                                                                                                                          <  1:56259/274‑1 (MQ=255)
                                                                                                                       gCCCCATGCTGCCCTGCCATCATAAATCCCCCTTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTtcaatatcaata                                                                                                                                          <  2:278178/274‑1 (MQ=255)
                                                                                                                        ccccATGCTGCCCTGCCATCGTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATaa                                                                                                                                         <  2:166252/274‑1 (MQ=255)
                                                                                                                          ccATGCTGCCCTGCCATCGTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCagag                                                                                                                                                                                                                                                >  4:375027/1‑169 (MQ=255)
                                                                                                                                    ctgccATCGTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGtt                                                                                                                                                                                                                                                           >  8:490922/1‑148 (MQ=255)
                                                                                                                                    cggccAACGTAAATCCCAACTAAATAAACAAAACGCATTGATCTGGCTTTGATTCATTTTCTGGGGCAGACTCCCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGGGACCTGCCCGTTGAGTTTCAGGGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCAt                                                                                                                            <  4:527336/273‑1 (MQ=255)
                                                                                                                                         aTCGTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGAttt                                                                                                                                                                                                                                                      >  8:700/1‑148 (MQ=255)
                                                                                                                                           cGTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTc                                                                                                                                                                                                                                                    <  7:490922/148‑1 (MQ=255)
                                                                                                                                           cGTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTc                                                                                                                                                                                                                                                    <  7:18176/148‑1 (MQ=255)
                                                                                                                                                   cccATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAgg                                                                                                                                                                                                                                            >  7:501897/1‑148 (MQ=255)
                                                                                                                                                      aTTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTACAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGACTGCAGCGCTTTTTCGGACACGTCTTCGTCGGCATCAATATCAATAACCCCATAACCCATCTGAGCGAAAGTTTGccg                                                                                                          >  4:124515/1‑273 (MQ=255)
                                                                                                                                                       ttAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGAt                                                                                                                                                                                                                                                        >  7:330033/1‑132 (MQ=255)
                                                                                                                                                          aataaaCACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTCCaga                                                                                                                                                                                                                                                 <  2:247566/136‑1 (MQ=255)
                                                                                                                                                              aaCACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTAc                                                                                                                                                                                                                                 <  7:700/148‑1 (MQ=255)
                                                                                                                                                                    aCGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTAca                                                                                                                                                                                                                                >  6:222710/1‑143 (MQ=255)
                                                                                                                                                                    aCGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTAca                                                                                                                                                                                                                                <  5:222710/143‑1 (MQ=255)
                                                                                                                                                                    aCGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGa                                                                                                                                                                                                                           <  8:501897/148‑1 (MQ=255)
                                                                                                                                                                       cATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGa                                                                                                                                                                                                                                                         >  7:968352/1‑115 (MQ=255)
                                                                                                                                                                       cATTGATCTGACTTTGATTCATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGa                                                                                                                                                                                                                                                         <  8:968352/115‑1 (MQ=255)
                                                                                                                                                                         ttGATCTGGCTTTAATTTGTTTTCTTGATCAGACTCGAAAAGTGGAAAGCTCAACGCGGGAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGGGACCTGCCCGTTGAGTTTCAGAGAAGGGGAATGAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCAGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACACATCTGGGCGGAAGTGTGCAGATATTGCGCGGCGATGTTg                                                                                       <  1:112823/275‑1 (MQ=255)
                                                                                                                                                                         ttGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAAtt                                                                                                                                                                                                                                      >  8:800180/1‑132 (MQ=255)
                                                                                                                                                                         ttGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAAtt                                                                                                                                                                                                                                      <  7:800180/132‑1 (MQ=255)
                                                                                                                                                                          tGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCCAAGTAGAATGTGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTTAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCGTTTTCGGCAACGTCTTCGTCGTCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGCTATGGCGGGGCGATGtttg                                                                                      >  4:109733/1‑273 (MQ=255)
                                                                                                                                                                             tCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGAcc                                                                                                                                                                                                                                                                   <  8:254212/99‑1 (MQ=255)
                                                                                                                                                                             tCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGAcc                                                                                                                                                                                                                                                                   >  7:254212/1‑99 (MQ=255)
                                                                                                                                                                                 aCTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTg                                                                               >  2:327216/1‑275 (MQ=255)
                                                                                                                                                                                    ttGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGaa                                                                                                                                                                                                                                        <  8:287645/119‑1 (MQ=255)
                                                                                                                                                                                    ttGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGaa                                                                                                                                                                                                                                        >  7:287645/1‑119 (MQ=255)
                                                                                                                                                                                         ttatttTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAGGTCACTGACTTGcc                                                                                                                                                                                                                                                               <  5:580022/91‑1 (MQ=255)
                                                                                                                                                                                         ttatttTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAGGTCACTGACTTGcc                                                                                                                                                                                                                                                               >  6:580022/1‑91 (MQ=255)
                                                                                                                                                                                         ttatttTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGcc                                                                                                                                                                                                                                                               >  5:448602/1‑91 (MQ=255)
                                                                                                                                                                                         ttatttTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGcc                                                                                                                                                                                                                                                               <  6:448602/91‑1 (MQ=255)
                                                                                                                                                                                         ttatttTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAg                                                                                                                                                                                                                                             <  5:603875/109‑1 (MQ=255)
                                                                                                                                                                                         ttatttTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAg                                                                                                                                                                                                                                             >  6:603875/1‑109 (MQ=255)
                                                                                                                                                                                               tCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGa                                                                                                                                                                                                                                                         >  5:827335/1‑91 (MQ=255)
                                                                                                                                                                                               tCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGa                                                                                                                                                                                                                                                         <  6:827335/91‑1 (MQ=255)
                                                                                                                                                                                                   gAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGaa                                                                                                                                                                                                                                        >  7:413580/1‑104 (MQ=255)
                                                                                                                                                                                                   gAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGaa                                                                                                                                                                                                                                        <  8:413580/104‑1 (MQ=255)
                                                                                                                                                                                                         aCTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCCGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGATTCAATATCAATAACCCCCTAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGAGATGTTGACGCCCTGCTCGGCGAAGCTATTGTTCAgcgc                                                       >  2:430001/1‑275 (MQ=255)
                                                                                                                                                                                                         aCTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGt                                                                                                                                                                                                                                   >  5:732948/1‑103 (MQ=255)
                                                                                                                                                                                                         aCTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGt                                                                                                                                                                                                                                   <  6:732948/103‑1 (MQ=255)
                                                                                                                                                                                                         aCTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGaa                                                                                                                                                                                                                 <  7:178972/121‑1 (MQ=255)
                                                                                                                                                                                                         aCTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGaa                                                                                                                                                                                                                 >  8:178972/1‑121 (MQ=255)
                                                                                                                                                                                                          cTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACg                                                                                                                                                                                                                         <  5:102218/112‑1 (MQ=255)
                                                                                                                                                                                                          cTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACg                                                                                                                                                                                                                         >  6:102218/1‑112 (MQ=255)
                                                                                                                                                                                                            cGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAGGTCACTGACTTGc                                                                                                                                                                                                                                                                <  6:503579/71‑1 (MQ=38)
                                                                                                                                                                                                            cGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAGGTCACTGACTTGc                                                                                                                                                                                                                                                                >  5:503579/1‑71 (MQ=38)
                                                                                                                                                                                                            cGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAg                                                                                                                                                                                                                                             >  8:441273/1‑90 (MQ=255)
                                                                                                                                                                                                            cGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAg                                                                                                                                                                                                                                             <  7:994839/90‑1 (MQ=255)
                                                                                                                                                                                                            cGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAg                                                                                                                                                                                                                                             <  7:441273/90‑1 (MQ=255)
                                                                                                                                                                                                            cGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTa                                                                                                                                                                                                                                     >  6:533622/1‑98 (MQ=255)
                                                                                                                                                                                                            cGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTa                                                                                                                                                                                                                                     <  5:533622/98‑1 (MQ=255)
                                                                                                                                                                                                                aaGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGt                                                                                                                                                                                                                                   >  6:618471/1‑96 (MQ=255)
                                                                                                                                                                                                                aaGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGt                                                                                                                                                                                                                                   <  5:618471/96‑1 (MQ=255)
                                                                                                                                                                                                                        tGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGAt                                                                                                                                                                                                                                                        >  8:1107523/1‑67 (MQ=255)
                                                                                                                                                                                                                        tGCGCAACGCGGCAACGGGGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGACCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACTGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTGGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCAGCGCAGTTAGCACGCCCg                                         <  2:109734/274‑1 (MQ=255)
                                                                                                                                                                                                                           gcAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTACAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATTACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCCAAGTTGGCGCCCTGCTCGGCGGCGCTTGTGTTCCGCCAAGTTTGCACTCCGGGCCg                                     >  4:487197/1‑275 (MQ=255)
                                                                                                                                                                                                                            cAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAgggg                                                                                                                                                                                                                                          <  6:444440/77‑1 (MQ=255)
                                                                                                                                                                                                                            cAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAgggg                                                                                                                                                                                                                                          >  5:444440/1‑77 (MQ=255)
                                                                                                                                                                                                                              aCGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGaa                                                                                                                                                                                                                                        >  5:138141/1‑77 (MQ=255)
                                                                                                                                                                                                                              aCGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGaa                                                                                                                                                                                                                                        <  6:138141/77‑1 (MQ=255)
                                                                                                                                                                                                                              aCGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAAt                                                                                                                                                                                                                >  6:569476/1‑101 (MQ=255)
                                                                                                                                                                                                                              aCGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAAt                                                                                                                                                                                                                <  5:569476/101‑1 (MQ=255)
                                                                                                                                                                                                                                     aaCGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCAGCGCAGTTAGCACGCCCGGACGGTTTTCGTgg                           <  2:498310/275‑1 (MQ=255)
                                                                                                                                                                                                                                      aCGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCtt                                                                                                                                                                                               <  7:1003963/110‑1 (MQ=255)
                                                                                                                                                                                                                                                 aaGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCAGCGCAGTTAGCACGCCCGGACGGTTTTCGTGGATGTGCATCAg                >  2:326403/1‑274 (MQ=255)
                                                                                                                                                                                                                                                         aaaaaaCGGGCAAGTCAGTGACGTGCCCGTTAATTTTCAGAGAAGGGGaa                                                                                                                                                                                                                                        <  7:274012/50‑1 (MQ=255)
                                                                                                                                                                                                                                                         aaaaaaCGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAg                                                                                                                                                                                                                                             >  8:375365/1‑45 (MQ=255)
                                                                                                                                                                                                                                                         aaaaaaCGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAg                                                                                                                                                                                                                                             <  7:375365/45‑1 (MQ=255)
                                                                                                                                                                                                                                                         aaaaaaCGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGaa                                                                                                                                                                                                                                        >  8:274012/1‑50 (MQ=255)
                                                                                                                                                                                                                                                          atcaaCGGGCAGGTCACTGACTTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTtcgt                                                                                                                                                           <  8:989248/123‑1 (MQ=255)
                                                                                                                                                                                                                                                          atcaaCGGGCAGGTCACTGACTTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTtcgt                                                                                                                                                           >  7:989248/4‑126 (MQ=255)
                                                                                                                                                                                                                                                           aaaaCGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGGACAGCAGACGGGCGCGAATGGTACCCGTAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCAGCGCAGTTAGCACGCCCGGACGGTTTTCGTGGATGTGCATCAGACGACGCCCa      <  2:127793/274‑1 (MQ=255)
                                                                                                                                                                                                                                                              aCGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCAGCGCAGTTAGCACGCCCGGACGGTTTTCGTGGATGTGCATCAGACGACGCCCACCgg  >  2:67847/1‑274 (MQ=255)
                                                                                                                                                                                                                                                               cGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGc                                                                                                                                                                                                 >  8:840701/1‑83 (MQ=255)
                                                                                                                                                                                                                                                               cGGGCAAGTCAGTGACCTGCCCGGTGATTTTCAGAGAAGGGGAATGAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGc                                                                                                                                                                                                 <  7:840701/83‑1 (MQ=255)
                                                                                                                                                                                                                                                                   cAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATaa                                                                                                                                 <  8:362577/143‑1 (MQ=255)
                                                                                                                                                                                                                                                                      gTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCttttt                                                                                                                                                                          >  5:309090/1‑99 (MQ=255)
                                                                                                                                                                                                                                                                      gTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCttttt                                                                                                                                                                          <  6:309090/99‑1 (MQ=255)
                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                 
GAGACTGACTGTTAATAAGCGCTGAAACTTATGAGTAACAGTACAATCAGTATGATGACAAGTCGCATCATAACCCTTCTCCTTCAAGCCCTCGCTTCGGTGAGGGCTTTACCGTTACAGCCCCATGCTGCCCTGCCATCGTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCAGCGCAGTTAGCACGCCCGGACGGTTTTCGTGGATGTGCATCAGACGACGCCCACCGT  >  NC_000913/3056879‑3057407

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: