Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,648,363 A→C 5.3% noncoding (70/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,3630AC 5.3% 450.3 / 10.9 149noncoding (70/84 nt)agrBantisense sRNA antitoxin for the DinQ toxin
Reads supporting (aligned to +/- strand):  ref base A (79/62);  new base C (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 = 7.95e-01

Only 100 of 153 total aligned reads displayed.

TATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCGAAGAGAGACACTACCTGCAAC  >  NC_000913/3648214‑3648512
                                                                                                                                                      |                                                                                                                                                     
tataCCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCaaa                                                                                                                                                       >  1:719684/1‑150 (MQ=255)
   aCCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTc                                                                                                                                                    >  2:1354556/1‑150 (MQ=255)
      gAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGt                                                                                                                                                 >  1:1539097/1‑150 (MQ=255)
      gAAGTGCTCCCTCCGCCAGCTCAAGAAATCGCTAATTCATGCAATGTTAGCCACTGTCTAATAGTATTGAGCTGTAAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGTCTGGAAAGTAAAAAACCGACGCAAAGTCGGt                                                                                                                                                 <  1:348758/150‑1 (MQ=37)
           gCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACTTTGCGTCGGtttttt                                                                                                                                            <  1:602636/150‑1 (MQ=38)
           gCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtttttt                                                                                                                                            >  2:1596268/1‑150 (MQ=255)
           gCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtttttt                                                                                                                                            >  2:410738/1‑150 (MQ=255)
            ctccctccGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACTTTGCGTCGGtttttta                                                                                                                                           <  1:1709848/150‑2 (MQ=38)
             tccctccGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAAACCGACTTTGCGTCGGtttttta                                                                                                                                           <  1:1509267/150‑2 (MQ=38)
              ccctccGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTAc                                                                                                                                         >  2:323148/1‑150 (MQ=255)
               cctccGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGCCTGGAAAGTATAAAACCGACGCAAAGTCGGTTTTTTTACg                                                                                                                                        <  1:1352635/150‑1 (MQ=255)
                  ccGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTa                                                                                                                                          >  2:1650537/1‑145 (MQ=255)
                  ccGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTa                                                                                                                                          <  1:1650537/145‑1 (MQ=255)
                   cGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGCCGCAAAGTCGGTTTTTTTACGTCCt                                                                                                                                    <  1:989401/150‑1 (MQ=255)
                   cGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGAAAATAGTATTGATCTGTTAGATAAGAACACTCTCACTCCAGCCAGAGCCACCAACCCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCt                                                                                                                                    <  1:765699/150‑1 (MQ=37)
                    gCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTg                                                                                                                                   >  1:1463098/1‑150 (MQ=255)
                       aGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGAtt                                                                                                                                >  2:620562/1‑150 (MQ=255)
                         cTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCa                                                                                                                              >  2:1428515/1‑150 (MQ=255)
                           gaagaaATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGa                                                                                                                            >  1:1256655/1‑150 (MQ=255)
                             agaaATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGCCGCAAAGTCGGTTTTTTTACGTCAAGATTCAGAcc                                                                                                                          >  2:1609355/1‑150 (MQ=255)
                                    gCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGATCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTc                                                                                                                   <  1:955547/150‑1 (MQ=255)
                                    gCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTc                                                                                                                   <  1:1686944/150‑1 (MQ=255)
                                      aaaTTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAAATCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCaa                                                                                                                 <  1:895982/149‑1 (MQ=255)
                                       aaTTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCaaa                                                                                                                >  1:1037729/1‑150 (MQ=255)
                                        atccTTGCAATCTTAGCCCCAGGCTAATAGTATTGAGCCGTTAGATAAGACCTCTCTCACTCCAGCCAGATCCACCAACTCAGGGCCGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAAt                                                                                                               <  1:1732497/147‑1 (MQ=25)
                                         ttCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATg                                                                                                              <  1:423238/150‑1 (MQ=255)
                                           cTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGaa                                                                                                            >  2:749078/1‑150 (MQ=255)
                                            ttGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAAt                                                                                                           >  1:142172/1‑150 (MQ=255)
                                             tGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATa                                                                                                          >  1:300840/1‑150 (MQ=255)
                                               tactattAGCCAGTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACTGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGc                                                                                                        >  1:1201943/6‑150 (MQ=255)
                                               cAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGc                                                                                                        >  2:738388/1‑150 (MQ=255)
                                                actattAGCCAGTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGcc                                                                                                       >  2:157778/5‑150 (MQ=255)
                                                actattAGCCAGTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATACAGACCTCCTTTCAAATGAATAGcc                                                                                                       >  1:32725/5‑150 (MQ=255)
                                                aaTGTTACCCACTGCCTAATAGTATTGCGCTGTACTATAATAACTATATCACTCCAGCCACAGCCACCACCTCAGGGATGGAATGTAAAAAACCGACGCAAAGTCGGTTTTTTTAGTTCCTGATTCAGACCTCATTTCAAATGAATAccc                                                                                                       <  2:1630686/150‑1 (MQ=18)
                                                   gTTAGCCACTGGCTATTAGTTTTGAGCTGTTAGATAAGAAGTCTCTCACTCCAGACAGAGCCACCAACTCATGTCTTTTAAGTAAAAAACCGACGCATAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAAc                                                                                                    <  1:510449/150‑1 (MQ=255)
                                                   gTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCTAc                                                                                                    <  1:539115/150‑1 (MQ=255)
                                                     tAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTc                                                                                                  >  1:1159206/1‑150 (MQ=255)
                                                           cTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAAtt                                                                                            >  1:1024044/1‑150 (MQ=255)
                                                             ggCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTca                                                                                          >  2:346744/1‑150 (MQ=255)
                                                              gCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTcac                                                                                         <  1:542371/150‑1 (MQ=255)
                                                                 aaTAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCAACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACAcc                                                                                      <  1:1446903/150‑1 (MQ=255)
                                                                     gTATTGATGTGTTAGATGATTGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTAtt                                                                                  <  2:168310/150‑1 (MQ=17)
                                                                       aTTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTAc                                                                                >  2:1202579/1‑150 (MQ=255)
                                                                       aTTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTAc                                                                                >  1:809317/1‑150 (MQ=255)
                                                                       aTTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTAc                                                                                <  1:496230/150‑1 (MQ=255)
                                                                          gAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCtt                                                                             >  1:1402073/1‑150 (MQ=255)
                                                                          gAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCtt                                                                             >  2:1324514/1‑150 (MQ=255)
                                                                           aGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTc                                                                            >  1:1590658/1‑150 (MQ=255)
                                                                              tGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACGAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCctc                                                                         >  2:205549/1‑150 (MQ=255)
                                                                                ttAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTg                                                                       <  2:1321235/150‑1 (MQ=255)
                                                                                  aGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCa                                                                     >  2:44812/1‑150 (MQ=255)
                                                                                   gATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCAc                                                                    <  1:323148/150‑1 (MQ=255)
                                                                                     tAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACtt                                                                  >  2:637726/1‑150 (MQ=255)
                                                                                     tAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACtt                                                                  <  2:767647/150‑1 (MQ=255)
                                                                                     tAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTAAAAATTCACCCCTATTACCTTCCTCTGCAcgt                                                                  >  1:1617004/1‑148 (MQ=255)
                                                                                      aaGAACTCTCTCACTCCAGCCAGAGCCACCAACTTAGGGCTGGATAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCTACTTAAAATTCACACCTATTACCTTCCTCTGCACTTa                                                                 <  1:1158444/150‑1 (MQ=255)
                                                                                      aaGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGCCCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTa                                                                 >  2:871282/1‑150 (MQ=255)
                                                                                       aGAACTCTCTCACTCCAGCCAGTGCCACCAACTCAGGGCTGGAAAGTATAAAACAGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTac                                                                <  1:1609355/150‑1 (MQ=255)
                                                                                          aCTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTAcaca                                                             >  1:201741/1‑150 (MQ=255)
                                                                                           ctctctCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCACCTCAAACTTCCCACCTAGTACCTTCCTGTGCGCGTAGACAt                                                            >  1:301987/1‑150 (MQ=21)
                                                                                           ctctctCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACAt                                                            >  1:373684/1‑150 (MQ=255)
                                                                                             ctctcACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGCCCTCCGTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACCCATTc                                                          >  1:1121064/1‑150 (MQ=255)
                                                                                             ctctcACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCTCTTACACAttt                                                          <  1:620562/150‑2 (MQ=255)
                                                                                              tctcACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCg                                                         <  1:1629222/150‑1 (MQ=255)
                                                                                              tctcACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCg                                                         >  2:689105/1‑150 (MQ=255)
                                                                                               ctcACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGt                                                        >  2:1363135/1‑150 (MQ=255)
                                                                                                 cACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTa                                                      <  2:719684/150‑1 (MQ=255)
                                                                                                 cACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAAATCACACCTATTACCTTCCTCTGCACTGACACATTCGTTa                                                      >  2:86821/1‑150 (MQ=255)
                                                                                                  aCTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCGCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTaa                                                     >  1:1142867/1‑150 (MQ=255)
                                                                                                     ccagccagAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGACTAGCCAACTCAAAATTCACACCTATTACCTTCCTATGCAGTTACACAGTCGTTAAGTc                                                  >  1:342665/1‑150 (MQ=255)
                                                                                                      cagccagAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCTCACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCa                                                 >  1:477631/1‑150 (MQ=255)
                                                                                                        gccagAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCata                                               >  1:719351/1‑150 (MQ=255)
                                                                                                                  ccaACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTg                                     <  2:1715906/150‑1 (MQ=255)
                                                                                                                    aaCTCAGGGCTGGAAAGTAAAAAACCGAGGCAAAGTCGGTTTTTTTACTTCCTGATTCAGACCTCCTTTCAAATTAATATCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGAc                                   <  2:548150/150‑1 (MQ=255)
                                                                                                                     aCTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACt                                  <  2:300840/150‑1 (MQ=255)
                                                                                                                     aCTCAGGGATGGAAAGTAAAAAACCGACTTTGCGTCGGTTTTTTTGCGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACt                                  <  2:1081897/150‑1 (MQ=255)
                                                                                                                        cAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGAATTAt                               >  2:515887/1‑150 (MQ=255)
                                                                                                                        cAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCATCCTCTGCACTTACACAATCGTTAAGTCATATATGTTTTTGACTTaa                               >  2:670647/1‑149 (MQ=255)
                                                                                                                          gggCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATcc                             >  2:422137/1‑150 (MQ=255)
                                                                                                                           ggCTGGAAATTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCg                            <  2:435664/150‑1 (MQ=255)
                                                                                                                           ggCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCTTATATGTTTTTGACTTATCCg                            >  1:1609463/1‑150 (MQ=255)
                                                                                                                             cTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCTCTTACACATTCGTTAAGTCATATATGTTTTTGACTTATACGCt                          >  2:1091542/1‑150 (MQ=255)
                                                                                                                             cTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCt                          <  2:1158828/150‑1 (MQ=255)
                                                                                                                                gAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCg                       >  2:770587/1‑150 (MQ=255)
                                                                                                                                gAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCg                       <  2:589529/150‑1 (MQ=255)
                                                                                                                                gAAAGTAAAAAAACGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTACTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCg                       <  2:813279/150‑1 (MQ=255)
                                                                                                                                   cgtaaaaaaaCCGACTTTGCGTCGGTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCGAag                    >  2:1423995/5‑150 (MQ=255)
                                                                                                                                   aGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCGAag                    <  2:1463098/150‑1 (MQ=255)
                                                                                                                                    gtaaaaaaaCCGACTTTGCGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACGTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCGaat                    >  2:1152644/4‑149 (MQ=255)
                                                                                                                                     tAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCAGCTTCGAagag                  >  2:1668205/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:157778/150‑1 (MQ=255)
                                                                                                                                                 cGCAAAGTCGGTTTTTTTACGTCCTGATTCAGAACTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTAATACCTTCCACTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCGAAGAGAGACACTACCt       <  2:308388/150‑1 (MQ=255)
                                                                                                                                                    aaaGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCGAAGAGAGACACTACCTGca    <  1:44812/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
TATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCGAAGAGAGACACTACCTGCAAC  >  NC_000913/3648214‑3648512

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: