Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,282,830 T→C 5.5% intergenic (+29/+110) rplY → / ← yejK 50S ribosomal subunit protein L25/spermidine nucleoid‑associated protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,282,8300TC 5.5% 487.1 / 10.9 146intergenic (+29/+110)rplY/yejK50S ribosomal subunit protein L25/spermidine nucleoid‑associated protein
Reads supporting (aligned to +/- strand):  ref base T (74/63);  new base C (5/3);  total (80/66)
Fisher's exact test for biased strand distribution p-value = 7.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 147 total aligned reads displayed.

CTACAGCGAAGTTCTGACCATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCG  >  NC_000913/2282684‑2282978
                                                                                                                                                  |                                                                                                                                                     
cTACAGCGAAGTTCTGACCATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTcc                                                                                                                                                     >  2:317250/1‑149 (MQ=255)
     gCGAAGTTCTGACCATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCcg‑gcgg                                                                                                                                               <  1:317250/149‑1 (MQ=255)
     gCGAAGTTCTGACCATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCcg‑gcgg                                                                                                                                               >  3:543457/1‑149 (MQ=255)
     gCGAAGTTCTGACCATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCcg‑gcgg                                                                                                                                               <  1:1315957/149‑1 (MQ=255)
      cGAAGTTCTGACCATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCggg                                                                                                                                              <  1:531901/149‑1 (MQ=255)
          gTTCTGACCATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGttt                                                                                                                                          <  1:1019564/149‑1 (MQ=255)
          gTTCTGACCATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGttt                                                                                                                                          <  2:336062/149‑1 (MQ=255)
           ttCTGACCATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGtttt                                                                                                                                         >  2:1088208/1‑149 (MQ=255)
             cTGACCATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGtttttt                                                                                                                                       <  3:434836/149‑1 (MQ=255)
             cTGACCATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGtttttt                                                                                                                                       <  3:1421068/149‑1 (MQ=255)
             cTGACCATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGtttttt                                                                                                                                       <  4:178531/149‑1 (MQ=255)
              tGACCATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCCGGA‑GCGGGGtttttta                                                                                                                                      <  3:915572/149‑2 (MQ=255)
                 ccATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTAt                                                                                                                                   >  4:20590/1‑149 (MQ=255)
                   aTCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCt                                                                                                                                 >  3:595152/1‑149 (MQ=255)
                      gttgttGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCa                                                                                                                              <  1:759384/149‑1 (MQ=255)
                       ttgttgACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAg                                                                                                                             >  2:1472492/1‑149 (MQ=255)
                              cGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGc                                                                                                                      <  3:1279229/149‑1 (MQ=255)
                                gTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGtttttt                                                                                                                                       <  1:17628/130‑1 (MQ=255)
                                gTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGtttttt                                                                                                                                       >  2:17628/1‑130 (MQ=255)
                                  aaaGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTCCAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGa                                                                                                                  <  3:131689/149‑1 (MQ=255)
                                  aaaGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACTGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGa                                                                                                                  >  2:1154071/1‑149 (MQ=255)
                                       aaTCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACg                                                                                                             <  1:1378805/149‑1 (MQ=255)
                                        aTCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACgg                                                                                                            >  3:269658/1‑149 (MQ=255)
                                        aTCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACgg                                                                                                            >  4:925206/1‑149 (MQ=255)
                                         tCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGt                                                                                                           >  1:1047577/1‑149 (MQ=255)
                                          cAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTa                                                                                                          >  3:1437664/1‑149 (MQ=255)
                                          cAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTa                                                                                                          <  3:820138/149‑1 (MQ=255)
                                          cAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTa                                                                                                          <  1:918218/149‑1 (MQ=255)
                                           aatgTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCTGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTAt                                                                                                         >  4:1296612/4‑149 (MQ=255)
                                           aaaGTTAAAGGTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTAt                                                                                                         >  4:938007/1‑149 (MQ=255)
                                            aaGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTAtt                                                                                                        <  3:807037/149‑1 (MQ=255)
                                               ttAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTAtt                                                                                                        <  4:1410678/146‑1 (MQ=255)
                                               ttAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCt                                                                                                     >  1:994453/1‑149 (MQ=255)
                                               ttAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCt                                                                                                     >  4:1303581/1‑149 (MQ=255)
                                                   aGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTAttt                                                                                                 >  2:1120284/1‑149 (MQ=255)
                                                    gCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTAtttt                                                                                                <  2:203264/149‑1 (MQ=255)
                                                    gCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTAtttt                                                                                                >  2:681924/1‑149 (MQ=255)
                                                     cTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTAttttt                                                                                               >  4:1136619/1‑149 (MQ=255)
                                                      tCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTa                                                                                                                                    <  4:724186/111‑1 (MQ=255)
                                                      tCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTAtttttt                                                                                              >  1:1018603/1‑149 (MQ=255)
                                                      tCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTAtttttt                                                                                              <  1:1074231/149‑1 (MQ=255)
                                                      tCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTAtttttt                                                                                              >  2:1317448/1‑149 (MQ=255)
                                                          gACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAg                                                                                          <  3:143839/149‑1 (MQ=255)
                                                             gTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGcc                                                                                       >  1:10794/1‑149 (MQ=255)
                                                              tACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCa                                                                                      <  3:925206/149‑1 (MQ=255)
                                                               aCAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAg                                                                                     >  1:143642/1‑149 (MQ=255)
                                                                  gCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATa                                                                                  >  3:283672/1‑149 (MQ=255)
                                                                   cGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATaa                                                                                 <  4:543457/149‑1 (MQ=255)
                                                                   cGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATaa                                                                                 <  4:661421/149‑1 (MQ=255)
                                                                      cACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGc                                                                              <  2:10794/149‑1 (MQ=255)
                                                                       aCCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCg                                                                             <  1:352688/149‑1 (MQ=255)
                                                                         ccGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGtt                                                                           >  2:1111875/1‑149 (MQ=255)
                                                                         ccGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGtt                                                                           >  2:87967/1‑149 (MQ=255)
                                                                              cAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCAcgc                                                                      <  2:866047/149‑1 (MQ=255)
                                                                                 aCCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCAcgccgc                                                                   >  1:775623/1‑149 (MQ=255)
                                                                                 aCCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCAcgccgc                                                                   <  1:490186/149‑1 (MQ=255)
                                                                                 aCCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTGTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCAcgccgc                                                                   >  4:316549/1‑149 (MQ=255)
                                                                                   cGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATGCGGCGTGAAcgccttat                                                                 <  3:721788/149‑5 (MQ=255)
                                                                                   cGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCAt                                                                 >  3:1422485/1‑149 (MQ=255)
                                                                                     aaGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCt                                                               <  1:1108152/149‑1 (MQ=255)
                                                                                      aGCTGCAGCACATCGACTTCGTTCGCGCTTAACTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTg                                                              <  3:434039/149‑1 (MQ=255)
                                                                                            atcacaTCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTc                                                        <  4:1207122/147‑1 (MQ=255)
                                                                                            agcaCATCGACTTCGTTCGCGCTTAATCTCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTc                                                        >  4:894677/1‑149 (MQ=255)
                                                                                               acaTCGACTTCTGTCGCGCTTTATTGCGCAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGAGGGGTTTGTAGGTGCAAATTATGCCTGATTACGGTATTTCTACTTTGTGGATACCAGATACAGCCTTCCAACCCCATCTGGCATTAAGa                                                     >  4:858612/1‑149 (MQ=255)
                                                                                               acaTCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATGCGGCGTGAACGCCTTATCTGGc                                                            >  3:316784/1‑142 (MQ=255)
                                                                                                caTCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCGTCTGGCATTCAGaa                                                    <  1:1111875/149‑1 (MQ=255)
                                                                                                caTCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGaa                                                    <  3:574001/149‑1 (MQ=255)
                                                                                                caTCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGaa                                                    <  3:585285/149‑1 (MQ=255)
                                                                                                 aTCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAac                                                   >  4:345981/1‑149 (MQ=255)
                                                                                                     aCTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGc                                                                                        <  3:1222172/108‑1 (MQ=255)
                                                                                                     aCTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACaa                                               >  2:129071/1‑149 (MQ=255)
                                                                                                      cTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACaaa                                              >  2:90142/1‑149 (MQ=255)
                                                                                                       ttcgttcgCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACaaaa                                             >  1:874116/1‑149 (MQ=255)
                                                                                                           ttcgCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACaaaa                                             >  4:663311/1‑145 (MQ=255)
                                                                                                           ttcgCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGa                                         <  3:405834/149‑1 (MQ=255)
                                                                                                           ttcgCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGa                                         <  1:681924/149‑1 (MQ=255)
                                                                                                                gcTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATGTGCAGATTATGCCTGATTACGGTATTGCTATTATTTGCAGGCCAGCTAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAAt                                    <  4:1388668/149‑1 (MQ=255)
                                                                                                                 cTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAAtt                                   <  1:122984/149‑1 (MQ=255)
                                                                                                                  ttAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGc                                                                              <  2:1186384/105‑1 (MQ=255)
                                                                                                                  ttAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGc                                                                              >  1:1186384/1‑105 (MQ=255)
                                                                                                                    aaTTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTgcc                                >  4:506825/1‑149 (MQ=255)
                                                                                                                     aTTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGaa                                                    >  4:1427421/1‑128 (MQ=255)
                                                                                                                          tGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCATTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGa                          >  2:418333/1‑149 (MQ=255)
                                                                                                                            aaTAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACg                        >  3:1037328/1‑149 (MQ=255)
                                                                                                                                 gTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGAGAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCgg                   >  2:585913/1‑149 (MQ=255)
                                                                                                                                  ttGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGgc                  >  4:73554/1‑149 (MQ=255)
                                                                                                                                      aaaaaaCCCCGCCGGA‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTg              >  3:1302110/1‑149 (MQ=255)
                                                                                                                                      aaaaaaCCCCGC‑CGGAGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTg              >  1:244874/1‑149 (MQ=39)
                                                                                                                                       aaaaaCCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAg                                                                                <  3:1102411/82‑1 (MQ=255)
                                                                                                                                       aaaaaCCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGc             >  1:581810/1‑149 (MQ=255)
                                                                                                                                       aaaaaCCCCGCCGGA‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGaa                                                    <  2:577355/110‑1 (MQ=38)
                                                                                                                                       aaaaaCCCCGCCGGA‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGc             >  3:1138113/1‑149 (MQ=255)
                                                                                                                                        aaaaCCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGGGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCa            >  2:407505/1‑149 (MQ=255)
                                                                                                                                         aaaCCCCGCCGGA‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCaa           >  3:553110/1‑149 (MQ=255)
                                                                                                                                          aaCCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAAt          <  2:994453/149‑1 (MQ=255)
                                                                                                                                           aCCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAAtt         <  3:1282808/149‑1 (MQ=255)
                                                                                                                                           aCCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAAtt         >  4:1079813/1‑149 (MQ=255)
                                                                                                                                             cccGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCAcg                                                                       >  3:958724/1‑85 (MQ=255)
                                                                                                                                             cccGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTgg       >  1:96236/1‑149 (MQ=255)
                                                                                                                                                  tCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTcgcg  >  3:723668/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
CTACAGCGAAGTTCTGACCATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCG  >  NC_000913/2282684‑2282978

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