Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,648,363 A→C 6.1% 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 6.1% 428.9 / 12.3 142noncoding (70/84 nt)agrBantisense sRNA antitoxin for the DinQ toxin
Reads supporting (aligned to +/- strand):  ref base A (72/60);  new base C (6/3);  total (78/64)
Fisher's exact test for biased strand distribution p-value = 7.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.08e-01

Only 100 of 144 total aligned reads displayed.

AAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCGAAGAGAGACACTACCTGCA  >  NC_000913/3648221‑3648510
                                                                                                                                              |                                                                                                                                                   
aaGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtt                                                                                                                                              <  2:636328/150‑1 (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:739381/1‑150 (MQ=255)
              ccAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGa                                                                                                                                <  1:1648119/150‑1 (MQ=255)
                    gaagaaATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTCAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGa                                                                                                                          <  1:414762/150‑1 (MQ=255)
                      agaaATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGAcc                                                                                                                        >  2:1719459/1‑150 (MQ=255)
                       gaaATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTATCTAACTCCAGCCAGAGCCACCAACTCAGGTCTGGAAAGTAAAAAACCGAAGAAAAGACGGTTTTTTTACGTCATGATTAAGAcaa                                                                                                                       >  1:731424/1‑148 (MQ=255)
                            tgCTAATTCTTGCAATATTATCACCTTGGTAATAGTATTGATGTGTTAGATCATCTCTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCttt                                                                                                                  <  1:937325/149‑1 (MQ=17)
                             gCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGTCGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTc                                                                                                                 >  2:1036031/1‑150 (MQ=255)
                             gCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTc                                                                                                                 <  1:1592219/150‑1 (MQ=255)
                             gCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTc                                                                                                                 <  1:665794/150‑1 (MQ=255)
                             gCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTc                                                                                                                 >  2:540339/1‑150 (MQ=255)
                              cTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCAATCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGTTTTTTGTACGTCCTGAGTCAGACCTCCTTTCa                                                                                                                >  1:253187/1‑150 (MQ=255)
                               tAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCaa                                                                                                               >  1:1580517/1‑150 (MQ=255)
                                aaTTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCaaa                                                                                                              <  1:1407170/150‑1 (MQ=255)
                                 tttCTTGCAATCTTATCCAATGTCTAATAGTATTGAGGGGTTAGATCATTCCTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCGGGAAAGTAAAAAACCGACGCAAAGGCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAAt                                                                                                             <  1:1485465/149‑1 (MQ=18)
                                 aTTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGTCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAAt                                                                                                             <  1:1157972/150‑1 (MQ=255)
                                  ttCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATg                                                                                                            <  1:1610151/150‑1 (MQ=255)
                                    cTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGaa                                                                                                          >  2:47656/1‑150 (MQ=255)
                                      tGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATa                                                                                                        >  1:1135118/1‑150 (MQ=255)
                                        cAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGc                                                                                                      >  1:1353880/1‑150 (MQ=255)
                                        cAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGc                                                                                                      >  1:860025/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)
                                              tAGCCACTGGATAATAGTTTAGAGCTGAAAGATAAGAACTCTCTTACTCCAGCCAGATCCAACAACTCAGGGTTGGAAAGTTAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTc                                                                                                <  2:269391/150‑1 (MQ=25)
                                               aGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCa                                                                                               <  1:1278354/150‑1 (MQ=255)
                                               aGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCa                                                                                               >  1:411740/1‑150 (MQ=255)
                                                gCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCATAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCaa                                                                                              >  1:1129720/1‑150 (MQ=255)
                                                   aCTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAAt                                                                                           <  2:1328715/150‑1 (MQ=255)
                                                   aCTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAAt                                                                                           <  2:1289487/150‑1 (MQ=255)
                                                      gtcTATTAGTATTGAGCTGTTAGATTAGTTCTCTCTCACTCCAGGGAGAGCCACCAACTCAGGGGTTGAAAGTAAAAAACCGACGCAATGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAATATTca                                                                                        <  1:754132/148‑1 (MQ=21)
                                                          aaTAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACAACAGAAAGTCGGTTTTTTTTCGTCCAGAATCACTCCTCCGTTCAAAAGTATATCCAACTCAAAATTCCCAcca                                                                                   >  2:1037329/1‑149 (MQ=21)
                                                            tAGTATTGGGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTTGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTa                                                                                  <  1:361696/150‑1 (MQ=255)
                                                               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)
                                                                     gCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTcc                                                                         <  1:691345/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)
                                                                          tAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGc                                                                    >  1:805880/1‑150 (MQ=255)
                                                                           aGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACCCCTATTACCTTCCTCTGCa                                                                   >  1:1670012/1‑150 (MQ=255)
                                                                            gATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCAc                                                                  >  1:807233/1‑150 (MQ=255)
                                                                            gATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCAc                                                                  >  1:707381/1‑150 (MQ=255)
                                                                             aTAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACt                                                                 >  2:790340/1‑150 (MQ=255)
                                                                             aTAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACt                                                                 <  2:860025/150‑1 (MQ=255)
                                                                             aTAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACt                                                                 <  1:1081918/150‑1 (MQ=255)
                                                                                    ctctctCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACAt                                                          >  1:830165/1‑150 (MQ=255)
                                                                                      ctctcACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTc                                                        >  2:947439/1‑150 (MQ=255)
                                                                                      ctctcACTAAAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGCTTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTc                                                        <  1:121492/150‑1 (MQ=255)
                                                                                       tctcACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCg                                                       <  1:763902/150‑1 (MQ=255)
                                                                                       tctcACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTATCAAATGAATAGCCAACTCAAAATTAACACCTATTACCTTCCTCTGCACTTACAAATTCg                                                       >  1:1492053/1‑150 (MQ=255)
                                                                                        ctcACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGATTAGCCAACTCAAATTTCACACCTATTACCTTCCTCTGCCCTTACACATTCGt                                                      >  2:743073/1‑150 (MQ=255)
                                                                                        ctcACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGt                                                      <  2:1135118/150‑1 (MQ=255)
                                                                                            cTCCAGCCAGAGACACCAACTCAGGGCTGGAAAGAAAAAAACCGACGCAAAGTCGGTTTTTTTACTTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTATGCACTTACACATTCGTTAAg                                                  >  2:619580/1‑150 (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)
                                                                                                   cagAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAACTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACAATCGTAAAGTAatata                                           >  1:12156/1‑150 (MQ=255)
                                                                                                       gCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGtt                                       >  1:1107298/1‑150 (MQ=255)
                                                                                                       gCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTACCTCTGCACTTACACATTCGTTAAGTCAAATATgat                                       >  2:422060/1‑148 (MQ=255)
                                                                                                            caACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTGCGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGa                                  <  1:1430552/150‑1 (MQ=255)
                                                                                                             aaCTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGAc                                 >  2:1090715/1‑150 (MQ=255)
                                                                                                              aCTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACt                                <  2:411740/150‑1 (MQ=255)
                                                                                                                  agtgCTGGAAAGTAAAAAAGCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATc                            <  2:12156/147‑1 (MQ=255)
                                                                                                                   gggCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATcc                           <  1:967088/150‑1 (MQ=255)
                                                                                                                      cTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCAAACCTATTACCTTCCTCAGCAATAACAAATTCGATAAGTCATATATGTTTTTGACTTATCCGCt                        >  1:291480/1‑150 (MQ=38)
                                                                                                                      cTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTATCAAATGAATAGCCAACTCAAAATTCACACCTATTACATTCCTCTCCACTTACACAGTCGTTAAGTCATATATGTTTTTGACTTATCCGCt                        >  2:73214/1‑150 (MQ=255)
                                                                                                                        ggAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTc                      <  2:1154217/150‑1 (MQ=255)
                                                                                                                        ggAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTc                      >  2:595347/1‑150 (MQ=255)
                                                                                                                         gacgtaaaaaaaCCGACTTTGCGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCg                     >  2:1357214/7‑150 (MQ=255)
                                                                                                                         gacgtaaaaaaaCCGACTTTGCGTCGGTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCGa                    >  2:723188/7‑150 (MQ=255)
                                                                                                                         gAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCg                     >  2:1642389/1‑150 (MQ=255)
                                                                                                                         gAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCg                     <  2:813654/150‑1 (MQ=255)
                                                                                                                         gAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCg                     <  2:1534102/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)
                                                                                                                             gTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCGAaga                 <  1:830228/150‑1 (MQ=255)
                                                                                                                             gTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCGAaga                 >  1:94502/1‑150 (MQ=255)
                                                                                                                              tAAATAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGATTAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCGAagag                <  2:1560071/150‑1 (MQ=255)
                                                                                                                              aaaaaaaCCGACTTTGCGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCGAagag                >  2:441172/2‑150 (MQ=255)
                                                                                                                              aaaaaaaCCGACTTTGCGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCGAagag                >  2:1139943/2‑150 (MQ=255)
                                                                                                                               aaaaaaCCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCGAagaga               <  1:678974/150‑1 (MQ=255)
                                                                                                                               aaaaaaCCGACGCAAAGTCGGTTTTGTGACGTCCTGAGTCAGCCCTCCTGTCAAATGAATAGCCCCCGCAAACTTCACACCTATTACCTTCCGCTGCACTTACACCTTCGGTAAGGCAGATATGTTTTTGACGTCTCCGCTGCGAagagc               >  1:938206/1‑149 (MQ=255)
                                                                                                                                   aaCCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCGAAGAGAGaca           <  1:540339/150‑1 (MQ=255)
                                                                                                                                      cGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACGTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCGAAGAGAGACACTa        >  1:855455/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)
                                                                                                                                              |                                                                                                                                                   
AAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATATATGTTTTTGACTTATCCGCTTCGAAGAGAGACACTACCTGCA  >  NC_000913/3648221‑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: