Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,282,832 C→G 7.1% intergenic (+31/+108) 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,8320CG 7.1% 540.5 / 26.4 155intergenic (+31/+108)rplY/yejK50S ribosomal subunit protein L25/spermidine nucleoid‑associated protein
Reads supporting (aligned to +/- strand):  ref base C (75/69);  new base G (4/7);  total (79/76)
Fisher's exact test for biased strand distribution p-value = 3.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

Only 100 of 160 total aligned reads displayed.

CTACAGCGAAGTTCTGACCATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTC  >  NC_000913/2282684‑2282975
                                                                                                                                                    |                                                                                                                                                
cTACAGCGAAGTTCTGACCATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTcc                                                                                                                                                  <  3:313414/149‑1 (MQ=255)
 tACAGCGAAGTTCTGACCATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCcg                                                                                                                                                 >  1:145618/1‑149 (MQ=255)
   cAGCGAAGTTCTGACCATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCcg‑gc                                                                                                                                              <  3:286458/149‑1 (MQ=255)
     gCGAAGTTCTGACCATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCcg‑gcgg                                                                                                                                            <  3:1630826/149‑1 (MQ=255)
      cGAAGTTCTGACCATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCggg                                                                                                                                           <  2:1331691/149‑1 (MQ=255)
       gAAGTTCTGACCATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCgggg                                                                                                                                          >  3:490156/1‑149 (MQ=255)
         aGTTCTGACCATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGtt                                                                                                                                        <  4:1069810/149‑1 (MQ=255)
           ttCTGACCATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGtttt                                                                                                                                      >  1:1321922/1‑149 (MQ=255)
           ttCTGACCATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGtttt                                                                                                                                      >  3:1329174/1‑149 (MQ=255)
             cTGACCATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGtttttt                                                                                                                                    <  4:928254/149‑1 (MQ=255)
                aCCATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTa                                                                                                                                 <  2:769024/149‑1 (MQ=255)
                 ccATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTAt                                                                                                                                >  3:534759/1‑149 (MQ=255)
                   aTCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGtttttt                                                                                                                                    <  1:880074/143‑1 (MQ=255)
                   aTCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCt                                                                                                                              >  1:1404932/1‑149 (MQ=255)
                    tCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTg                                                                                                                             >  3:1563804/1‑149 (MQ=255)
                       ttgttgACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAg                                                                                                                          >  1:18965/1‑149 (MQ=255)
                        tgttgACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGa                                                                                                                         <  3:140771/149‑1 (MQ=255)
                           tgACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTa                                                                                                                      <  1:222837/149‑1 (MQ=255)
                              cGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGc                                                                                                                   <  2:537830/149‑1 (MQ=255)
                               ggTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGcc                                                                                                                  <  3:403555/149‑1 (MQ=255)
                               ggTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGcc                                                                                                                  <  3:937775/149‑1 (MQ=255)
                                gTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCt                                                                                                                 <  1:1488132/149‑1 (MQ=255)
                                 tAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTg                                                                                                                >  3:378238/1‑149 (MQ=255)
                                   aaGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGAt                                                                                                              >  2:1025775/1‑149 (MQ=255)
                                    aGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCCGGA‑GCGGGGtttt                                                                                                                                      <  1:1091386/124‑1 (MQ=255)
                                      aaaTCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTAc                                                                                                           <  4:711456/149‑1 (MQ=255)
                                       aaTCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACg                                                                                                          >  2:1318955/1‑149 (MQ=255)
                                         tCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGGACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGt                                                                                                        >  3:597751/1‑149 (MQ=255)
                                         tCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGt                                                                                                        >  4:1637784/1‑149 (MQ=255)
                                         tCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGt                                                                                                        >  4:1295031/1‑149 (MQ=255)
                                         tCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGt                                                                                                        >  2:16447/1‑149 (MQ=255)
                                          cAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTa                                                                                                       >  2:451368/1‑149 (MQ=255)
                                           aaaGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTGTTTGTATCTACAGATTATGCCTGATTACGGTAt                                                                                                      >  4:238634/1‑149 (MQ=255)
                                               ttAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCCGGA‑GCGGGGttttttacaa                                                                                                                                <  2:975683/119‑5 (MQ=38)
                                               ttAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCCGGA‑GCGGGGttttttacaa                                                                                                                                >  1:975683/1‑115 (MQ=38)
                                                  aaGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTAtt                                                                                               <  2:581446/149‑1 (MQ=255)
                                                   aGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTAttt                                                                                              <  4:1337030/149‑1 (MQ=255)
                                                   aGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTAttt                                                                                              >  3:996739/1‑149 (MQ=255)
                                                    gCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTAtttg                                                                                             <  2:93225/149‑2 (MQ=255)
                                                      tCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTg                                                                                                    >  2:1217476/1‑140 (MQ=255)
                                                      tCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTg                                                                                                    <  1:1217476/140‑1 (MQ=255)
                                                      tCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTAtttttt                                                                                           <  1:358303/149‑1 (MQ=255)
                                                      tCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTAtttttt                                                                                           <  3:187220/149‑1 (MQ=255)
                                                         ggACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCa                                                                                        >  1:546687/1‑149 (MQ=255)
                                                            cGTACTGCGTCACCCGTACAAACAGAGGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAACCCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGc                                                                                     <  2:877593/149‑1 (MQ=255)
                                                               aCAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAg                                                                                  >  4:593163/1‑149 (MQ=255)
                                                                 aGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGAt                                                                                >  2:1630736/1‑149 (MQ=255)
                                                                    gTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAg                                                                             <  3:905455/149‑1 (MQ=255)
                                                                    gTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAg                                                                             >  3:308969/1‑149 (MQ=255)
                                                                       aCCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCg                                                                          >  3:595562/1‑149 (MQ=255)
                                                                       aCCCGCACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCg                                                                          >  3:1624963/1‑149 (MQ=255)
                                                                         ccGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGa                                                                                 >  4:1187849/1‑140 (MQ=255)
                                                                         ccGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGa                                                                                 <  3:1187849/140‑1 (MQ=255)
                                                                             aCAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTTCAGATTATGCCTGATTACGGTAtt                                                                                                     >  1:5519/1‑116 (MQ=255)
                                                                               aaaCCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCAcgcc                                                                  >  3:1667089/1‑149 (MQ=255)
                                                                                      aGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATGCGGCGTGAACGCCTTATCTg                                                           <  1:16447/149‑1 (MQ=255)
                                                                                               acaTCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACTCCGCATCTGGCATTCAGa                                                  >  3:1508813/1‑149 (MQ=255)
                                                                                                caTCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGaa                                                 <  4:476363/149‑1 (MQ=255)
                                                                                                caTCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGaa                                                 <  4:1624963/149‑1 (MQ=255)
                                                                                                caTCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGaa                                                 <  3:378749/149‑1 (MQ=255)
                                                                                                  tCGACTTCGTTCGCGCTTACTTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAaca                                               <  4:1150008/149‑1 (MQ=255)
                                                                                                     aCTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCt                                                                                                                 <  3:88927/80‑1 (MQ=255)
                                                                                                     aCTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACaa                                            >  3:68932/1‑149 (MQ=255)
                                                                                                      cTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAg                                                                             >  4:862673/1‑115 (MQ=255)
                                                                                                       ttcgttcgCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACaaaa                                          >  4:717599/1‑149 (MQ=255)
                                                                                                       ttcgttcgCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACaaaa                                          >  2:34087/1‑149 (MQ=255)
                                                                                                       ttcgttcgCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACaaaa                                          >  1:1300249/1‑149 (MQ=255)
                                                                                                           ttcgCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGa                                      <  4:1329174/149‑1 (MQ=255)
                                                                                                           ttcgCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGa                                      <  4:378238/149‑1 (MQ=255)
                                                                                                           ttcgCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGa                                      <  3:1284828/149‑1 (MQ=255)
                                                                                                                gcTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAAt                                 <  2:1414720/149‑1 (MQ=255)
                                                                                                                gcTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAAt                                 >  3:1223696/1‑149 (MQ=255)
                                                                                                                 cTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAAtt                                <  4:308969/149‑1 (MQ=255)
                                                                                                                 cTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAAtt                                <  4:917046/149‑1 (MQ=255)
                                                                                                                  ttAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGaa                                                 >  1:1075125/1‑131 (MQ=255)
                                                                                                                  ttAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGaa                                                 <  2:1075125/131‑1 (MQ=255)
                                                                                                                           gAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGAc                      <  1:34087/149‑1 (MQ=255)
                                                                                                                             aTAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGt                    >  4:1419518/1‑149 (MQ=255)
                                                                                                                              tAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTc                   <  3:1516533/149‑1 (MQ=255)
                                                                                                                                 gatacaaaaaaCCCCGCCGGA‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGaa                                                 <  3:11528/111‑1 (MQ=38)
                                                                                                                                 gatacaaaaaaCCCCGC‑CGGAGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGaa                                                 >  4:11528/6‑116 (MQ=38)
                                                                                                                                 gTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCgg                >  1:1022240/1‑149 (MQ=255)
                                                                                                                                     tAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTc                   >  1:294393/1‑142 (MQ=255)
                                                                                                                                     tAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTc                   <  2:294393/142‑1 (MQ=255)
                                                                                                                                     tAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCt            >  4:610732/1‑149 (MQ=255)
                                                                                                                                       aaaaaCCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGc          >  4:786816/1‑149 (MQ=255)
                                                                                                                                       aaaaaCCCCGCTCCG‑GCGGGGTTTTTTGTATATGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGAGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGc          <  4:972261/149‑1 (MQ=255)
                                                                                                                                       aaaaaCCCCGCCGGA‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGa                                      >  2:1583950/1‑121 (MQ=38)
                                                                                                                                       aaaaaCCCCGCCGGA‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGa                                      <  1:1583950/121‑1 (MQ=255)
                                                                                                                                         aaaCCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCaa        <  1:771539/149‑1 (MQ=255)
                                                                                                                                            ccccGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTg     <  2:832021/149‑1 (MQ=255)
                                                                                                                                             cccGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTgg    >  4:1485597/1‑149 (MQ=255)
                                                                                                                                             cccGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTgg    >  4:79543/1‑149 (MQ=255)
                                                                                                                                              ccGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCa                                                    >  4:151961/1‑100 (MQ=255)
                                                                                                                                              ccGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCa                                                    <  3:151961/100‑1 (MQ=255)
                                                                                                                                              ccGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGt   >  2:785837/1‑149 (MQ=255)
                                                                                                                                               cGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGTGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTc  <  3:1295031/149‑1 (MQ=255)
                                                                                                                                               cGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTc  <  2:1300249/149‑1 (MQ=255)
                                                                                                                                               cGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTc  <  2:546687/149‑1 (MQ=255)
                                                                                                                                               cGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTc  <  1:1630736/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
CTACAGCGAAGTTCTGACCATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTC  >  NC_000913/2282684‑2282975

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