Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,648,362 A→G 6.4% noncoding (69/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,3620AG 6.4% 404.9 / 13.9 142noncoding (69/84 nt)agrBantisense sRNA antitoxin for the DinQ toxin
Reads supporting (aligned to +/- strand):  ref base A (71/61);  new base G (6/3);  total (78/64)
Fisher's exact test for biased strand distribution p-value = 5.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.94e-01

Only 100 of 145 total aligned reads displayed.

ATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCGAAGAGAGACACTACCTGCA  >  NC_000913/3648213‑3648510
                                                                                                                                                     |                                                                                                                                                    
        aaGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtt                                                                                                                                              <  2:636328/150‑1 (MQ=255)
             ctccctccGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACTTTGCGTCGGtttttta                                                                                                                                         <  1:959673/150‑2 (MQ=38)
              tccctccGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATATGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTa                                                                                                                                        >  2:230281/1‑150 (MQ=255)
              tccctccGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTa                                                                                                                                        >  1:682895/1‑150 (MQ=255)
               ccctccGCGAGCTGAAGAAATCGCTAAATCTTGCAATGTTAGCCACTTGCTAATAGTATTGAGCTGTTAGATAAGATCTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACTTTGCGTCGGttttttact                                                                                                                                       <  2:309857/150‑4 (MQ=37)
                cctccGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACg                                                                                                                                      >  1:1154217/1‑150 (MQ=255)
                   ccGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTATCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACTTTGCGTCGGttttttactttcc                                                                                                                                   <  1:1650837/150‑8 (MQ=37)
                   ccGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTcc                                                                                                                                   >  2:1278354/1‑150 (MQ=255)
                   ccGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTcc                                                                                                                                   >  2:739381/1‑150 (MQ=255)
                      ccAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGa                                                                                                                                <  1:1648119/150‑1 (MQ=255)
                            gaagaaATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTCAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGa                                                                                                                          <  1:414762/150‑1 (MQ=255)
                               gaaATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGAcct                                                                                                                       <  1:662901/150‑1 (MQ=255)
                               gaaATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTATCTAACTCCAGCCAGAGCCACCAACTCAGGTCTGGAAAGTAAAAAACCGAAGAAAAGACGGTTTTTTTACGTCATGATTAAGAcaa                                                                                                                       >  1:731424/1‑148 (MQ=255)
                                    tgCTAATTCTTGCAATATTATCACCTTGGTAATAGTATTGATGTGTTAGATCATCTCTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCttt                                                                                                                  <  1:937325/149‑1 (MQ=17)
                                    cGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCttt                                                                                                                  >  2:1592219/1‑150 (MQ=255)
                                     gCTAATTCTTGCAATGTTAGCCTCTGTCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCAAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTc                                                                                                                 <  1:1481655/150‑1 (MQ=255)
                                     gCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTc                                                                                                                 <  1:1592219/150‑1 (MQ=255)
                                     gCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTc                                                                                                                 >  2:540339/1‑150 (MQ=255)
                                     gCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTc                                                                                                                 <  2:1679854/150‑1 (MQ=255)
                                       tAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCaa                                                                                                               >  1:1580517/1‑150 (MQ=255)
                                        aaTTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCaaa                                                                                                              <  1:1407170/150‑1 (MQ=255)
                                         tttCTTGCAATCTTATCCAATGTCTAATAGTATTGAGGGGTTAGATCATTCCTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCGGGAAAGTAAAAAACCGACGCAAAGGCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAAt                                                                                                             <  1:1485465/149‑1 (MQ=18)
                                          ttCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATg                                                                                                            <  1:1610151/150‑1 (MQ=255)
                                            cTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGaa                                                                                                          <  1:47656/150‑1 (MQ=255)
                                            cTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGaa                                                                                                          >  1:823757/1‑150 (MQ=255)
                                            cTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGaa                                                                                                          >  2:47656/1‑150 (MQ=255)
                                              tGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATa                                                                                                        >  1:1135118/1‑150 (MQ=255)
                                                cAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGc                                                                                                      >  1:860025/1‑150 (MQ=255)
                                                cAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGc                                                                                                      >  1:1353880/1‑150 (MQ=255)
                                                 actattAGCCAGTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGcc                                                                                                     >  2:1676570/5‑150 (MQ=255)
                                                 aaTGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGcc                                                                                                     >  2:897158/1‑150 (MQ=255)
                                                   tGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTGCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCaa                                                                                                   >  2:1172143/1‑150 (MQ=255)
                                                       aGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCa                                                                                               <  1:1278354/150‑1 (MQ=255)
                                                        gCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCATAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCaa                                                                                              >  1:1129720/1‑150 (MQ=255)
                                                           aCTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAAt                                                                                           <  2:1328715/150‑1 (MQ=255)
                                                               gCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTcac                                                                                       <  2:1353880/150‑1 (MQ=255)
                                                                cTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTcaca                                                                                      <  1:1578906/150‑1 (MQ=255)
                                                                  aaTAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACAACAGAAAGTCGGTTTTTTTTCGTCCAGAATCACTCCTCCGTTCAAAAGTATATCCAACTCAAAATTCCCAcca                                                                                   >  2:1037329/1‑149 (MQ=21)
                                                                       tATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTa                                                                               <  2:217813/150‑1 (MQ=255)
                                                                        aTTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTAc                                                                              <  1:1150567/150‑1 (MQ=255)
                                                                        aTTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTAc                                                                              >  2:691345/1‑150 (MQ=255)
                                                                         ttGAGCTGTTAGATAAGAACTCTCTCACTCCCGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGCCCGCATTTCAAATGATTAGCCAACTCAAAATTCACACCTATTAca                                                                             >  1:286774/1‑149 (MQ=25)
                                                                          tGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAATTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCt                                                                            >  1:395854/1‑150 (MQ=255)
                                                                           gAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCtt                                                                           >  2:1386221/1‑150 (MQ=255)
                                                                             gCTGTTAGATATGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTcc                                                                         <  1:230281/150‑1 (MQ=255)
                                                                             gCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTcc                                                                         >  2:1390728/1‑150 (MQ=255)
                                                                              cTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCct                                                                        >  1:1043033/1‑150 (MQ=255)
                                                                               tGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCctc                                                                       <  1:171976/150‑1 (MQ=255)
                                                                               tGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCctc                                                                       <  1:1174769/150‑1 (MQ=255)
                                                                                   aGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACCCCTATTACCTTCCTCTGCa                                                                   >  1:1670012/1‑150 (MQ=255)
                                                                                    taTAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAATTAAAAAACCGACTAAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCTACTCAAAATTCACACCTATTACCTTCCTCTGCAc                                                                  <  1:287700/149‑1 (MQ=255)
                                                                                    gATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCAc                                                                  >  1:807233/1‑150 (MQ=255)
                                                                                    gATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCAc                                                                  >  1:707381/1‑150 (MQ=255)
                                                                                     aTAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACt                                                                 >  2:790340/1‑150 (MQ=255)
                                                                                     aTAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACt                                                                 <  1:1081918/150‑1 (MQ=255)
                                                                                     aTAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACt                                                                 <  2:860025/150‑1 (MQ=255)
                                                                                        aGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAATAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAATTAATTGCCAACGCAAAATTGACCCCTATTTCCTTACTCTGCACTTtca                                                             >  1:765689/1‑150 (MQ=21)
                                                                                              ctctcACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTc                                                        >  2:947439/1‑150 (MQ=255)
                                                                                              ctctcACTAAAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGCTTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTc                                                        <  1:121492/150‑1 (MQ=255)
                                                                                               tctcACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCg                                                       <  1:763902/150‑1 (MQ=255)
                                                                                                ctcACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGATTAGCCAACTCAAATTTCACACCTATTACCTTCCTCTGCCCTTACACATTCGt                                                      >  2:743073/1‑150 (MQ=255)
                                                                                                ctcACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGt                                                      <  2:1135118/150‑1 (MQ=255)
                                                                                                 tcACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGtt                                                     <  1:897158/150‑1 (MQ=255)
                                                                                                    cTCCAGCCAGAGACACCAACTCAGGGCTGGAAAGAAAAAAACCGACGCAAAGTCGGTTTTTTTACTTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTATGCACTTACACATTCGTTAAg                                                  >  2:619580/1‑150 (MQ=255)
                                                                                                    atCCAGCCAGAGCAACCAACTCAGTGCTGGAAAGTAAAAAACCGACGCAAATTCGGTTTTTTTACGTCCTGATTCAGACGTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAg                                                  <  2:731424/149‑1 (MQ=255)
                                                                                                      ccagccagAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAAGAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTc                                                >  1:140793/1‑150 (MQ=255)
                                                                                                      ccagccagAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCGTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTc                                                >  1:1314083/1‑150 (MQ=255)
                                                                                                       cagccagAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCa                                               >  2:1277415/1‑150 (MQ=255)
                                                                                                         gccagAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTGAGCCCTCCTTTCAAATGAATAGCCAAATCAAAATTCACACCTATTACCTACCGCTGCACTTACACATTGGTTAAGTCata                                             >  2:269553/1‑150 (MQ=37)
                                                                                                         gccagAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCata                                             <  1:1386221/150‑1 (MQ=255)
                                                                                                         gccagAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCata                                             >  1:718183/1‑150 (MQ=255)
                                                                                                           cagAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAACTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACAATCGTAAAGTAatata                                           >  1:12156/1‑150 (MQ=255)
                                                                                                               gCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGtt                                       >  1:1107298/1‑150 (MQ=255)
                                                                                                               gCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTACCTCTGCACTTACACATTCGTTAAGTCAAATATgat                                       >  2:422060/1‑148 (MQ=255)
                                                                                                                     aaCTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGAc                                 <  1:383145/150‑1 (MQ=255)
                                                                                                                      aCTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACt                                <  2:411740/150‑1 (MQ=255)
                                                                                                                       cTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATGCAGACCTCCTTGCAAATGAATAGCCAACTCAAAATTCACCCCTATTACCTTCCTCTGCACTTACCCATTCGTTAAGTCATATATGTTTTTGACtt                               >  2:426917/1‑150 (MQ=255)
                                                                                                                              cTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCAAACCTATTACCTTCCTCAGCAATAACAAATTCGATAAGTCATATATGTTTTTGACTTATCCGCt                        >  1:291480/1‑150 (MQ=38)
                                                                                                                                ggAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTc                      >  2:595347/1‑150 (MQ=255)
                                                                                                                                ggAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTc                      <  2:1154217/150‑1 (MQ=255)
                                                                                                                                 gacgtaaaaaaaCCGACTTTGCGTCGGTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCGa                    >  2:723188/7‑150 (MQ=255)
                                                                                                                                 gAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCg                     <  2:1534102/150‑1 (MQ=255)
                                                                                                                                 gAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCg                     <  2:813654/150‑1 (MQ=255)
                                                                                                                                  aaaGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCGa                    >  1:1059704/1‑150 (MQ=255)
                                                                                                                                   aaGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCGaa                   >  1:986580/1‑150 (MQ=255)
                                                                                                                                    aGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCGAag                  <  2:395854/150‑1 (MQ=255)
                                                                                                                                     gTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGAGATCCTTTCAAATGAATAGCCAACTCAAAATTAACACCTATTTCCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCGAaga                 >  1:841705/1‑150 (MQ=255)
                                                                                                                                     gTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCGAaga                 <  1:830228/150‑1 (MQ=255)
                                                                                                                                     gTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCGAaga                 >  1:94502/1‑150 (MQ=255)
                                                                                                                                      tAAATAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGATTAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCGAagag                <  2:1560071/150‑1 (MQ=255)
                                                                                                                                      tAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCATTCAAATGAATAACAAACTCAAAATTCACACCTATTGCCTTCCTCTGCACTTAAACATTCGTTAAGTCATATATGTTTTTGCCTTATGCGCTTCGAagag                >  2:1175456/1‑150 (MQ=39)
                                                                                                                                      aaaaaaaCCGACTTTGCGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCGAagag                >  2:1139943/2‑150 (MQ=255)
                                                                                                                                       aaaaaaCCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCGAagaga               <  1:678974/150‑1 (MQ=255)
                                                                                                                                       aaaaaaCCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCGAagaga               <  2:101255/150‑1 (MQ=255)
                                                                                                                                             ccGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCGAAGAGAGACACt         >  2:1089211/1‑150 (MQ=255)
                                                                                                                                               gACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCGAAGAGAGACACTAc       >  1:1101113/1‑150 (MQ=255)
                                                                                                                                                aatcAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTAAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCGAAGAGAGACACTAcc      <  2:917005/147‑1 (MQ=255)
                                                                                                                                                aCGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCGAAGAGAGACACTAcc      >  1:572554/1‑150 (MQ=255)
                                                                                                                                                    aaaGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCGAAGAGAGACACTACCTGca  >  2:546931/1‑150 (MQ=255)
                                                                                                                                                    aaaGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCGAAGAGAGACACTACCTGca  >  2:38459/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCGAAGAGAGACACTACCTGCA  >  NC_000913/3648213‑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: