Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,648,361 A→T 5.2% noncoding (68/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,3610AT 5.2% 447.4 / 6.8 149noncoding (68/84 nt)agrBantisense sRNA antitoxin for the DinQ toxin
Reads supporting (aligned to +/- strand):  ref base A (79/62);  new base T (4/4);  total (83/66)
Fisher's exact test for biased strand distribution p-value = 7.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

Only 100 of 154 total aligned reads displayed.

TATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCGAAGAGAGACACTACCTGCA  >  NC_000913/3648212‑3648510
                                                                                                                                                      |                                                                                                                                                     
  tataCCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCaaa                                                                                                                                                     >  1:719684/1‑150 (MQ=255)
     aCCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAACCCGCCGCaagggg                                                                                                                                                  >  1:553199/1‑146 (MQ=255)
     aCCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTc                                                                                                                                                  >  1:1705265/1‑150 (MQ=255)
      ccGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCg                                                                                                                                                 >  1:589529/1‑150 (MQ=255)
        gAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGt                                                                                                                                               >  1:1539097/1‑150 (MQ=255)
        gAAGTGCTCCCTCCGCCAGCTCAAGAAATCGCTAATTCATGCAATGTTAGCCACTGTCTAATAGTATTGAGCTGTAAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGTCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGt                                                                                                                                               <  1:348758/150‑1 (MQ=37)
             gCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGtttttt                                                                                                                                          >  2:410738/1‑150 (MQ=255)
              ctccctccGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGttttttt                                                                                                                                         <  2:575132/150‑1 (MQ=255)
               tccctccGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAAACCGACTTTGCGTCGGtttttta                                                                                                                                         <  1:1509267/150‑2 (MQ=38)
                ccctccGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTAc                                                                                                                                       >  2:323148/1‑150 (MQ=255)
                ccctccGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTAc                                                                                                                                       >  2:1629222/1‑150 (MQ=255)
                 cctccGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGCCTGGAAAGT‑ATAAAACCGACGCAAAGTCGGTTTTTTTACg                                                                                                                                      <  1:1352635/150‑1 (MQ=255)
                    ccGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTa                                                                                                                                        >  2:1650537/1‑145 (MQ=255)
                     cGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGCCGCAAAGTCGGTTTTTTTACGTCCt                                                                                                                                  <  1:989401/150‑1 (MQ=255)
                      gCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTg                                                                                                                                 >  1:1463098/1‑150 (MQ=255)
                         aGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGAtt                                                                                                                              >  2:620562/1‑150 (MQ=255)
                           cTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCa                                                                                                                            >  2:1428515/1‑150 (MQ=255)
                             gaagaaATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGa                                                                                                                          >  1:1256655/1‑150 (MQ=255)
                               agaaATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGCCGCAAAGTCGGTTTTTTTACGTCAAGATTCAGAcc                                                                                                                        >  2:1609355/1‑150 (MQ=255)
                                      gCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTc                                                                                                                 <  1:1686944/150‑1 (MQ=255)
                                       cTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCa                                                                                                                <  1:341463/150‑1 (MQ=255)
                                        aaaTTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAAATCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCaa                                                                                                               <  1:895982/149‑1 (MQ=255)
                                         aaTTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCaaa                                                                                                              >  1:1037729/1‑150 (MQ=255)
                                          atccTTGCAATCTTAGCCCCAGGCTAATAGTATTGAGCCGTTAGATAAGACCTCTCTCACTCCAGCCAGATCCACCAACTCAGGGCCGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAAt                                                                                                             <  1:1732497/147‑1 (MQ=25)
                                           ttCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATg                                                                                                            <  1:423238/150‑1 (MQ=255)
                                             cTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGaa                                                                                                          >  2:749078/1‑150 (MQ=255)
                                              ttGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAAt                                                                                                         >  1:142172/1‑150 (MQ=255)
                                               tGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAAtg                                                                                                        >  2:1654181/1‑149 (MQ=255)
                                               tGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATa                                                                                                        >  1:300840/1‑150 (MQ=255)
                                                atactattAGCCAGTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAg                                                                                                       >  1:989585/7‑150 (MQ=255)
                                                 tactattAGCCAGTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACTGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGc                                                                                                      >  1:1201943/6‑150 (MQ=255)
                                                 cAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGc                                                                                                      >  2:738388/1‑150 (MQ=255)
                                                  actattAGCCAGTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGcc                                                                                                     >  2:157778/5‑150 (MQ=255)
                                                     gTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCTAc                                                                                                  <  1:539115/150‑1 (MQ=255)
                                                     gTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAAc                                                                                                  <  1:539112/150‑1 (MQ=255)
                                                       tAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTc                                                                                                >  1:1159206/1‑150 (MQ=255)
                                                             cTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAAtt                                                                                          >  1:1024044/1‑150 (MQ=255)
                                                               ggCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTca                                                                                        >  2:346744/1‑150 (MQ=255)
                                                               ggCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTca                                                                                        >  1:935262/1‑150 (MQ=255)
                                                                      aGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTAt                                                                                 >  1:1632537/1‑150 (MQ=255)
                                                                       gTATTGATGTGTTAGATGATTGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTAtt                                                                                <  2:168310/150‑1 (MQ=17)
                                                                         aTTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTAc                                                                              >  1:809317/1‑150 (MQ=255)
                                                                          ttGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTAcc                                                                             <  2:525846/150‑1 (MQ=255)
                                                                            gAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCtt                                                                           >  1:1402073/1‑150 (MQ=255)
                                                                            gAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCtt                                                                           >  2:1324514/1‑150 (MQ=255)
                                                                             aGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTc                                                                          >  1:1590658/1‑150 (MQ=255)
                                                                                tGTTAGATAATTCCTATCTCACTCCAGGCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCctc                                                                       <  2:1326157/150‑1 (MQ=25)
                                                                                tGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACGAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCctc                                                                       >  2:205549/1‑150 (MQ=255)
                                                                                tGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCctc                                                                       >  1:499072/1‑150 (MQ=255)
                                                                                    aGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCa                                                                   >  2:44812/1‑150 (MQ=255)
                                                                                     gATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCAc                                                                  <  1:323148/150‑1 (MQ=255)
                                                                                       tAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACtt                                                                >  2:158610/1‑150 (MQ=255)
                                                                                       tAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACtt                                                                >  2:637726/1‑150 (MQ=255)
                                                                                       tAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACtt                                                                <  2:767647/150‑1 (MQ=255)
                                                                                        aaGAACTCTCTCACTCCAGCCAGAGCCACCAACTTAGGGCTGGATAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCTACTTAAAATTCACACCTATTACCTTCCTCTGCACTTa                                                               <  1:1158444/150‑1 (MQ=255)
                                                                                        aaGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGCCCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTa                                                               >  2:871282/1‑150 (MQ=255)
                                                                                            aCTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTAcaca                                                           >  1:201741/1‑150 (MQ=255)
                                                                                             ctctctCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCACCTCAAACTTCCCACCTAGTACCTTCCTGTGCGCGTAGACAt                                                          >  1:301987/1‑150 (MQ=21)
                                                                                             ctctctCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACAt                                                          >  1:373684/1‑150 (MQ=255)
                                                                                               ctctcACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGCCCTCCGTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACCCATTc                                                        >  1:1121064/1‑150 (MQ=255)
                                                                                               ctctcACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCTCTTACACAttt                                                        <  1:620562/150‑2 (MQ=255)
                                                                                               ctctcACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACGCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTc                                                        >  2:1701843/1‑150 (MQ=255)
                                                                                                tctcACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCg                                                       <  1:1629222/150‑1 (MQ=255)
                                                                                                tctcACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCg                                                       >  2:689105/1‑150 (MQ=255)
                                                                                                 ctcACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGt                                                      >  2:1363135/1‑150 (MQ=255)
                                                                                                   cACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTa                                                    <  2:719684/150‑1 (MQ=255)
                                                                                                   cACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAAATCACACCTATTACCTTCCTCTGCACTGACACATTCGTTa                                                    >  2:86821/1‑150 (MQ=255)
                                                                                                    aCTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCGCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTaa                                                   >  1:1142867/1‑150 (MQ=255)
                                                                                                       ccagccagAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGACTAGCCAACTCAAAATTCACACCTATTACCTTCCTATGCAGTTACACAGTCGTTAAGTc                                                >  1:342665/1‑150 (MQ=255)
                                                                                                        cagccagAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCTCACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCa                                               >  1:477631/1‑150 (MQ=255)
                                                                                                        cagccagAGCCACCAACGCAGGGCTGGAAAGT‑AAAAACCCGCCGCAAAGTCGGGTTTTTTACGTCCAGATTCAGCCCTCCTTTCAAATGACTAGCCAACTCAAAATTCACACCTAGTACCTTCCTCTCCACATCCACAAGCGTTAAggaa                                               >  2:1217197/1‑147 (MQ=255)
                                                                                                                    aataCTCAGGGCTGGAAAGT‑AAAAAACCGACGAAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTg                                   <  1:318123/147‑1 (MQ=255)
                                                                                                                      aaCTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTCTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGAc                                 <  2:1124691/150‑1 (MQ=255)
                                                                                                                       aCTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACt                                <  2:300840/150‑1 (MQ=255)
                                                                                                                          cAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGAATTAt                             >  2:515887/1‑150 (MQ=255)
                                                                                                                          cAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCATCCTCTGCACTTACACAATCGTTAAGTCATATATGTTTTTGACTTaa                             >  2:670647/1‑149 (MQ=255)
                                                                                                                            gggCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATcc                           >  2:422137/1‑150 (MQ=255)
                                                                                                                             ggCTGGAAATT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCg                          <  2:435664/150‑1 (MQ=255)
                                                                                                                             ggCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCTTATATGTTTTTGACTTATCCg                          >  1:1609463/1‑150 (MQ=255)
                                                                                                                               gttgAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCt                        <  2:1158833/147‑1 (MQ=255)
                                                                                                                               cTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCTCTTACACATTCGTTAAGTCATATATGTTTTTGACTTATACGCt                        >  2:1091542/1‑150 (MQ=255)
                                                                                                                               cTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCt                        <  2:1158828/150‑1 (MQ=255)
                                                                                                                                  gAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCg                     <  1:1596268/150‑1 (MQ=255)
                                                                                                                                  gAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCg                     >  2:770587/1‑150 (MQ=255)
                                                                                                                                  gAAAGT‑AAAAAAACGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTACTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCg                     <  2:813279/150‑1 (MQ=255)
                                                                                                                                     cgtaaaaaaaCCGACTTTGCGTCGG‑TTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCGAag                  >  2:1423995/5‑150 (MQ=255)
                                                                                                                                     aGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCGAag                  <  1:1557285/150‑1 (MQ=255)
                                                                                                                                     aGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCGAag                  <  2:1463098/150‑1 (MQ=255)
                                                                                                                                       taaaaaaaCCGACTTTGCGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATTTTTTTGACTTATCCGCTTCGAaga                 >  1:1610767/3‑150 (MQ=255)
                                                                                                                                         aaaaaaCCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCGAagaga               >  1:763333/1‑150 (MQ=255)
                                                                                                                                          aaaaaCCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTCCCTTCCTCTGCACTTACACATTCGTTAAGTCATATAAGTTTTTGACTTATCCGCTTCGAagagag              >  1:503485/1‑150 (MQ=255)
                                                                                                                                          aaaaaCCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCGAagagag              <  1:1202579/150‑1 (MQ=255)
                                                                                                                                           aaaaCCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCGAAgagaga             >  2:1204390/1‑150 (MQ=255)
                                                                                                                                             aaCCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCGAAGAGAGaca           <  2:1479846/150‑1 (MQ=255)
                                                                                                                                               ccGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCATTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCGAAGAGAGACACt         <  2:553199/150‑1 (MQ=255)
                                                                                                                                                cGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCGAAGAGAGACACTa        <  1:205549/150‑1 (MQ=255)
                                                                                                                                                 gACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCGAAGAGAGACACTAc       <  1:1431729/150‑1 (MQ=255)
                                                                                                                                                  aCGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCGAAGAGAGACACTAcc      >  1:655015/1‑150 (MQ=255)
                                                                                                                                                   cGCAAAGTCGGTTTTTTTACGTCCTGATTCAGAACTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTAATACCTTCCACTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCGAAGAGAGACACTACCt     <  2:308388/150‑1 (MQ=255)
                                                                                                                                                      aaaGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCGAAGAGAGACACTACCTGca  <  1:44812/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
TATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCGAAGAGAGACACTACCTGCA  >  NC_000913/3648212‑3648510

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: ‑