Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,648,359 G→T 6.0% noncoding (66/84 nt) agrB → antisense sRNA antitoxin for the DinQ toxin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,648,3590GT 6.0% 416.0 / 11.8 153noncoding (66/84 nt)agrBantisense sRNA antitoxin for the DinQ toxin
Reads supporting (aligned to +/- strand):  ref base G (82/61);  new base T (5/5);  total (87/66)
Fisher's exact test for biased strand distribution p-value = 7.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.98e-01

Only 100 of 153 total aligned reads displayed.

TATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCGAAGAGAGACACTACCT  >  NC_000913/3648212‑3648507
                                                                                                                                                    |                                                                                                                                                     
tatataCCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCa                                                                                                                                                     >  1:31023/1‑150 (MQ=255)
  tataCCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCaaa                                                                                                                                                   >  1:719684/1‑150 (MQ=255)
     aCCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAACCCGCCGCaagggg                                                                                                                                                >  1:553199/1‑146 (MQ=255)
     aCCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTc                                                                                                                                                >  1:1705265/1‑150 (MQ=255)
     aCCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTc                                                                                                                                                >  2:1354556/1‑150 (MQ=255)
        gAAGTGCTCCCTCCGCCAGCTCAAGAAATCGCTAATTCATGCAATGTTAGCCACTGTCTAATAGTATTGAGCTGTAAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGTCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑t                                                                                                                                            <  1:348758/150‑1 (MQ=37)
             gCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAATAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑tttttt                                                                                                                                       <  1:477191/150‑1 (MQ=255)
             gCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACTTTGCGTCGG‑tttttt                                                                                                                                       <  1:602636/150‑1 (MQ=38)
             gCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑tttttt                                                                                                                                       >  2:1596268/1‑150 (MQ=255)
             gCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑tttttt                                                                                                                                       >  2:410738/1‑150 (MQ=255)
               tccctccGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAAACCGACTTTGCGTCGG‑tttttta                                                                                                                                      <  1:1509267/150‑2 (MQ=38)
                ccctccGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTTAc                                                                                                                                    >  2:323148/1‑150 (MQ=255)
                 cctccGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGCCTGGAAAGT‑ATAAAACCGACGCAAAGTCGG‑TTTTTTTACg                                                                                                                                   <  1:1352635/150‑1 (MQ=255)
                    ccGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTTa                                                                                                                                     >  2:1650537/1‑145 (MQ=255)
                     cGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGCCGCAAAGTCGG‑TTTTTTTACGTCCt                                                                                                                               <  1:989401/150‑1 (MQ=255)
                     cGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGAAAATAGTATTGATCTGTTAGATAAGAACACTCTCACTCCAGCCAGAGCCACCAACCCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTTACGTCCt                                                                                                                               <  1:765699/150‑1 (MQ=37)
                         aGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTTACGTCCTGAtt                                                                                                                           >  2:620562/1‑150 (MQ=255)
                           cTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTTACGTCCTGATTCa                                                                                                                         >  2:1428515/1‑150 (MQ=255)
                             gaagaaATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTTACGTCCTGATTCAGa                                                                                                                       >  1:1256655/1‑150 (MQ=255)
                               agaaATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGCCGCAAAGTCGG‑TTTTTTTACGTCAAGATTCAGAcc                                                                                                                     >  2:1609355/1‑150 (MQ=255)
                                      gCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGATCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTTACGTCCTGATTCAGACCTCCTTTc                                                                                                              <  1:955547/150‑1 (MQ=255)
                                      gCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTTACGTCCTGATTCAGACCTCCTTTc                                                                                                              <  1:1686944/150‑1 (MQ=255)
                                        aaaTTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAAATCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTTACGTCCTGATTCAGACCTCCTTTCaa                                                                                                            <  1:895982/149‑1 (MQ=255)
                                         aaTTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTTACGTCCTGATTCAGACCTCCTTTCaaa                                                                                                           >  1:1037729/1‑150 (MQ=255)
                                          atccTTGCAATCTTAGCCCCAGGCTAATAGTATTGAGCCGTTAGATAAGACCTCTCTCACTCCAGCCAGATCCACCAACTCAGGGCCGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTTACGTCCTGATTCAGACCTCCTTTCAAAt                                                                                                          <  1:1732497/147‑1 (MQ=25)
                                           ttCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATg                                                                                                         <  1:423238/150‑1 (MQ=255)
                                             cTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGaa                                                                                                       >  2:749078/1‑150 (MQ=255)
                                              ttGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAAt                                                                                                      >  1:142172/1‑150 (MQ=255)
                                               tGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACTCAAAGTCGG‑TTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATa                                                                                                     >  1:24147/1‑150 (MQ=255)
                                                 tactattAGCCAGTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACTGACGCAAAGTCGG‑TTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGc                                                                                                   >  1:1201943/6‑150 (MQ=255)
                                                 cAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGc                                                                                                   >  2:738388/1‑150 (MQ=255)
                                                  actattAGCCAGTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGcc                                                                                                  >  2:157778/5‑150 (MQ=255)
                                                  actattAGCCAGTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTTACGTCCTGATACAGACCTCCTTTCAAATGAATAGcc                                                                                                  >  1:32725/5‑150 (MQ=255)
                                                  aaTGTTACCCACTGCCTAATAGTATTGCGCTGTACTATAATAACTATATCACTCCAGCCACAGCCACCACCTCAGGGATGGAATGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTTAGTTCCTGATTCAGACCTCATTTCAAATGAATAccc                                                                                                  <  2:1630686/150‑1 (MQ=18)
                                                     gTTAGCCACTGGCTATTAGTTTTGAGCTGTTAGATAAGAAGTCTCTCACTCCAGACAGAGCCACCAACTCATGTCTTTTAAGT‑AAAAAACCGACGCATAGTCGG‑TTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAAc                                                                                               <  1:510449/150‑1 (MQ=255)
                                                     gTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAAc                                                                                               <  1:539112/150‑1 (MQ=255)
                                                       tAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTc                                                                                             >  1:1159206/1‑150 (MQ=255)
                                                             cTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAAtt                                                                                       >  1:1024044/1‑150 (MQ=255)
                                                               ggCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTca                                                                                     >  2:346744/1‑150 (MQ=255)
                                                               ggCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTca                                                                                     >  1:1530801/1‑150 (MQ=255)
                                                                gCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTcac                                                                                    <  1:542371/150‑1 (MQ=255)
                                                                   aaTAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCAACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACAcc                                                                                 <  1:1446903/150‑1 (MQ=255)
                                                                      aGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTAt                                                                              >  1:1632537/1‑150 (MQ=255)
                                                                       gTATTGATGTGTTAGATGATTGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTAtt                                                                             <  2:168310/150‑1 (MQ=17)
                                                                         aTTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTAc                                                                           >  2:1202579/1‑150 (MQ=255)
                                                                         aTTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTAc                                                                           >  1:809317/1‑150 (MQ=255)
                                                                         aTTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTAc                                                                           <  1:496230/150‑1 (MQ=255)
                                                                            gAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCtt                                                                        >  2:1324514/1‑150 (MQ=255)
                                                                            gAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCtt                                                                        >  1:1402073/1‑150 (MQ=255)
                                                                             aGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTc                                                                       >  1:1590658/1‑150 (MQ=255)
                                                                                tGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACGAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCctc                                                                    >  2:205549/1‑150 (MQ=255)
                                                                                tGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCctc                                                                    >  1:499072/1‑150 (MQ=255)
                                                                                  ttAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTg                                                                  <  2:1321235/150‑1 (MQ=255)
                                                                                    aGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCa                                                                >  2:44812/1‑150 (MQ=255)
                                                                                     gATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCAc                                                               <  1:323148/150‑1 (MQ=255)
                                                                                       tAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACtt                                                             <  2:767647/150‑1 (MQ=255)
                                                                                       tAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACtt                                                             >  2:637726/1‑150 (MQ=255)
                                                                                        aaGAACTCTCTCACTCCAGCCAGAGCCACCAACTTAGGGCTGGATAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCTACTTAAAATTCACACCTATTACCTTCCTCTGCACTTa                                                            <  1:1158444/150‑1 (MQ=255)
                                                                                        aaGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTTACGTCCTGATTCAGCCCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTa                                                            >  2:871282/1‑150 (MQ=255)
                                                                                             ctctctCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCACCTCAAACTTCCCACCTAGTACCTTCCTGTGCGCGTAGACAt                                                       >  1:301987/1‑150 (MQ=21)
                                                                                             ctctctCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACAt                                                       >  1:373684/1‑150 (MQ=255)
                                                                                               ctctcACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTTACGTCCTGATTCAGCCCTCCGTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACCCATTc                                                     >  1:1121064/1‑150 (MQ=255)
                                                                                               ctctcACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCTCTTACACAttt                                                     <  1:620562/150‑2 (MQ=255)
                                                                                                tctcACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCg                                                    <  1:1629222/150‑1 (MQ=255)
                                                                                                tctcACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCg                                                    >  2:689105/1‑150 (MQ=255)
                                                                                                 ctcACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGt                                                   >  2:1363135/1‑150 (MQ=255)
                                                                                                   cACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTa                                                 <  2:719684/150‑1 (MQ=255)
                                                                                                   cACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAAATCACACCTATTACCTTCCTCTGCACTGACACATTCGTTa                                                 >  2:86821/1‑150 (MQ=255)
                                                                                                    aCTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTTACGTCCTGATTCAGACCGCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTaa                                                >  1:1142867/1‑150 (MQ=255)
                                                                                                       ccagccagAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGACTAGCCAACTCAAAATTCACACCTATTACCTTCCTATGCAGTTACACAGTCGTTAAGTc                                             >  1:342665/1‑150 (MQ=255)
                                                                                                          gccagAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCata                                          >  1:719351/1‑150 (MQ=255)
                                                                                                             agagCCAACAACTCATGGCTGGAAAGT‑AAAAAACCGAATCAAATTCGGTTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCatata                                        <  1:193051/150‑1 (MQ=255)
                                                                                                                    ccaACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTg                                <  2:1715906/150‑1 (MQ=255)
                                                                                                                      aaCTCAGGGCTGGAAAGT‑AAAAAACCGAGGCAAAGTCGG‑TTTTTTTACTTCCTGATTCAGACCTCCTTTCAAATTAATATCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGAc                              <  2:548150/150‑1 (MQ=255)
                                                                                                                       aCTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACt                             <  2:300840/150‑1 (MQ=255)
                                                                                                                       aCTCAGGGATGGAAAGT‑AAAAAACCGACTTTGCGTCGG‑TTTTTTTGCGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACt                             <  2:1081897/150‑1 (MQ=255)
                                                                                                                          cAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGAATTAt                          >  2:515887/1‑150 (MQ=255)
                                                                                                                          cAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCATCCTCTGCACTTACACAATCGTTAAGTCATATATGTTTTTGACTTaa                          >  2:670647/1‑149 (MQ=255)
                                                                                                                            gggCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATcc                        >  2:422137/1‑150 (MQ=255)
                                                                                                                             ggCTGGAAATT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCg                       <  2:435664/150‑1 (MQ=255)
                                                                                                                             ggCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCTTATATGTTTTTGACTTATCCg                       >  1:1609463/1‑150 (MQ=255)
                                                                                                                               cTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCTCTTACACATTCGTTAAGTCATATATGTTTTTGACTTATACGCt                     >  2:1091542/1‑150 (MQ=255)
                                                                                                                               cTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCt                     <  2:1158828/150‑1 (MQ=255)
                                                                                                                                  gAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCg                  <  1:1596268/150‑1 (MQ=255)
                                                                                                                                  gAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCg                  <  2:589529/150‑1 (MQ=255)
                                                                                                                                  gAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCg                  >  2:770587/1‑150 (MQ=255)
                                                                                                                                  gAAAGT‑AAAAAAACGACGCAAAGTCGG‑TTTTTTTACGTCCTGATTCAGACCTACTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCg                  <  2:813279/150‑1 (MQ=255)
                                                                                                                                     cgtaaaaaaaCCGACTTTGCGTCGG‑‑TTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCGAag               >  2:1423995/5‑150 (MQ=255)
                                                                                                                                     aGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCGAag               <  1:1557285/150‑1 (MQ=255)
                                                                                                                                     aGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCGAag               <  2:1463098/150‑1 (MQ=255)
                                                                                                                                      gtaaaaaaaCCGACTTTGCGTCGG‑TTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACGTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCGaat               >  2:1152644/4‑149 (MQ=255)
                                                                                                                                       taaaaaaaCCGACTTTGCGTCGG‑TTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATTTTTTTGACTTATCCGCTTCGAaga              >  1:1610767/3‑150 (MQ=255)
                                                                                                                                       t‑AAAAAACCGACGCAAAGTCGG‑TTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCAGCTTCGAagag             >  2:1668205/1‑150 (MQ=255)
                                                                                                                                          aaaaaCCGACGCAAAGTCGG‑TTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCGAagagag           <  1:1202579/150‑1 (MQ=255)
                                                                                                                                           aaaaCCGACGCAAAGTCGG‑TTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCGAAgagaga          >  2:1204390/1‑150 (MQ=255)
                                                                                                                                             aaCCGACGCAAAGTCGG‑TTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCGAAGAGAGaca        <  2:1479846/150‑1 (MQ=255)
                                                                                                                                               ccGACGCAAAGTCGG‑TTTTTTTACGTCCTGATTCAGACCTCATTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCGAAGAGAGACACt      <  2:553199/150‑1 (MQ=255)
                                                                                                                                                cGACGCAAAGTCGG‑TTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCGAAGAGAGACACTa     <  1:205549/150‑1 (MQ=255)
                                                                                                                                                 gACGCAAAGTCGG‑TTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCGAAGAGAGACACTAc    <  1:1431729/150‑1 (MQ=255)
                                                                                                                                                   cGCAAAGTCGG‑TTTTTTTACGTCCTGATTCAGAACTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTAATACCTTCCACTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCGAAGAGAGACACTACCt  <  2:308388/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCGAAGAGAGACACTACCT  >  NC_000913/3648212‑3648507

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑