Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,199,681 A→T 5.5% *376R (TGA→AGA)  intE ← e14 prophage; putative integrase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,199,6810AT 5.5% 261.3 / 2.4 93*376R (TGA→AGA) intEe14 prophage; putative integrase
Reads supporting (aligned to +/- strand):  ref base A (46/41);  new base T (2/4);  total (48/45)
Fisher's exact test for biased strand distribution p-value = 4.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

TTTCATCGATTTTAGGTGATCTTCTCTGTGATATTTCACGTCTTACCAGTAACTGATATGGCTGTCCGCCGCTCGCTTAAAGTGGACTTTTTAGTTTTTATCATGTGCGGTGAGAAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTTGTGTGTCTAACCCCTGTTCTTTATACAGACGTTCAGATAAAGACCTTTGCTCATGGAATGTCGCA  >  NC_000913/1199533‑1199829
                                                                                                                                                    |                                                                                                                                                    
tttCATCGATTTTAGGTGATCTTCTCTGTGATATTTCACGTCTTACCAGTAACTGATATGGCTGTCCGCCGCTCGCTTAAAGTGGACTTTTTAGTTTTTATCATGTGCGGTGAGAAATTCAATGTGGCGTTGAGATGCTTAAAGGTTcac                                                                                                                                                     <  1:919780/150‑1 (MQ=255)
   cATCGATTTTAGGTGATCTTCTCTGTGATATTTCACGTCTTACCAGTAACTGATATGGCTGTCCGCCGCTCGCTTAAAGTGGACTTTTTAGTTTTTATCATGTGCGGTGAGAAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAAc                                                                                                                                                  <  2:1560480/150‑1 (MQ=255)
   cATCGATTTTAGGTGATCTTCTCTGTGATATTTCACGTCTTACCAGTAACTGATATGGCTGTCCGCCGCTCGCTTAAAGTGGACTTTTTAGTTTTTATCATGTGCGGTGAGAAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAAc                                                                                                                                                  >  2:408660/1‑150 (MQ=255)
    aTCGATTTTAGGTGATCTTCTCTGTGATATTTCACGTCTTACCAGTAACTGATATGGCTGTCCGCCGCTCGCTTAAAGTGGACTTTTTAGTTTTTATCATGTGCGGTGAGAAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACg                                                                                                                                                 <  1:535901/150‑1 (MQ=255)
      cGATTTTAGGTGATCTTCTCTGTGATATTTCACGTCTTACCAGTAACTGATATGGCTGTCCGCCGCTCGCTTAAAGTGGACTTTTTAGTTTTTATCATGTGCGGTGAGAAATTCAATGTAGCGTTGTGAACCTTTAAGCATCTCAACGcc                                                                                                                                               >  2:447887/1‑149 (MQ=255)
       gATTTTAGGTGATCTTCTCTGTGATATTTCACGTCTTACCAGTAACTGATATGGCTGTCCGCCGCTCGCTTAAAGTGGACTTTTTAGTTTTTATCATGTGCGGTGAGAAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTa                                                                                                                                              >  2:707569/1‑150 (MQ=255)
        aTTTTAGGTGATCTTCTCTGTGATATTTCACGTCTTACCAGTAACTGATATGGCTGTCCGCCGCTCGCTTAAAGTGGACTTTTTAGTTTTTATCATGTGCGGTGAGAAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTAc                                                                                                                                             >  1:487834/1‑150 (MQ=255)
            tAGGTGATCTTCTCTGTGATATTTCACGTCTTACCAGTAACTGATATGGCTGTCCGCCGCTCGCTTAAAGTGGACTTTTTAGTTTATATCATGTGCGGTGAGAAATTCAATGTGGAGTTGAGATGCTTAAAGGATCACAACGCTACtttt                                                                                                                                         >  1:1039011/1‑149 (MQ=255)
              ggTGCTCTTCTCTGTGATATTTCACGTCTTACCAGTAACTGATATGGCTGGCCGCCGCTCACGTAAAGTGGACTTTTTAGTGTTTATCATGTGCGGTGAGATATTCAATGTGGCGGTGAGTTGCTTTAAGGTTCACAACGCTACTTTGCt                                                                                                                                       >  2:281114/1‑150 (MQ=255)
              ggTGATCTTCTCTGTGATATTTCACGTCTTACCAGTAACTGATATGGCTGTCCGCCGCTCGCTTAAAGTGGACTTTTTAGTTTTTATCATGTGCGGTGAGAAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAATGCTACTTTGCt                                                                                                                                       >  2:839941/1‑150 (MQ=255)
              ggTGATCTTCTCTGTGATATTTCACGTCTTACCAGTAACTGATATGGCTGTCCGCCGCTCGCTTAAAGTGGACTTTTTAGTTTTTATCATGTGCGGTGAGAAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCt                                                                                                                                       >  1:1460086/1‑150 (MQ=255)
               gTGATCTTCTCTGTGATATTTCACGTCTTACCAGTAACTGATATGGCTGTCCGCCGCTCGCTTAAAGTGGACTTTTTAGTTTTTATCATGTGCGGTGAGAAATTCAATGTAGCGTTGTGAACCTTTAAGCATCTCAACGCCACATTgaat                                                                                                                                      <  1:447887/150‑4 (MQ=37)
               gTGATCTTCTCTGTGATATTTCACGTCTTACCAGTAACTGATATGGCTGTCCGCCGCTCGCTTAAAGTGGACTTTTTAGTTTTTATCATGTGCAGTGAGAAATTCAATGTAGCGTTGTGAACCTTTAAGCATCTCAACGCCACATTgaat                                                                                                                                      <  1:363575/150‑4 (MQ=37)
                tGATCTTCTCTGTGATATTTCACGTCTTACCAGTAACTGATATGGCTGGCCGCCGCTCGCTTAAAGTCGACTTTTTAGTTTTTATCATGTGCGGTGAGAAAATCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTcc                                                                                                                                     <  1:918476/150‑1 (MQ=255)
                 gATCTTCTCTGTGATATTTCACGTCTTACCAGTAACTGATATGGCTGTCCGCCGCTCGCTTAAAGTGGACTTTTTAGTTTTTATCATGTGCGGTGAGAAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTcaa                                                                                                                                    >  1:1035553/1‑148 (MQ=255)
                   tcttctCTGTGATATTTCACGTCTTACCAGTAACTGATATGGCTGTCCGCCGCTCGCTTAAAGTGGACTTTTTAGTTTTTATCATGTGCGGTGAGAAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATc                                                                                                                                  >  1:1338781/1‑150 (MQ=255)
                       ctctGTGATATTTCACGTCTTACCAGTAACTGATATGGCTGTCCGCCGCTCGCTTAAAGTGGACTTTTTAGTTTTTATCATGTGCGGTGAGAAATTCTATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCttt                                                                                                                              >  2:736294/1‑150 (MQ=255)
                         ctGTGATATTTCACGTCTTACCAGTAACTGATATGGCTGTCCGCCGCTCGCTTAAAGTGGACTTTTTAGTTTTTATCATGTGCGGTGAGAAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTAc                                                                                                                            <  1:1109560/150‑1 (MQ=255)
                          tgtgATATTTCACGTCTTACCAGTAACTGATATGGCTGTCCGCCGCTCGCTTAAAGTGGACTTTTTAGTTTTTATCATGTGCGGTGAGAAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTAcc                                                                                                                           >  1:720913/1‑150 (MQ=255)
                             gATATTTCACGTCTTACCAGTAACTGATATGTCTGTCCGCCGCTCGCTTAAAGTGGACTTTTTAGTTTTTATCATGTGCGGTGAGAAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCg                                                                                                                        >  2:1295045/1‑150 (MQ=255)
                             gATATTTCACGTCTTACCAGTAACTGATATGGCTGTCCGCCGCTCGCTTAAAGTGGACTTTTTAGTTTTTATCATGTGCGGTGAGAAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCg                                                                                                                        >  1:976891/1‑150 (MQ=255)
                              atatTTCACGTCTTACCAGAAAATGATATGGCTGTCCGACGCTCGCTTAAAGTGGACTTTTTAGTTTTTATCATGTGCGGAGAGAAATTCAATGTAGCTTTAAGATGCTTAAAGGTTCACAACGCTACTTTGATCCATCCTTTAcataga                                                                                                                       >  2:1518629/1‑145 (MQ=255)
                                atTTCACGTCTTACCAGTAACTGATATGGCTGTCCGCCGCTCGCTTAAAGTGGACTTTTTAGTTTTTATCATGTGCGGTGAGAAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGatc                                                                                                                     <  2:358916/150‑1 (MQ=255)
                                 tttCACGTCTTACCAGTAACTGATATGGCTGTCCGCCGCTCGCTTAAAGTGGACTTTTTAGTTTTTATCATGTGCGGTGAGAAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGatca                                                                                                                    >  1:1443106/1‑150 (MQ=255)
                                  ttCACGTCTTACCAGTAACTGATATGGCTGTCCGCCGCTCGCTTAAAGTGGACTTTTTAGTTTTTATCATGTGCGGTGAGAAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGatcat                                                                                                                   >  2:808804/1‑150 (MQ=255)
                                     aCGTCTTACCAGTAACTGATATGGCTGTCCGCCGCTCGCTTAAAGTGGACTTTTTAGTTTTTATCATGTGCGGTGAGAAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATcatcat                                                                                                                <  2:1338781/150‑1 (MQ=255)
                                        tCTTACCAGTAACTGATATGGCTGTCCGCCGCTCGCTTAAAGTTGACTTTTTAGTTTTTATCATGTGCGGTGAGAAATTCTATGTGGCGTTGAGATGCTTAAATGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGAt                                                                                                             <  1:1515885/150‑1 (MQ=255)
                                           tACCAGTAACTGATATGGCTGTCCGCCGCTCGCTTAAAGTGGACTTTTTAGTTTTTATCATGTGCGGTGAGAAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAAc                                                                                                          <  2:376438/150‑1 (MQ=255)
                                               aGTAACTGATATGGCTGTCCGCCGCTCGCTTAAAGTGGACTTTTTAGTTTTTATCATGTGCGGTGAGAAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATc                                                                                                      >  2:790121/1‑150 (MQ=255)
                                                   aCTGATATGGCTGTCCGCCGCTCGCTTAAAGTGGACTTTTTAGTTTTTATCATGTGCGGTGAGAAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGtt                                                                                                  <  2:407725/150‑1 (MQ=255)
                                                   aCTGATATGGCTGTCCGACGCTCGCGTCAAGTGGACTTTTTAGTTTTTATCATGTGCGGTGAGAAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGtt                                                                                                  <  1:1136804/150‑1 (MQ=255)
                                                    cTGATATGGCTGTCCGCCGATCGCTTAAAGTGGACTTTTTAGTTTTTATCATGTGCGGTGAGAAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGATACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGttt                                                                                                 <  1:1687847/150‑1 (MQ=255)
                                                     tGATATGGCTGTCCGCCGCTCGCTTAAAGTGGACTTTTTAGTTTTTATCATGTTCGGTGTGAAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGGTCAAGCCTTTACCTCGATCATCATGATAACGATCGGTttg                                                                                                >  1:1126210/1‑150 (MQ=255)
                                                     tGATATGGCTGTCCGCCGCTCGCTTAAAGTGGACTTTTTAGTTTTTATCATGTGCGGTGAGAAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTttg                                                                                                <  2:1541827/150‑1 (MQ=255)
                                                      gATATGGCTGTCCGCCGCTCGCTTAAAGTGGACTTTTTAGTTTTTATCATGTGCGGTGAGAAATTCAATGTGGCGTTCAGATGCTTAAAGTTTCACAACTCTACTTTGCTGCCTCCTTTACATCGATCAACATGATAACGATCGGTttgt                                                                                               >  2:493242/1‑150 (MQ=255)
                                                       atatGGCTGTCCGCCGCTCGCTTAAAGTGGACTTTTTAGTTTTTATCATGTGCGGTGAGAAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTttgtt                                                                                              <  2:976891/150‑1 (MQ=255)
                                                       atatGGCTGTCCGCCGCTCGCTTAAAGTGGACTTTTTAGTTTTTATCATGTGCGGTGAGAAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTttgtt                                                                                              >  1:1709978/1‑150 (MQ=255)
                                                        tatGGCTGTCCGCCGCTCGCTTAAAGTGGACTTTTTAGTTTTTATCATGTGCGGTGAGAAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTttgttg                                                                                             >  2:126915/1‑150 (MQ=255)
                                                         atGGCTGTCCGCCGCTCGCTTAAAGTGGACTTTTTAGTTTTTATCATGTGCGGTGAGAAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTttgttgt                                                                                            <  1:1403196/150‑1 (MQ=255)
                                                          tGGCTGTCCGCCGCTCGCTTAAAGTGGACTTTTTAGTTTTTATCATGTGCGGTGAGAAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTttgttgtt                                                                                           <  2:465736/150‑1 (MQ=255)
                                                           ggCTGTCGGCCGCTCGCTTATAGTGGACTTTTTAGTTTTTATCATGTGCGGTGAGAAATTCATTGTGGCGTTGAGATGCTTAAAGGTTCAGAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTttgttgttg                                                                                          <  1:1713422/150‑1 (MQ=255)
                                                               gTCCGCCGCTCGCTTAAAGTGGACTTTTTAGTTTTTATCATGTGCGGTGAGAAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGtt                                                                                      >  1:400118/1‑150 (MQ=255)
                                                                 ccgccgCTCGCTTAAAGTGGACTTTTTAGTTTTTATCATGTGCGGTGAGCAATTCAATGTGGCGGTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGtttt                                                                                    >  1:236597/1‑150 (MQ=255)
                                                                 ccgccgCTCGCTTAAAGTGGACTTTTTAGTTTTTATCATGTGCGGTGAGAAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGtttt                                                                                    <  1:419827/150‑1 (MQ=255)
                                                                 ccgccgCTCGCTTAAAGTGGACTTTTTAGTTTTTATCATGTGCGGTGAGAAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGtttt                                                                                    >  1:653394/1‑150 (MQ=255)
                                                                  cgccgcTCGCTTAAAGTGGACTTTTTAGTTTTTATCATGTGCGGTGAGAAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTa                                                                                   <  1:707569/150‑1 (MQ=255)
                                                                   gccgcTCGCTTATAGTGGACTTTTTAGTTTTTATCATGTGCGGTGAGAAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTAt                                                                                  <  1:1684125/150‑1 (MQ=255)
                                                                   gccgcTCGCTTAAAGTGGCCTTTTTAGTTTTTATCATGTGCGGTGAGAAATTCAAAGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTAt                                                                                  >  2:888873/1‑150 (MQ=255)
                                                                   gccgcTCGCTTAAAGTGGACTTTTTAGTTTTTATCATGTGCGGTGAGAAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTAt                                                                                  <  1:1013487/150‑1 (MQ=255)
                                                                    ccgcTCGCTTAAAGTGGACTTTTTAGTTTTTATCATGTGCGGTGAGAAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTATg                                                                                 <  2:1460086/150‑1 (MQ=255)
                                                                    ccgcTCGCTTAAAGTGGACTTTTTAGTTTTTATCATGTGCGGTGAGAAATTCAATGTGGCGTTGAGATGCTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTATGa                                                                                >  1:555129/1‑150 (MQ=255)
                                                                    ccgcTCGCTTAAACGGGACCTTTTAGTTTTTATCATGCGCGGTGAGAAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTAATTTGCTCCATCCTATACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTATg                                                                                 <  2:1054947/150‑1 (MQ=255)
                                                                    ccgcTCGATTAAAGTGGACTTTTTAGTTTTTATCATGTGCGGTGAGAAATTCAATGTGGCGTTGATATGCTTAAAGGTTCAAAAAGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTATg                                                                                 <  2:310047/150‑1 (MQ=255)
                                                                       ctcgctTAAAGTGGACTTTTTAGTTTTTATCATGTGCGGTGAGAAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTATGAcc                                                                              >  2:1350771/1‑150 (MQ=255)
                                                                       ctcgctTAAAGTGGACTTTTTAGTTTTTATCATGTGCGGTGAGAAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTAAATCCTATACCTCGATCATCATGATAACGATCGGTTTGGTGGTGGTTTTTATGAca                                                                              >  2:301791/1‑149 (MQ=255)
                                                                               aaGTGGACTTTTTAGTTTTTATCATGTGCGGTGAGAAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGtt                                                                      >  2:1387136/1‑150 (MQ=255)
                                                                                aGTGGACTTTTTAGTTTTTATCATGTGCGGTGAGAAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGttt                                                                     <  2:1443106/150‑1 (MQ=255)
                                                                                  tGGACTTTTTAGTTTTTATCATGTGCGGTGAGAAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGGCTGTTGTTGTGTTTTATGACCAAGTAGttttt                                                                   <  2:1059213/150‑1 (MQ=255)
                                                                                   tgACTTTTTAGTTTTTTTCATGTGCGGTGAGAAATTGAATGTTTCGTTTAGATGCTTTATGGTTCTCTACGGTACTTTGCTCCATCCTTTACCTCGATGATCATGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTtg                                                                  <  1:307468/149‑1 (MQ=255)
                                                                                   gtaCTTTTTAGTTTATATCATGTGCGGTGAGAAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACTACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTtg                                                                  <  1:90376/148‑1 (MQ=255)
                                                                                   ggACTTTTTAGTTTTTATCATGTGCGGTGAGAAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTtg                                                                  >  1:20570/1‑150 (MQ=255)
                                                                                     tcTTTTTAGTTTTTATCATGTTCTGTGAGAAATTCATTGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCCTGAGAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTtgtg                                                                <  2:815406/149‑1 (MQ=255)
                                                                                       tttttAGTTTTTATCATGTGCGGTGAGAAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTtgtgtg                                                              <  1:408660/150‑1 (MQ=255)
                                                                                       tttttAGTTTTTATCATGTGCGGTGAGAAATTCAATGTGGCGTTGAGATGATTAAAGGTTCACAACGCTACTTTGCTACATCCTTTAAATCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTtgtgtt                                                              >  1:1280074/1‑149 (MQ=255)
                                                                                             gTTTTTATCATGTGCGGTGAGAAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACATCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTTGTGTGTCTAAc                                                        >  1:352082/1‑150 (MQ=255)
                                                                                              tttttATCATGTGCGGTGAGAAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTATGACAAAGTAGTTTTTGTGTGTCTAAcc                                                       >  2:607300/1‑150 (MQ=255)
                                                                                                  tATCATGTGCGGTGAGAAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTTGTGTGTCTAACCCCTg                                                   >  1:1304698/1‑150 (MQ=255)
                                                                                                   aTCATGTGCGGTGAGAAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTTGTGTGTCTAACCCCTGt                                                  >  2:1511404/1‑150 (MQ=255)
                                                                                                      aTGTGCGGTGAGAAATTCTATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTTGTGTGTCTAACCCCTGTTCt                                               >  2:939742/1‑150 (MQ=255)
                                                                                                      aTGTGCGGTGAGAAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTTGTGTGTCTAACCCCTGTTCt                                               <  1:1387136/150‑1 (MQ=255)
                                                                                                         tgCGGTGAGAAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTTGTGTGTCTAACCCCTGTTCTTta                                            <  2:400118/150‑1 (MQ=255)
                                                                                                          gCGGTGAGAAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTTGTGTGTCTAACCCCTGTTCTTtat                                           >  1:182544/1‑150 (MQ=255)
                                                                                                          gCGGTGAGAAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTTGTGTGTCTAACCCCTGTTCTTtat                                           <  2:487834/150‑1 (MQ=255)
                                                                                                              tGAGAAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTTGTGTGTCTAACCCCTGTTCTTTATACAg                                       <  1:808804/150‑1 (MQ=255)
                                                                                                               gagaAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTTGTGTGTCTAACCCCTGTTCTTtata                                          >  1:809735/1‑146 (MQ=255)
                                                                                                               gagaAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTTGTGTGTCTAACCCCTGTTCTTtata                                          <  2:809735/146‑1 (MQ=255)
                                                                                                               gagaAATTCAATGCGGCGTTGAGATGCTTAAAGGTGCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGAAAACGATCGGTTTGTTGGTGTGTTTTATGACCAAGTAGTTTTTGTGTGTCTAACCCCTGTTCTTTATACAGa                                      <  2:1582208/150‑1 (MQ=255)
                                                                                                                 gaAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTTGTGTGTCTAACCCCTGTTCTTTATACAGACg                                    <  1:1350771/150‑1 (MQ=255)
                                                                                                                   aaTTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTTGTGTGTCTAACCCCTGTTCTTTATACAGACGtt                                  >  2:85521/1‑150 (MQ=255)
                                                                                                                   aaTTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCATCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTTGTGTGTCTAACCCCTGTTCTTTATACAGACGTTc                                 <  2:1537426/150‑1 (MQ=255)
                                                                                                                     ttCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTTGTGTGTCTAACCCCTGTTCTTTATACAGACGTTCa                                <  2:653394/150‑1 (MQ=255)
                                                                                                                     ttCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTTGTGTGTCTAACCCCTGTTCTTTATACAGACGTTCa                                <  2:1135547/150‑1 (MQ=255)
                                                                                                                         aagtagCGTTGTGAACCTTTAAGCATCTCAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTTGTGTGTCTAACCCCTGTTCTTTATACAGACGTTCAGATa                            >  1:1056884/6‑150 (MQ=255)
                                                                                                                         aTGTGGCGTTGAGATGCTTAAAGGTTCTCAACGCTACTTTGCTCCATCCTTTACCTCGATCATCCTGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTGTTTGTGTGTCTAACCCCTGTTCTTTATACAGACGTTCAGATa                            <  2:1039011/150‑1 (MQ=255)
                                                                                                                           gtgGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATACTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTTGTGTGTCTAACCCCTGTTCTTTATACAGACGTTCAGATaaa                          >  1:137342/1‑150 (MQ=255)
                                                                                                                                     gaTGCTTAAAGGTTCACAACGCTACTTTTCTCCATCCTTTTCCTCGATCATCATGATTACGATCGGTGTGTTGTTGTGTTTTATGACCAAGTAGTTTTTGTGTGTCTAACCCCTGTTCTTTATACAGACGTTCAGATAAAGACATTTGCt                <  2:1719564/150‑1 (MQ=255)
                                                                                                                                     gaTGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTTGTGTGTCTAACCCCTGTTCTTTATACAGACGTTCAGATAAAGACCTTTGCt                >  2:1514472/1‑150 (MQ=255)
                                                                                                                                     gaTGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACAATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTTGTGTGTCTAACCCCTGTTCTTTATACAGACGTTCAGATAAAGACCTTTGCt                >  2:1214516/1‑150 (MQ=255)
                                                                                                                                        gCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTTGTGTGTCTAACCCCTGTTCTTTATACAGACGTTCAGATAAAGACCTTTGCTCAt             >  1:935848/1‑150 (MQ=255)
                                                                                                                                         cTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTTGTGTGTCTAACCCCTGTTCTTTATACAGACGTTCAGATAAAGACCTTTGCTCATg            >  2:924002/1‑150 (MQ=255)
                                                                                                                                           tAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTTGTGTGTCTAACCCCTGTTCTTTATACAGACGTTCAGATAAAGACCTTTGCTCATGGa          <  1:1295045/150‑1 (MQ=255)
                                                                                                                                             aaGGTTCACTACGCTACTTTGGTCCATCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTTGTGTGTCTAACCCCTGTTCTTTATACAGACGTTCAGATTAAGACCTTTGCTCATGGAAt        <  1:1518629/150‑1 (MQ=255)
                                                                                                                                                gTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTTGTGTGTCTAACCCCTGTTCTTTATACAGACGTTCAGATAAAGAACTTTGCTCATGGAATGTc     >  1:1178478/1‑150 (MQ=255)
                                                                                                                                                  tCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTTGTGTGTCTAACCCCTGTTCTTTATACAGACGTTCAGATAAAGACCTTTGCTCATGGAATGTCGc   >  1:1091747/1‑150 (MQ=255)
                                                                                                                                                   cacaACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTTGTGTGTCTAACCCCTGTTCTTTATACAGACGTTCAGATAAAGACATTTGCTCATGGAATGTCGCa  >  1:647499/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TTTCATCGATTTTAGGTGATCTTCTCTGTGATATTTCACGTCTTACCAGTAACTGATATGGCTGTCCGCCGCTCGCTTAAAGTGGACTTTTTAGTTTTTATCATGTGCGGTGAGAAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTTGTGTGTCTAACCCCTGTTCTTTATACAGACGTTCAGATAAAGACCTTTGCTCATGGAATGTCGCA  >  NC_000913/1199533‑1199829

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