Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,764,693 C→T 5.8% intergenic (‑307/+20) sufA ← / ← rydB Fe‑S cluster assembly protein/novel sRNA, function unknown

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,764,6930CT 5.8% 660.6 / 19.8 205intergenic (‑307/+20)sufA/rydBFe‑S cluster assembly protein/novel sRNA, function unknown
Reads supporting (aligned to +/- strand):  ref base C (94/99);  new base T (7/5);  total (101/104)
Fisher's exact test for biased strand distribution p-value = 5.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.64e-01

Only 100 of 205 total aligned reads displayed.

ATTAGGTATCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAAT  >  NC_000913/1764547‑1764839
                                                                                                                                                   |                                                                                                                                                  
aTTAGGTATCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGa                                                                                                                                                  >  4:1262210/1‑149 (MQ=255)
aTTAGGTATCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGa                                                                                                                                                  <  1:927196/149‑1 (MQ=255)
    ggTATCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGAcaca                                                                                                                                              >  3:761834/1‑149 (MQ=255)
    ggTATCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGAcaca                                                                                                                                              <  1:744695/149‑1 (MQ=255)
     gTATCTCATTGTTAGATAATGGAAATTATGTAAGGAGCATTAAGGCTGTGATAAAAAAGAAAAAGTTGTATTGAAAAAGACTATTTAAGAGATGGGTAAGAAGGTGCAGCGTTCAGTAATAAGAAAACCCTTAAGTCTGTGCGacacgg                                                                                                                                             <  4:164589/149‑3 (MQ=255)
         ctcATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCtt                                                                                                                                         <  4:912047/149‑1 (MQ=255)
         ctcATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCtt                                                                                                                                         <  1:44933/149‑1 (MQ=255)
         ctcATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCtt                                                                                                                                         >  3:1204140/1‑149 (MQ=255)
          tcATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAgg                                                                                                                                            <  3:997550/145‑1 (MQ=255)
          tcATGGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTa                                                                                                                                        >  1:1493058/1‑149 (MQ=255)
            aTTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAg                                                                                                                                      >  4:608359/1‑149 (MQ=255)
            aTTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAg                                                                                                                                      >  1:1526944/1‑149 (MQ=255)
            aTTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAg                                                                                                                                      <  3:157835/149‑1 (MQ=255)
            aTTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGCCTGTGTCGCACAGACTTAAg                                                                                                                                      <  2:538333/149‑1 (MQ=255)
            aTTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGATAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAg                                                                                                                                      <  1:1101479/149‑1 (MQ=255)
             ttGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAgg                                                                                                                                     <  2:484266/149‑1 (MQ=255)
              tGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAggg                                                                                                                                    <  1:481369/149‑1 (MQ=255)
                ttAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGt                                                                                                                                   >  4:735773/1‑148 (MQ=255)
                ttAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGttt                                                                                                                                 >  3:893397/1‑149 (MQ=255)
                      aaTGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTAcc                                                                                                                            >  3:1547342/1‑149 (MQ=255)
                         ggTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCa                                                                                                                         >  3:1309266/1‑149 (MQ=255)
                          gTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTAccc                                                                                                                           >  2:1525798/1‑146 (MQ=255)
                                atCTAAGGAGCGTTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGgcg                                                                                                                  >  3:165702/1‑149 (MQ=255)
                                 tCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGgcgc                                                                                                                 >  4:1199069/1‑149 (MQ=255)
                                   tAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAgg                                                                                                                                     >  2:976468/1‑127 (MQ=255)
                                   tAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGCCTGTGTCGCACAGACTTAAg                                                                                                                                      <  2:1609341/126‑1 (MQ=38)
                                   tAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGCCTGTGTCGCACAGACTTAAg                                                                                                                                      >  1:1609341/1‑126 (MQ=38)
                                            ttAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGCCTGTGTCGCACAGACTTAAGGGttt                                                                                                                                 >  4:1480026/1‑122 (MQ=38)
                                                      gatgaAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTg                                                                                            <  4:1001666/149‑1 (MQ=255)
                                                         gaaaaagaaaaaGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCt                                                                                         >  2:546176/1‑149 (MQ=255)
                                                           aaaagaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTa                                                                                       >  3:1472176/1‑149 (MQ=255)
                                                            aaagaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAg                                                                                      >  4:1438457/1‑149 (MQ=255)
                                                             aagaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGc                                                                                     <  2:1317098/149‑1 (MQ=255)
                                                               gaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGctc                                                                                   >  1:278578/1‑149 (MQ=255)
                                                                     gTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAg                                                                             >  1:767405/1‑149 (MQ=255)
                                                                      ttGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAgg                                                                            <  2:676605/149‑1 (MQ=255)
                                                                      ttGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAgg                                                                            <  4:1309266/149‑1 (MQ=255)
                                                                      ttGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAgg                                                                            >  4:184325/1‑149 (MQ=255)
                                                                       tGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAggg                                                                           <  4:1469141/149‑1 (MQ=255)
                                                                       tGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAggg                                                                           <  2:89067/149‑1 (MQ=255)
                                                                        gTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAgggg                                                                          <  4:292505/149‑1 (MQ=255)
                                                                          aTTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCg                                                                        >  3:947026/1‑149 (MQ=255)
                                                                           ttGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTAc                                                                                                                             >  3:934186/1‑95 (MQ=255)
                                                                           ttGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGa                                                                       >  1:1697147/1‑149 (MQ=255)
                                                                           ttGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGa                                                                       >  4:390671/1‑149 (MQ=255)
                                                                              aaaaTGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATaa                                                                    <  2:178259/149‑1 (MQ=255)
                                                                                 ctGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGaa                                                                 <  3:1495414/148‑1 (MQ=255)
                                                                                  tGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGaat                                                                >  3:1539366/1‑149 (MQ=255)
                                                                                  tGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCAACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGaat                                                                <  4:1181294/149‑1 (MQ=255)
                                                                                      tATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTAc                                                                                                                             >  1:21080/1‑84 (MQ=255)
                                                                                      tATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTAc                                                                                                                             <  2:21080/84‑1 (MQ=255)
                                                                                      tATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAAtc                                                            >  2:424673/1‑149 (MQ=255)
                                                                                         ttAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCa                                                         <  4:837487/149‑1 (MQ=255)
                                                                                           aaGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATa                                                       >  3:1594001/1‑149 (MQ=255)
                                                                                               gaTAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGc                                                   <  4:161581/149‑1 (MQ=255)
                                                                                                  aGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTaa                                                >  1:993091/1‑149 (MQ=255)
                                                                                                   ggTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAAc                                               <  4:1419382/149‑1 (MQ=255)
                                                                                                    gTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGaa                                                                              >  2:1191243/1‑117 (MQ=255)
                                                                                                    gTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGaa                                                                              <  1:1191243/117‑1 (MQ=255)
                                                                                                         aaaGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACccgccg                                         >  1:603172/1‑149 (MQ=255)
                                                                                                          aaGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGt                                        >  2:867723/1‑149 (MQ=255)
                                                                                                           aGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGtt                                       <  4:1594001/149‑1 (MQ=255)
                                                                                                               cAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTaa                                   <  3:1421986/149‑1 (MQ=255)
                                                                                                               cAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTaa                                   <  3:23975/149‑1 (MQ=255)
                                                                                                               cAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTaa                                   <  3:1262210/149‑1 (MQ=255)
                                                                                                                aGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAgggg                                                                          <  1:867603/109‑1 (MQ=255)
                                                                                                                 gCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAaca                                 >  4:1372459/1‑149 (MQ=255)
                                                                                                                    ttCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGTTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAaca                                 <  3:1064303/146‑1 (MQ=255)
                                                                                                                    ttCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAaca                                 >  4:1064303/1‑146 (MQ=255)
                                                                                                                    ttCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACaa                              >  4:1061261/1‑149 (MQ=255)
                                                                                                                     tCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACaaa                             <  3:588332/149‑1 (MQ=255)
                                                                                                                     tCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACaaa                             <  3:1061261/149‑1 (MQ=255)
                                                                                                                        gAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAGGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTg                          >  3:1209874/1‑149 (MQ=255)
                                                                                                                         aaaTAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGc                         <  4:947026/149‑1 (MQ=255)
                                                                                                                         aaaTAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGc                         >  3:33710/1‑149 (MQ=255)
                                                                                                                           aTAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGa                       <  4:890518/149‑1 (MQ=255)
                                                                                                                           aTAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGa                       <  2:1256960/149‑1 (MQ=255)
                                                                                                                            tAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGAt                      >  3:868831/1‑149 (MQ=255)
                                                                                                                            tAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCCACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCACAAAGCTAACCCGCCGTTTTAACACAAACTGCGAt                      <  4:1222075/149‑1 (MQ=255)
                                                                                                                              aGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTa                    <  3:184325/149‑1 (MQ=255)
                                                                                                                               gAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGa                       <  3:509516/145‑1 (MQ=255)
                                                                                                                               gAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAg                   >  1:711200/1‑149 (MQ=255)
                                                                                                                               gAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAg                   >  3:452799/1‑149 (MQ=255)
                                                                                                                                aaaaCCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGt                  <  4:1539366/149‑1 (MQ=255)
                                                                                                                                     ccTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGc                         >  4:53545/1‑137 (MQ=255)
                                                                                                                                     ccTTAAGCCTGTGTCGCACAGACTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGtattat             >  2:352703/1‑149 (MQ=255)
                                                                                                                                      cTTAAGCCTGTGTCGCACAGACTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGaa                                                                              <  2:1076171/83‑1 (MQ=25)
                                                                                                                                      cTTAAGCCTGTGTCGCACAGACTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGaa                                                                              >  1:1076171/1‑83 (MQ=25)
                                                                                                                                       ttAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTaa                                                <  4:509383/112‑1 (MQ=255)
                                                                                                                                       ttAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTAttt           >  3:15040/1‑149 (MQ=255)
                                                                                                                                       ttAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTAttt           >  4:41801/1‑149 (MQ=255)
                                                                                                                                       ttAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTAttt           >  2:267217/1‑149 (MQ=255)
                                                                                                                                       ttAAGGCTGTGCGACACAGGCTTAAGGTTTTCGACCCCAACCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAAACACATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTAttt           <  4:27550/149‑1 (MQ=255)
                                                                                                                                       ttAAGCCTGTGTCGCACAGACTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTAttt           >  4:93071/1‑149 (MQ=255)
                                                                                                                                        tAAGCCTGTGTCGCACAGACTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTAtttt          <  1:352703/149‑1 (MQ=255)
                                                                                                                                         aaGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTAttttt         >  2:608092/1‑149 (MQ=255)
                                                                                                                                         aaGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTAttttt         <  4:1575065/149‑1 (MQ=255)
                                                                                                                                            tCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGaa      <  4:15040/149‑1 (MQ=255)
                                                                                                                                            tCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGaa      >  4:821809/1‑149 (MQ=255)
                                                                                                                                               gtgCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACaa   <  3:706539/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
ATTAGGTATCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAAT  >  NC_000913/1764547‑1764839

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