Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913286,5690GA55.9% 113.9 / inf 356S392N (AGC→AAC) yagGCP4‑6 prophage; putative sugar transporter
Reads supporting (aligned to +/- strand):  ref base G (72/85);  new base A (90/109);  total (162/194)
Fisher's exact test for biased strand distribution p-value = 9.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.66e-01

Only 100 of 382 total aligned reads displayed.

TTCTGACGTGGTGGACTACGAGGAGAGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAGCAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCTCACCGATGCCCGCGTGGAGGCCATCAGCCGGCAGCTGATTAA  >  NC_000913/286420‑286719
                                                                                                                                                     |                                                                                                                                                      
    gACGTGGTGGACTACGAGGAGAGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCcaacagc                                                                                                                                                    >  3:1563731/1‑150 (MQ=255)
    gACGTGGTGGACTACGAGGAGAGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCcaaca                                                                                                                                                      <  2:929793/148‑1 (MQ=255)
     ccGTGGTGGACTACGAGGAGAGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCcagcagcag                                                                                                                                                  <  7:1504933/150‑1 (MQ=255)
     aCGTGGTGGACTACGAGGAGAGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCcagcagcag                                                                                                                                                  <  8:766796/151‑1 (MQ=255)
     aCGTGGTGGACTACGAGGAGAGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCcaacagc                                                                                                                                                    >  7:1073217/1‑149 (MQ=255)
       gtggtgGACTACGAGGAGAGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAGCAGCAGCg                                                                                                                                                <  2:1062035/151‑1 (MQ=255)
       gtggtgGACTACGAGGAGAGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCg                                                                                                                                                <  5:1410795/151‑1 (MQ=255)
          gtggACTACGTGGAGAGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAGCAGCAGCGTGc                                                                                                                                             <  5:890930/151‑1 (MQ=255)
          gtggACTACGAGGAGAGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAGCAGCAGCGTGCa                                                                                                                                            <  5:1089808/151‑1 (MQ=255)
            ggACTACGAGGAGAGCCGCAGCGGTCGCCGCCTCGAGGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCGTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTGAACTATTCCGCCAACAGCAGCGTGCAg                                                                                                                                           <  2:840182/151‑1 (MQ=255)
            ggACTACGAGGAGAGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAGCAGCAGCGTGCAg                                                                                                                                           <  6:1068213/151‑1 (MQ=255)
             gACTACGAGGAGAGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGc                                                                                                                                          <  5:1582700/151‑1 (MQ=255)
              aCTACGAGGAGAGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGcc                                                                                                                                         <  2:53917/151‑1 (MQ=255)
                tACGAGGAGAGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCgg                                                                                                                                       >  6:400409/1‑151 (MQ=255)
                 aCGAGGAGAGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAGCAGCAGCGTGCAGCCGGt                                                                                                                                      >  5:341682/1‑151 (MQ=255)
                 aCGAGGAGAGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCgg                                                                                                                                       <  1:334750/150‑1 (MQ=255)
                 aCGAGGAGAGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCg                                                                                                                                        <  6:211248/149‑1 (MQ=255)
                 aCGAGGAGAGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGt                                                                                                                                      <  3:193193/151‑1 (MQ=255)
                  cGAGGAGAGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGtt                                                                                                                                     <  3:1429808/151‑1 (MQ=255)
                    aggagAGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGa                                                                                                                                   <  3:766426/151‑1 (MQ=255)
                     ggagAGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAg                                                                                                                                  <  6:583012/151‑1 (MQ=255)
                         agCCGCAGCGGTCGCCGCCTCGCCGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGc                                                                                                                              >  3:1019807/1‑151 (MQ=255)
                             gCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAGCAGCAGCGTGCAGCCGGTTGAGGTGCTcac                                                                                                                          <  5:28320/151‑1 (MQ=255)
                               aGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCCCCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAGCAGCAGCGTGCAGCCGGTTGAGGTGCTcaccc                                                                                                                        >  8:613803/1‑150 (MQ=255)
                               aGCGGTCGCCGCCTCGACGGGCTGGTGGTCGCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAGCAGCAGCGGGCAGCCGGTTGAGGGGCTcacca                                                                                                                        <  3:1114374/151‑1 (MQ=255)
                                  ggTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCaccacc                                                                                                                      >  8:1013814/1‑150 (MQ=255)
                                   gTCGCCGCCTCGACGGGGTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGGGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAGCAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCAt                                                                                                                    <  4:620584/151‑1 (MQ=255)
                                        cgccTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAGCAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCaaa                                                                                                                >  8:1231053/1‑150 (MQ=255)
                                         gccTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAAt                                                                                                              <  6:1654723/151‑1 (MQ=255)
                                          ccTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAGCAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAAtt                                                                                                             <  1:1347213/151‑1 (MQ=255)
                                          ccTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAAtt                                                                                                             >  3:407394/1‑151 (MQ=255)
                                           cTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAAttc                                                                                                            >  8:401800/1‑151 (MQ=255)
                                               aCGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAAttctgtt                                                                                                        <  7:1578117/151‑1 (MQ=255)
                                                cGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAAttctgttc                                                                                                       >  4:284461/1‑151 (MQ=255)
                                                 gggCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGAGCCTGGCGTACGTCAACTATTCCGCCAGCAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAAttctgttct                                                                                                      <  1:288063/151‑1 (MQ=255)
                                                  ggCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAGCAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAAttctgttctg                                                                                                     <  1:1069941/151‑1 (MQ=255)
                                                   gCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAGCAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGc                                                                                                    <  6:467054/151‑1 (MQ=255)
                                                   gCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAGCAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGc                                                                                                    <  3:574782/151‑1 (MQ=255)
                                                      ggTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCg                                                                                                   <  2:1422386/149‑1 (MQ=255)
                                                        ggtgctcCACCTACCTGTGCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCgtggt                                                                                               <  4:1312174/147‑1 (MQ=255)
                                                        ggtTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAGCAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCgtggt                                                                                               >  1:1141178/2‑151 (MQ=255)
                                                           tctcCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGcc                                                                                            >  6:422402/1‑151 (MQ=255)
                                                              ccACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAGCAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCggt                                                                                         >  8:325959/1‑151 (MQ=255)
                                                              ccACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAGCAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCgg                                                                                          <  4:1019867/150‑1 (MQ=255)
                                                              ccACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCg                                                                                           <  8:482298/149‑1 (MQ=255)
                                                                  ctacctGTTCAGCCTGAAGATTGGCCTGGCGATCGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGgtggtg                                                                                     >  1:530222/1‑151 (MQ=255)
                                                                   tacctGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAGCAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGc                                                                                    >  8:1144732/1‑151 (MQ=255)
                                                                      ctGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAGCAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGctct                                                                                 >  5:1547559/1‑151 (MQ=255)
                                                                      ctGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGctct                                                                                 >  3:895759/1‑151 (MQ=255)
                                                                       tGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAGCAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGctct                                                                                 >  6:1339380/1‑150 (MQ=255)
                                                                        gTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAGCAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTAc                                                                               <  5:322274/151‑1 (MQ=255)
                                                                        gTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAGCAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTAc                                                                               >  5:758154/1‑151 (MQ=255)
                                                                        gTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGctct                                                                                 <  7:316466/149‑1 (MQ=255)
                                                                        gTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGctct                                                                                 <  8:85895/149‑1 (MQ=255)
                                                                          tCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTAcgc                                                                             >  1:693174/1‑151 (MQ=255)
                                                                           cAGCCTGACGATTGGCCTGGCGCTTGGCGGGGCGGTGGTTGTCTGGTTCCTGGCGTAAGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTAcgc                                                                             <  2:1590841/150‑1 (MQ=255)
                                                                           cAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTAcgc                                                                             <  2:1549818/150‑1 (MQ=255)
                                                                            aGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTCCTGCGTGGTGCCGGTGGTGCTCTACGCgg                                                                           <  6:823544/151‑1 (MQ=255)
                                                                             gCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCCGTTGAGGTGCTCACCCCCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCcgg                                                                          >  1:1712862/1‑151 (MQ=255)
                                                                              ccTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGc                                                                         >  7:431244/1‑151 (MQ=255)
                                                                               cTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAGCAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGc                                                                         >  5:127328/1‑150 (MQ=255)
                                                                               cTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGc                                                                         >  5:562965/1‑150 (MQ=255)
                                                                                tGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAGCAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCAt                                                                       >  1:92575/1‑151 (MQ=255)
                                                                                tGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCAt                                                                       >  2:1708774/1‑151 (MQ=255)
                                                                                 gAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATg                                                                      >  7:1012477/1‑151 (MQ=255)
                                                                                  aaGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGt                                                                     >  2:25815/1‑151 (MQ=255)
                                                                                     aTTGGCCTGGCGATTGGCGGGGCGGTGGTGGTCTGGATCCTGGCGTACGTCAACTATTCCGCCAGCAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTtca                                                                  <  6:294576/151‑1 (MQ=255)
                                                                                       tGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTtcat                                                                 <  8:1461916/150‑1 (MQ=255)
                                                                                        ggCCTGGCGATTGGCGGGGCGGTGGTAGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAGCAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGGGTGGTGCCGGTGGTGCTCTACGCGGGCATGTtcatc                                                                <  4:382432/150‑1 (MQ=255)
                                                                                          ccTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATg                                                             >  1:842082/1‑151 (MQ=255)
                                                                                               cGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAGCAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTc                                                        <  8:1582242/151‑1 (MQ=255)
                                                                                               cGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAGCAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTc                                                        >  7:746374/1‑151 (MQ=255)
                                                                                               cGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAGCAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTc                                                        >  2:228132/1‑151 (MQ=255)
                                                                                               cGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTc                                                        >  3:1302968/1‑151 (MQ=255)
                                                                                                gATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCg                                                       <  1:1468508/151‑1 (MQ=255)
                                                                                                 aTTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAGCAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGc                                                      >  8:1283185/1‑151 (MQ=255)
                                                                                                 aTTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGACGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGc                                                      >  4:475828/1‑151 (MQ=255)
                                                                                                   tGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGctc                                                    >  2:295036/1‑151 (MQ=255)
                                                                                                   tGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGctc                                                    >  2:1261756/1‑151 (MQ=255)
                                                                                                      cggggcggTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAGCAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTAc                                                 >  1:233407/1‑151 (MQ=255)
                                                                                                      cggggcggTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAGCAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTAc                                                 >  5:1711048/1‑151 (MQ=255)
                                                                                                        gggcggTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAGCAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACaa                                               <  2:74203/151‑1 (MQ=255)
                                                                                                         ggcggTGGTGGGGTGGATCCTGGCGTAGGTCAACTATTCCGCCAGCAGCAGCGTGCAGCCGGTTTAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCTGGCATGTTCATCATGCTGTCGCTCTACAAg                                              <  7:1373039/151‑1 (MQ=255)
                                                                                                               ggtggGCTGGATCCTGGCGTACGTCAACTATTCCGCCAGCAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCTCAc                                         <  3:462564/150‑1 (MQ=255)
                                                                                                                 tggGCTGGATCCTGGCGTACGTCAACTATTCCGCCAGCAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCTCACCg                                       >  7:1341671/1‑150 (MQ=255)
                                                                                                                   ggCTGGATCCTGGCGTACGTCAACTATTCCGACAGCAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCTCACCGATg                                    <  7:1705150/151‑1 (MQ=255)
                                                                                                                      gggATCCTGGCGTACGGCAACTATTCCGCCAGCAGCAGCGGGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCTCACCGATGccc                                 <  8:189509/150‑1 (MQ=255)
                                                                                                                          tCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCTCACCGATGCCCGCGt                             >  4:1498012/1‑151 (MQ=255)
                                                                                                                              ggCGTACGTCAACTATTCCGCCAGCAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCTCACCGATGCCCGCGTgg                           <  6:932855/149‑1 (MQ=255)
                                                                                                                               gCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCTCACCGATGCCCGCGTGGAgg                        >  7:1184622/1‑151 (MQ=255)
                                                                                                                               gCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTGCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCTCACCGATGCCCGCGTGGAgg                        <  6:11928/151‑1 (MQ=255)
                                                                                                                                    cGTCAACTATTCCGCCAGCAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCTCACCGATGCCCGCGTGGAGGCCATc                   >  1:457805/1‑151 (MQ=255)
                                                                                                                                      tCAACTATTCCGCCAGCAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCTCACCGATGCCCGCGTGGAGGCCATc                   >  3:1167055/1‑149 (MQ=255)
                                                                                                                                       cAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCTCACCGATGCCCGCGTGGAGGCCATCAGc                >  4:918916/1‑151 (MQ=255)
                                                                                                                                            aTTCCGCCAGCAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCTCACCGATGCCCGCGTGGAGGCCATCAGCCGGc            <  4:1609945/150‑1 (MQ=255)
                                                                                                                                              tCCGCCAGCAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCTCACCGATGCCCGCGTGGAGGCCATCAGCCGGCAg          <  4:1126987/150‑1 (MQ=255)
                                                                                                                                              tCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCTCACCGATGCCcgc                               >  3:265727/1‑129 (MQ=255)
                                                                                                                                                cGCCAACAGCAGCGTGCAGGCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCTCACCGATGCCCGCGTGGAGGCCATCAGCCGGCAGCTg       <  2:61673/151‑1 (MQ=255)
                                                                                                                                                 gCCAGCAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCTCACCGATGCCCGCGTGGAGGCCATCAGCCGGCAGCTg       >  4:87816/1‑150 (MQ=255)
                                                                                                                                                   cagcagcagCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCTCACCGATGCCCGCGTGGAGGCCATCAGCCGGCAGCTGAtt    >  4:245304/1‑151 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
TTCTGACGTGGTGGACTACGAGGAGAGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAGCAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCTCACCGATGCCCGCGTGGAGGCCATCAGCCGGCAGCTGATTAA  >  NC_000913/286420‑286719

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