Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,411,411 A→C 8.2% intergenic (‑59/+28) yfbS ← / ← yfbT putative transporter/sugar phosphatas

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,411,4110AC 8.2% 335.8 / 16.8 119intergenic (‑59/+28)yfbS/yfbTputative transporter/sugar phosphatas
Reads supporting (aligned to +/- strand):  ref base A (54/54);  new base C (7/3);  total (62/57)
Fisher's exact test for biased strand distribution p-value = 3.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

Only 100 of 122 total aligned reads displayed.

GCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCGGGGTATCTGCCGGAGCGTTAACCGCAATGACAT  >  NC_000913/2411264‑2411560
                                                                                                                                                   |                                                                                                                                                     
gCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCAATAAAAAAGCCCCGACgaga                                                                                                                                                     >  2:809976/1‑150 (MQ=255)
  gACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACg                                                                                                                                                   >  1:439183/1‑150 (MQ=255)
     cGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACgggg                                                                                                                                                <  1:45411/150‑1 (MQ=255)
     cGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAAGTCACCGTTCACACCTTATCCTGGTTAATTGAGGGGGATGACTTGCTTCTGCCATAAAAAAGCCCCGCCGAGACgggg                                                                                                                                                >  1:1206906/1‑150 (MQ=255)
       cATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCt                                                                                                                                              >  1:1589515/1‑150 (MQ=255)
        aTCCATACGCACTCTGCCCGTCTCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTGGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTa                                                                                                                                             <  2:1331035/150‑1 (MQ=255)
         tCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGTCTCGTCGGGGCttt                                                                                                                                            <  1:1515343/150‑3 (MQ=255)
           cATACGCACTCTGCCCGTCGCAATCAAGACGATGGCAACGGCCAGTATTGCAAGAACCCATATCAATTCACCGTTCACACCTTATCCTTGTGAATTGAGGGGGCTGCCTTGCTTCTGCCATAAAAAAGCCCCGACGTGACGGGGCTTAAt                                                                                                                                          >  2:291025/1‑150 (MQ=255)
           cATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGTCTCGTCGGGGCttttt                                                                                                                                          <  1:542643/150‑5 (MQ=255)
           cATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTACACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGTCTCGTCGGGGCttttt                                                                                                                                          <  1:1006259/150‑5 (MQ=255)
            aTACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATc                                                                                                                                         >  2:1245253/1‑150 (MQ=255)
               cGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATg                                                                                                                                      <  1:1574860/150‑1 (MQ=255)
               cGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATg                                                                                                                                      >  1:183195/1‑150 (MQ=255)
               cGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAAATGTGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATg                                                                                                                                      >  1:1280904/1‑150 (MQ=255)
                   ctctGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGACCCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCa                                                                                                                                  <  1:1089861/150‑1 (MQ=255)
                   ctctGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCa                                                                                                                                  >  1:1583161/1‑150 (MQ=255)
                   ctctGCCCGACCCAAACAAGACTCTGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCa                                                                                                                                  <  2:248130/150‑1 (MQ=255)
                    tctGCCGGTCGCAACCAAGACGATGGCAACGGCCAGTAATGAAAGACCCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAg                                                                                                                                 <  2:107204/150‑1 (MQ=255)
                     ctGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAgg                                                                                                                                <  1:1191832/150‑1 (MQ=255)
                       gCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGgtg                                                                                                                              <  2:1523113/150‑1 (MQ=255)
                        cccGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGgtgt                                                                                                                             <  2:480200/150‑1 (MQ=255)
                        cccGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGgtgt                                                                                                                             >  1:1188328/1‑150 (MQ=255)
                         ccGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGtt                                                                                                                            <  1:1001113/150‑1 (MQ=255)
                             cGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATGCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCAc                                                                                                                        <  2:1293794/150‑1 (MQ=255)
                              gCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGTGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACt                                                                                                                       <  1:99876/150‑1 (MQ=255)
                               cAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTg                                                                                                                      <  1:545342/150‑1 (MQ=255)
                               cAAACAACACGATGGCAACGGCCAGTAATGAAAGAACCAAAATCAATTCACCGTACACAACTTATCCTTGTAAAAAGAGGGGGATGACTAGATTCTGCCATAAAAAAGCCACGACCAGACGGCGCTAAATCATGATAAGGTGTTACACTg                                                                                                                      >  1:1661896/1‑150 (MQ=255)
                                aaaCAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGa                                                                                                                     >  1:993827/1‑150 (MQ=255)
                                 aaCAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGaa                                                                                                                    <  2:1269595/150‑1 (MQ=255)
                                  aCAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGaat                                                                                                                   <  1:232984/150‑1 (MQ=255)
                                   cAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGaata                                                                                                                  <  1:1000066/150‑1 (MQ=255)
                                      gACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAataata                                                                                                               >  2:1569228/1‑150 (MQ=255)
                                       aCGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCTTTCACAACTTATCCTTGTTAATTGAGGGGTATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAAtaataa                                                                                                              <  1:1459143/150‑1 (MQ=255)
                                        cGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAAc                                                                                                             <  2:535818/150‑1 (MQ=255)
                                        cGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAAc                                                                                                             >  1:555786/1‑150 (MQ=255)
                                         gATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAATCCCCGACTCGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACa                                                                                                            >  1:1154315/1‑150 (MQ=255)
                                             gCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGc                                                                                                        <  2:1100500/150‑1 (MQ=255)
                                                    ccAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTgg                                                                                                 >  2:890672/1‑150 (MQ=255)
                                                          aTGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGttt                                                                                           >  2:877739/1‑150 (MQ=255)
                                                             aaaGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGt                                                                                        <  2:439183/150‑1 (MQ=255)
                                                             aaaGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGt                                                                                        >  2:387599/1‑150 (MQ=255)
                                                              aaGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTc                                                                                       <  1:1408986/150‑1 (MQ=255)
                                                               aGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTca                                                                                      >  2:1097499/1‑150 (MQ=255)
                                                                 aaCCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTcaca                                                                                    >  1:1095336/1‑150 (MQ=255)
                                                                   cccAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGt                                                                                  >  2:1103718/1‑150 (MQ=255)
                                                                    ccAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTa                                                                                 <  1:809976/150‑1 (MQ=255)
                                                                     cAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTaa                                                                                >  2:852196/1‑150 (MQ=255)
                                                                     cAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTaa                                                                                >  1:901042/1‑150 (MQ=255)
                                                                       aaTCATTTCACCGTTCAAAAATTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAATAGCCCCGACGAGACGGGTCTAAATCATGATCAGGTGTTTCTCTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAAtt                                                                              <  2:849571/150‑1 (MQ=255)
                                                                          aaaTTCACCGTTCACAAATTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGc                                                                           <  2:354342/149‑1 (MQ=255)
                                                                           aaTTCACCGTTCACAACTTATCCTTGTAAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCt                                                                          <  1:225719/150‑1 (MQ=255)
                                                                            aTTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGATACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTc                                                                         >  2:322050/1‑150 (MQ=255)
                                                                             ttCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTcc                                                                        <  1:556755/150‑1 (MQ=255)
                                                                             ttCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTcc                                                                        >  1:1189860/1‑150 (MQ=255)
                                                                                aCCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGa                                                                     >  2:71856/1‑150 (MQ=255)
                                                                                 ccGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGAc                                                                    <  2:1083185/150‑1 (MQ=255)
                                                                                  cGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACt                                                                   >  1:601524/1‑150 (MQ=255)
                                                                                   gTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACtg                                                                  >  1:429763/1‑150 (MQ=255)
                                                                                      cacaACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTgg                                                               >  1:1245757/1‑150 (MQ=255)
                                                                                       acaACTAATCCTTGTTAATTGATGGGGATGAGTTGATGCTTCCATAAAATAGCAACTTCGATACTGGGATATATAATGATTAGGTGTTGCACTGAAAAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTgga                                                              <  2:547178/150‑1 (MQ=255)
                                                                                         aaCTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTggagg                                                            <  1:1501347/150‑1 (MQ=255)
                                                                                         aaCTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTggagg                                                            <  2:555786/150‑1 (MQ=255)
                                                                                          aCTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTggagga                                                           <  2:677141/150‑1 (MQ=255)
                                                                                           attATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGAc                                                          <  2:1618343/149‑1 (MQ=255)
                                                                                             tATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCa                                                        >  2:307037/1‑150 (MQ=255)
                                                                                                ccTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGAc                                                          >  2:1258837/1‑145 (MQ=255)
                                                                                                ccTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGAc                                                          <  1:1258837/145‑1 (MQ=255)
                                                                                                   tGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGCCTGTGGAGGACCAAATCGa                                                  >  1:1313059/1‑150 (MQ=255)
                                                                                                     ttAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGAcc                                                <  2:194220/150‑1 (MQ=255)
                                                                                                      tAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCt                                               <  2:1188328/150‑1 (MQ=255)
                                                                                                       aaTTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTc                                              >  1:838253/1‑150 (MQ=255)
                                                                                                       aaTTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTATGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTc                                              >  1:460412/1‑150 (MQ=255)
                                                                                                           gAGGGGGATGACTTGATTCTGCTATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTAATTc                                          >  1:216007/1‑150 (MQ=255)
                                                                                                              ggggATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAgg                                       <  1:890672/150‑1 (MQ=255)
                                                                                                              ggggATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAgg                                       >  2:1457747/1‑150 (MQ=255)
                                                                                                                   tGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCgg                                  <  2:828422/150‑1 (MQ=255)
                                                                                                                     aCTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCgggg                                <  1:852196/150‑1 (MQ=255)
                                                                                                                     aCTTGATAGTGCCATAAAATAGCACAGCCGAGACGGGGATAAATCATGATGAGGTGTTTCACTTAATAATAACATCGCCATTTGGCTGTTTGGTAAAAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCgggg                                <  2:314405/150‑1 (MQ=255)
                                                                                                                       ttGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCGGGGTa                              <  2:1637619/150‑1 (MQ=255)
                                                                                                                             cTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCGGGGTATCTGcc                        >  1:707312/1‑150 (MQ=255)
                                                                                                                              tGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCGGGGTATCTGCCg                       <  2:138709/150‑1 (MQ=255)
                                                                                                                                 cATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGCTCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGCCCAAATCGACCTCATTCAGGCGCGGGGTATCTGCCGGAg                    >  1:82624/1‑150 (MQ=255)
                                                                                                                                 cATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCGGGGTATCTGCCGGAg                    >  2:592898/1‑150 (MQ=255)
                                                                                                                                   tgatttaGCCCCGTCTCGTCGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCGGGGTATCTGCCGGAGCg                  >  1:909624/7‑150 (MQ=255)
                                                                                                                                   tAAAAAAGCCCCGCCGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCGGGGTATCTGCCGGAGCg                  >  2:1072523/1‑150 (MQ=255)
                                                                                                                                      tttaGCCCCGTCTCGTCGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATAGACCTCATTCAGGCGCGGGGTATCTGCCGGAGCGTTa               >  2:686417/4‑150 (MQ=255)
                                                                                                                                      aaaaGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCGGGGTATCTGCCGGAGCGTTa               <  2:1322876/150‑1 (MQ=255)
                                                                                                                                      aaaaGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCGGGGTATCTGCCGGAGCGTTa               >  2:659000/1‑150 (MQ=255)
                                                                                                                                        aaGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCGGGGTATCTGCCGGAGCGTTAAc             <  2:1095336/150‑1 (MQ=255)
                                                                                                                                         aGCCCCGTCTCGTCGGGGCTAAATCATGATCAGGTGTTTCACTGTATAATAACATCGCAATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGCGGCCCCAATCGACCGCAGTCAGGCGCGGGGTATCTGCCGGAGCGTTAAcc            >  2:836703/1‑150 (MQ=255)
                                                                                                                                         aGCCCCGTCTCGTCGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCGGGGTATCTGCCGGAGCGTTAAcc            >  2:1346904/1‑150 (MQ=255)
                                                                                                                                         aGCCCCGTCTCGTCGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCGGGGTATCTGCCGGAGCGTTAAcc            >  2:473646/1‑150 (MQ=255)
                                                                                                                                         aGCCCCGTCTCGTCGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCGGGGTATCTGCCGGAGCGTTAAcc            >  2:1308382/1‑150 (MQ=255)
                                                                                                                                          gCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCGGGGTATCTGCCGGAGCGTTAACCg           >  1:717349/1‑150 (MQ=255)
                                                                                                                                           ccccGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCGGGGTATCTGCCGGAGCGTTAACCGc          >  2:34515/1‑150 (MQ=255)
                                                                                                                                           ccccGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCGGGGTATCTGCCGGAGCGTTAACCGc          >  2:764477/1‑150 (MQ=255)
                                                                                                                                           ccccGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCGGGGTATCTGCCGGAGCGTTAACCGc          >  1:1609931/1‑150 (MQ=255)
                                                                                                                                             ccGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCGGGGTATCTGCCGGAGCGTTAACCGCaa        >  2:1269505/1‑150 (MQ=255)
                                                                                                                                                   agaCGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCGGGGTATCTGCCGGAGCGTTAACCGCAatgacat  >  1:1439494/1‑150 (MQ=255)
                                                                                                                                                   agaCGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCGGGGTATCTGCCGGAGCGTTAACCGCAatgacat  <  2:1589515/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCGGGGTATCTGCCGGAGCGTTAACCGCAATGACAT  >  NC_000913/2411264‑2411560

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: