Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,199,678 T→A 5.3% intergenic (+90/+1) lit → / ← intE T4 phage exclusion protein; cell death peptidase, e14 prophage/e14 prophage; putative integrase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,199,6780TA 5.3% 278.0 / 3.9 95intergenic (+90/+1)lit/intET4 phage exclusion protein; cell death peptidase, e14 prophage/e14 prophage; putative integrase
Reads supporting (aligned to +/- strand):  ref base T (45/45);  new base A (3/2);  total (48/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01

TCCTTTTCATCGATTTTAGGTGATCTTCTCTGTGATATTTCACGTCTTACCAGTAACTGATATGGCTGTCCGCCGCTCGCTTAAAGTGGACTTTTTAGTTTTTATCATGTGCGGTGAGAAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTTGTGTGTCTAACCCCTGTTCTTTATACAGACGTTCAGATAAAGACCTTTGCTCATGGAATGTC  >  NC_000913/1199529‑1199826
                                                                                                                                                     |                                                                                                                                                    
tCCTTTTCATCGATTTTAGGTGATCTTCTCTGTGATATTTCATGTCTTACCAGTAACTGATATGGCTGTCCGCCGCTCGCTTAAAGTGGACTTTTTAGTTTTTATCATGTGCGGTGAGAAATTCAATGTGGCGTTGAGATGCTTAAAGGt                                                                                                                                                      <  1:1257440/150‑1 (MQ=255)
  cTTTTCATCGATTTTAGGTGATCTTCTCTGTGATATTTCACGTCTTACCAGTAACTGATATGGCTGTCCGCCGCTCGCTTAAAGTGGACTTTTTAGTTTTTATCATGTGCGGTGAGAAATTCAATGTGGCGTTGAGATGCTTAAAGGTTc                                                                                                                                                    <  1:1029061/150‑1 (MQ=255)
  cTTTTCATCGATTTTAGGTGATCTTCTCTGTGATATTTCACGTCTTACCAGTAACTGATATGGCTGTCCGCCGCTCGCTTAAAGTGGACTTTTTAGTATTTATCATGTGAGGTGAGAAATTCAATGTGGCGTTGAGATGCTTAAAGGTTc                                                                                                                                                    >  2:1515885/1‑150 (MQ=255)
    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                                                                                           >  1:1709978/1‑150 (MQ=255)
                                                           atatGGCTGTCCGCCGCTCGCTTAAAGTGGACTTTTTAGTTTTTATCATGTGCGGTGAGAAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTttgtt                                                                                           <  2:976891/150‑1 (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:653394/1‑150 (MQ=255)
                                                                     ccgccgCTCGCTTAAAGTGGACTTTTTAGTTTTTATCATGTGCGGTGAGAAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGtttt                                                                                 <  1:419827/150‑1 (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                                        <  2:487834/150‑1 (MQ=255)
                                                                                                              gCGGTGAGAAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTTGTGTGTCTAACCCCTGTTCTTtat                                        >  1:182544/1‑150 (MQ=255)
                                                                                                                  tGAGAAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTTGTGTGTCTAACCCCTGTTCTTTATACAg                                    <  1:808804/150‑1 (MQ=255)
                                                                                                                   gagaAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTTGTGTGTCTAACCCCTGTTCTTtata                                       <  2:809735/146‑1 (MQ=255)
                                                                                                                   gagaAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTTGTGTGTCTAACCCCTGTTCTTtata                                       >  1:809735/1‑146 (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:1135547/150‑1 (MQ=255)
                                                                                                                         ttCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTTGTGTGTCTAACCCCTGTTCTTTATACAGACGTTCa                             <  2:653394/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)
                                                                                                                                                     |                                                                                                                                                    
TCCTTTTCATCGATTTTAGGTGATCTTCTCTGTGATATTTCACGTCTTACCAGTAACTGATATGGCTGTCCGCCGCTCGCTTAAAGTGGACTTTTTAGTTTTTATCATGTGCGGTGAGAAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTTGTGTGTCTAACCCCTGTTCTTTATACAGACGTTCAGATAAAGACCTTTGCTCATGGAATGTC  >  NC_000913/1199529‑1199826

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: