Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 579,651 C→T 100% intergenic (‑57/‑229) ybcV ← / → ybcW DLP12 prophage; uncharacterized protein/DLP12 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913579,6510CT99.0% 13561.8 / 26.1 4689intergenic (‑57/‑229)ybcV/ybcWDLP12 prophage; uncharacterized protein/DLP12 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (21/30);  new base T (2110/2500);  total (2146/2543)
Fisher's exact test for biased strand distribution p-value = 5.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 4818 total aligned reads displayed.

GTGACGTTTTAGTTCAGAATAAAACAATTCACAGTTTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGCTTCAATCTGGTCTGACCTCTTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCAATTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACCGACAGATGTATGTAAGG  >  NC_000913/579505‑579796
                                                                                                                                                  |                                                                                                                                                 
 tgtcGTTTTATTTCAGAATAAAAGAATTCTCATTCTAAATCTTTTCGCACTTGATCGAATATTTGTTTAAAAATGGCAACCTGAGCCATTGGTAAAACGTTCCATGTGATACGAGGGCGCGTTGTTTGCATTATCGTTTTTATCGTTTc                                                                                                                                                <  2:56317/146‑1 (MQ=255)
 tgACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACTAGGGCGCGTAGTTTGCAttatcgtttttatcgttt                                                                                                                                                 >  1:246661/1‑148 (MQ=255)
 tgACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCAttatcgtttttatcgttt                                                                                                                                                 <  2:391608/148‑1 (MQ=255)
   aCGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGATGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCaa                                                                                                                                              >  1:218887/1‑149 (MQ=255)
    cGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAAt                                                                                                                                             <  1:276552/149‑1 (MQ=255)
        ttAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTg                                                                                                                                          <  8:247844/148‑1 (MQ=255)
           gTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCAATTGGTAAAACCTTCCATGTGATACGAGGGCGCGAAGTTTGCAATATCGTTGTTATCGTTTCAATCTGGACt                                                                                                                                      >  7:312650/1‑149 (MQ=255)
            ttCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCt                                                                                                                                      >  2:361374/1‑148 (MQ=255)
            ttCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCt                                                                                                                                      >  7:113370/1‑148 (MQ=255)
            ttCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGAAACCTGAGCCATTGGTAAAACCTTCCATGTGATACTAGGTCGCGTAGTTTGCATTATAGTTTTTATCGTTTCAATCTGGTCt                                                                                                                                      >  4:261121/1‑148 (MQ=255)
              ctgAATTATACAATTCACAGTCTTATTCTTTTGGCAGTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCTTGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGa                                                                                                                                    <  8:326882/146‑1 (MQ=255)
                gAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACct                                                                                                                                 >  4:350962/1‑149 (MQ=255)
                gAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACct                                                                                                                                 >  7:164838/1‑149 (MQ=255)
                  aTAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACctcc                                                                                                                               >  7:81544/1‑148 (MQ=255)
                   tAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACctcct                                                                                                                              >  3:122414/1‑147 (MQ=255)
                     aaaCAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCATTCTGGTCTGACCTCCTtg                                                                                                                            >  2:255147/1‑149 (MQ=255)
                     aaaCAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCctt                                                                                                                             <  7:288583/148‑1 (MQ=255)
                     aaaCAATTCACAGTCTAAATCTGTTCGCACTTGATCGAAGATTTCTTTAAAAATGGCAACCTGAGCCATTGGTATCACCTTCGATGTGAGGCGAGGGCGCGTAGTTGGCGTTATCGTTTTTATCGGTTCAATCTGGTCTGACCTCCTtg                                                                                                                            >  6:89480/1‑149 (MQ=255)
                       aCAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTtgtg                                                                                                                          >  2:17263/1‑149 (MQ=255)
                         aaTTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTTGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGtt                                                                                                                        <  2:229375/149‑1 (MQ=255)
                          aTTCACAGTTTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATAGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATGGCTTCAATCTGGTCTGACCTCTTTGTGttt                                                                                                                       >  2:278345/1‑149 (MQ=255)
                             cacaGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTttg                                                                                                                     <  3:8978/148‑1 (MQ=255)
                              acaGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAAt                                                                                                                                             <  8:15920/123‑1 (MQ=255)
                              acaGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTttgtt                                                                                                                   <  7:224737/149‑1 (MQ=255)
                              acaGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTttgtt                                                                                                                   >  1:134205/1‑149 (MQ=255)
                              acaGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTttgt                                                                                                                    <  8:271697/148‑1 (MQ=255)
                                    tAAATCTTTTCGCTCTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAAACTTCCTTGTGATACGAGGGCGCGTAGTTTGCAGTATAGGTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTgatgat                                                                                                             >  6:155369/1‑149 (MQ=255)
                                    tAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTgatgat                                                                                                             >  2:83665/1‑149 (MQ=255)
                                      aaTCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGAttt                                                                                                           <  4:183337/149‑1 (MQ=255)
                                      aaTCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGAtt                                                                                                            >  5:14885/1‑148 (MQ=255)
                                         cTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTtga                                                                                                                 <  1:298915/140‑1 (MQ=255)
                                          ttttCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCACGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGt                                                                                                       >  7:190963/1‑149 (MQ=255)
                                            ttCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGt                                                                                                       <  4:73115/147‑1 (MQ=255)
                                               gCACTTGATCGAATATTTCTTTAATAATGGCAACCTGAGCCATTGGTAAAACCTTGCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCaaa                                                                                                   <  2:113710/148‑1 (MQ=255)
                                                cACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAata                                                                                                 <  6:306779/149‑1 (MQ=255)
                                                     gATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTa                                                                                              <  3:59174/147‑1 (MQ=255)
                                                     gATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAg                                                                                             >  2:303818/1‑148 (MQ=255)
                                                     gATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAg                                                                                             >  6:27623/1‑148 (MQ=255)
                                                     gATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGCTTCAATCTGGTCTGACCTCTTTGTGTTTTGTTGATGATTTATGTCAAATATTa                                                                                              <  7:330258/147‑1 (MQ=255)
                                                      aTCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAgg                                                                                            <  1:100192/148‑1 (MQ=255)
                                                      aTCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAgg                                                                                            >  1:355519/1‑148 (MQ=255)
                                                      aTCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAgg                                                                                            <  8:106532/148‑1 (MQ=255)
                                                       tCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTTGTAAAACCTTCCATGTGATACGAGGGCGAGTAGTTTGCATTTTCGTTTTTATCGTTGCAATCTGGTCTGCCCTCCTTGAGATATGTTGATGATTTATGTCAAATATTTGGaa                                                                                          >  4:294973/1‑149 (MQ=255)
                                                         gAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATg                                                                                        >  2:202166/1‑149 (MQ=255)
                                                           atatTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGtt                                                                                      <  3:349611/149‑1 (MQ=255)
                                                              tttCTTTATATATGGGAACCTGAGCCATTGGTAAAACCTTCCTTGTGTTACGAGGTCGCGTATTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGGTGATGATTTATGTCAAATTTTAGGAATGtttt                                                                                    <  5:98104/148‑1 (MQ=255)
                                                              tttCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGtttt                                                                                    <  3:175348/148‑1 (MQ=255)
                                                               ttCTTTAAAAATGGCAAGCTGATCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGGTTTTATCTTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTc                                                                                   <  3:290865/148‑1 (MQ=255)
                                                               ttCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCa                                                                                  >  3:199394/1‑149 (MQ=255)
                                                                tCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCa                                                                                  <  5:309046/148‑1 (MQ=255)
                                                                 cTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGa                                                                                           >  5:274215/1‑138 (MQ=255)
                                                                  tttAAAAATGGCAACCTGAGCCATTGGTAACACCTTCCATGTGATACGAGGGCGCGCAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCGCCTTGTGTTTTGTTGAAGAATTATGTCAAATATTAGGAATGTTTTCact                                                                                >  6:239519/1‑146 (MQ=255)
                                                                  tttAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTa                                                                                              >  4:234375/1‑134 (MQ=255)
                                                                    tAAAAATGGCAACCTGAGCCATTGGTAAACCCTTCCATGTGATACAAGGGCGCGAAGATTGCATTATAGTTTTAATCGTTTCAATCTGGTATGCCCTCCTGGTGTTTTGTTGATGCTTTATGTCAAATATTAGGAAGGATTTCACTTaa                                                                             >  6:212483/1‑149 (MQ=255)
                                                                      aaaaTGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAAt                                                                            <  8:221128/148‑1 (MQ=255)
                                                                            gCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTAtt                                                                      <  5:49111/148‑1 (MQ=255)
                                                                            gCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATCAGGAATGTTTTCACTTAATAGTAtt                                                                      <  1:188688/148‑1 (MQ=255)
                                                                             cAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTgg                                                                    <  3:14680/149‑1 (MQ=255)
                                                                               aCCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGCGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGtt                                                                  <  8:188412/149‑1 (MQ=255)
                                                                                 cTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGTTGc                                                                <  5:237865/149‑1 (MQ=255)
                                                                                     gCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTaa                                                            <  7:58686/149‑1 (MQ=255)
                                                                                        aTTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTAACa                                                          >  1:161189/1‑148 (MQ=255)
                                                                                        aTTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTAACa                                                          >  5:176289/1‑148 (MQ=255)
                                                                                         ttGGTAAAACCTTCCATGTTTTACGAGTGCGCGTAGTTTTCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTAACaaa                                                        <  1:53187/149‑1 (MQ=255)
                                                                                         ttGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTAACaaa                                                        <  3:318775/149‑1 (MQ=255)
                                                                                             tAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTAACAAAGTGc                                                    >  8:132742/1‑149 (MQ=255)
                                                                                              aaaaCCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCaaa                                                                                                   <  4:44353/101‑1 (MQ=255)
                                                                                                       cATTTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTAACAAAGTGCGGTCctgct                                           <  7:142612/148‑1 (MQ=255)
                                                                                                       cATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTAACAAAGTGCGGTCctgct                                           <  8:70518/148‑1 (MQ=255)
                                                                                                        aTGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTgg                                         <  2:335860/149‑1 (MQ=255)
                                                                                                         tgtgATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTgg                                         >  8:302569/1‑148 (MQ=255)
                                                                                                            gATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTa                                                             <  8:376646/125‑1 (MQ=255)
                                                                                                            gATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTAACAAAGTGCgg                                                  <  6:333331/136‑1 (MQ=255)
                                                                                                              tACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCt                                   >  8:44287/1‑149 (MQ=255)
                                                                                                                 gAGGGCGCGTAGTTTGCATTTTCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATATTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGa                                <  4:322590/149‑1 (MQ=255)
                                                                                                                 gAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTg                                  >  3:231062/1‑147 (MQ=255)
                                                                                                                  aGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGGCAAATATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAg                               <  5:118048/149‑1 (MQ=255)
                                                                                                                   gggCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGt                                                              <  4:158068/117‑1 (MQ=255)
                                                                                                                   gggCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAgg                              <  5:72110/149‑1 (MQ=255)
                                                                                                                     gcgcGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAggg                             >  4:127609/1‑148 (MQ=255)
                                                                                                                     gcgcGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGAAGGGa                            <  6:180956/149‑1 (MQ=255)
                                                                                                                        cgTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAAt                         <  1:75406/149‑1 (MQ=255)
                                                                                                                         gTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATa                        >  7:257504/1‑149 (MQ=255)
                                                                                                                           aGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGtttt                                                                                    >  4:97283/1‑87 (MQ=255)
                                                                                                                             tttGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAAc                    <  6:64481/149‑1 (MQ=255)
                                                                                                                             ttggCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGGGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAAc                    <  6:36240/146‑1 (MQ=255)
                                                                                                                               tGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCGTTGTGTGTTGTTGATGATTAATGTCAAATATTAGGAATGGTTTCACTTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAAcc                   <  1:329263/148‑1 (MQ=255)
                                                                                                                               tGCATTATCGTTTTTATCGTTTCAATCTGGGCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAAcc                   <  6:123387/148‑1 (MQ=255)
                                                                                                                                 cATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACCGa                 <  2:147415/148‑1 (MQ=255)
                                                                                                                                  aTTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACCGa                 >  5:368712/1‑147 (MQ=255)
                                                                                                                                   ttatcgtttttatcgtttCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACACCCGa                 >  6:40008/1‑146 (MQ=255)
                                                                                                                                    tATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGaaa                          >  4:208548/1‑136 (MQ=255)
                                                                                                                                      tCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGTTGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACCGACAGat            <  5:208893/148‑1 (MQ=255)
                                                                                                                                      tCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACCGACAGat            <  7:356578/148‑1 (MQ=255)
                                                                                                                                      tCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACCGACAGat            <  1:375770/148‑1 (MQ=255)
                                                                                                                                        gTTTTTATAGTTTGAATGTGGTCTGACCTCCTTGTGTTTTTTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATa                                                                           <  5:344467/83‑1 (MQ=255)
                                                                                                                                            ttATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGTGGGAAATACAACCGACAGatgtatgt      >  3:38935/1‑148 (MQ=255)
                                                                                                                                            ttATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACCGACAGatgtatgt      <  3:259280/148‑1 (MQ=255)
                                                                                                                                              aTCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAAcc                   <  7:385255/133‑1 (MQ=255)
                                                                                                                                               tCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACCGACAGATGTATGTAAg   <  8:354276/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GTGACGTTTTAGTTCAGAATAAAACAATTCACAGTTTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGCTTCAATCTGGTCTGACCTCTTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCAATTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACCGACAGATGTATGTAAGG  >  NC_000913/579505‑579796

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: