Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,199,673 A→T 5.7% intergenic (+85/+6) 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,6730AT 5.7% 257.1 / 2.3 94intergenic (+85/+6)lit/intET4 phage exclusion protein; cell death peptidase, e14 prophage/e14 prophage; putative integrase
Reads supporting (aligned to +/- strand):  ref base A (44/44);  new base T (3/3);  total (47/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 = 1.00e+00

GGGAGTCCTTTTCATCGATTTTAGGTGATCTTCTCTGTGATATTTCACGTCTTACCAGTAACTGATATGGCTGTCCGCCGCTCGCTTAAAGTGGACTTTTTAGTTTTTATCATGTGCGGTGAGAAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTTGTGTGTCTAACCCCTGTTCTTTATACAGACGTTCAGATAAAGACCTTTGCTCATGGA  >  NC_000913/1199524‑1199821
                                                                                                                                                     |                                                                                                                                                    
gggAGTCCTGTTCATCGATTTTAGGTGATCTTCTCCGTGATATTTCACGTCTTACCAGTAACTGATATGGCTGTCCGCCGCTCGCTTAAAGTGGACTTTATAGTTTTTATCATGTGCGCTGAGAAATTCAATGTGGCGTTGAGATGCTTa                                                                                                                                                      >  1:310047/1‑150 (MQ=255)
  gAGTCCTTTTCATCGATTTTAGGTGATCTTCTCTGTGATATTTCACGTCTTACCAGTAACTGATATGGCTGTCCGCCGCTCGCTTAAAGTGGACTTTTTAGTTTTTATCATGTGCGGTGAGAAATTCAATGTGGCGTTGAGATGCTTaaa                                                                                                                                                    <  1:616845/150‑1 (MQ=255)
  gAGTCCTTTTCATCGATTTTAGGTGATCTTCTCTGTGATATTTCACGTCTTACCAGTAACTGATATGGCTGTCCGCCGCTCGCTTAAAGTGGACTTTTTAGTTTTTATCATGTGCGGTGAGAAATTCAATGTGGCGTTGAGATGCTTaaa                                                                                                                                                    <  1:586004/150‑1 (MQ=255)
     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                                                                                      <  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: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)
                                                                             ccgcTCGCTTAAAGTGGACTTTTTAGTTTTTATCATGTGCGGTGAGAAATTCAATGTGGCGTTGAGATGC‑TAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTATGa                                                                        >  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                                  >  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)
                                                                                                                            aaTTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCT‑CATCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTTGTGTGTCTAACCCCTGTTCTTTATACAGACGTTc                         <  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)
                                                                                                                                                     |                                                                                                                                                    
GGGAGTCCTTTTCATCGATTTTAGGTGATCTTCTCTGTGATATTTCACGTCTTACCAGTAACTGATATGGCTGTCCGCCGCTCGCTTAAAGTGGACTTTTTAGTTTTTATCATGTGCGGTGAGAAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTTGTGTGTCTAACCCCTGTTCTTTATACAGACGTTCAGATAAAGACCTTTGCTCATGGA  >  NC_000913/1199524‑1199821

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: ‑