Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,411,408 A→T 8.0% intergenic (‑56/+31) 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,4080AT 8.0% 311.3 / 11.4 120intergenic (‑56/+31)yfbS/yfbTputative transporter/sugar phosphatas
Reads supporting (aligned to +/- strand):  ref base A (55/52);  new base T (6/4);  total (63/57)
Fisher's exact test for biased strand distribution p-value = 7.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

Only 100 of 121 total aligned reads displayed.

ACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCGGGGTATCTGCCGGAGCGTTAACCGCAA  >  NC_000913/2411259‑2411554
                                                                                                                                                     |                                                                                                                                                  
     gCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCAATAAAAAAGCCCCGACgaga                                                                                                                                               >  2:809976/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)
                    cGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAAATGTGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATg                                                                                                                                >  1:1280904/1‑150 (MQ=255)
                        ctctGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGACCCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCa                                                                                                                            <  1:1089861/150‑1 (MQ=255)
                        ctctGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCa                                                                                                                            >  2:556755/1‑150 (MQ=255)
                        ctctGCCCGACCCAAACAAGACTCTGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCa                                                                                                                            <  2:248130/150‑1 (MQ=255)
                         tctGCCGGTCGCAACCAAGACGATGGCAACGGCCAGTAATGAAAGACCCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAg                                                                                                                           <  2:107204/150‑1 (MQ=255)
                         tctGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAg                                                                                                                           <  2:446423/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)
                                       aCAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGaat                                                                                                             <  1:1532609/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)
                                              gATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAATCCCCGACTCGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACa                                                                                                      >  1:1154315/1‑150 (MQ=255)
                                                  gCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGc                                                                                                  <  2:1100500/150‑1 (MQ=255)
                                                         ccAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTgg                                                                                           >  2:890672/1‑150 (MQ=255)
                                                              aaTGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGtt                                                                                      >  1:1637619/1‑150 (MQ=255)
                                                               aTGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGttt                                                                                     >  2:877739/1‑150 (MQ=255)
                                                                  aaaGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGt                                                                                  >  2:387599/1‑150 (MQ=255)
                                                                  aaaGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGt                                                                                  >  2:1342715/1‑150 (MQ=255)
                                                                      aaCCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTcaca                                                                              >  1:1095336/1‑150 (MQ=255)
                                                                       aCCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAg                                                                             >  1:1322876/1‑150 (MQ=255)
                                                                        cccAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGt                                                                            >  2:1103718/1‑150 (MQ=255)
                                                                          cAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTaa                                                                          >  1:901042/1‑150 (MQ=255)
                                                                          cAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTaa                                                                          >  2:852196/1‑150 (MQ=255)
                                                                            aaTCATTTCACCGTTCAAAAATTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAATAGCCCCGACGAGACGGGTCTAAATCATGATCAGGTGTTTCTCTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAAtt                                                                        <  2:849571/150‑1 (MQ=255)
                                                                               cAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGc                                                                     <  2:1583161/150‑1 (MQ=255)
                                                                               aaaTTCACCGTTCACAAATTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGc                                                                     <  2:354342/149‑1 (MQ=255)
                                                                                aaTTCACCGTTCACAACTTATCCTTGTTAATTGCGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCt                                                                    <  2:1076647/150‑1 (MQ=255)
                                                                                aaTTCACCGTTCACAACTTATCCTTGTAAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCt                                                                    <  1:225719/150‑1 (MQ=255)
                                                                                  ttCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTcc                                                                  >  1:1189860/1‑150 (MQ=255)
                                                                                  ttCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTcc                                                                  <  1:556755/150‑1 (MQ=255)
                                                                                  atCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTcc                                                                  <  1:637681/149‑1 (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)
                                                                                               aCTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTggagga                                                     <  2:677141/150‑1 (MQ=255)
                                                                                                attATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGAc                                                    <  2:1618343/149‑1 (MQ=255)
                                                                                                  tATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCa                                                  >  2:307037/1‑150 (MQ=255)
                                                                                                     ccTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGAc                                                    <  1:1258837/145‑1 (MQ=255)
                                                                                                      cTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATc                                              >  1:465734/1‑150 (MQ=255)
                                                                                                        tGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGCCTGTGGAGGACCAAATCGa                                            >  1:1313059/1‑150 (MQ=255)
                                                                                                         tttAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACCGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGAc                                           <  2:1280904/149‑1 (MQ=255)
                                                                                                           tAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCt                                         <  1:322050/150‑1 (MQ=255)
                                                                                                           tAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCt                                         <  2:1188328/150‑1 (MQ=255)
                                                                                                           tAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATAAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCt                                         >  2:337346/1‑150 (MQ=255)
                                                                                                            taTTGTGGTGGATGACTTTATTCTGCCATAAAAAAGCGCCGACGAGTCGGGGCTAAATCTTGATCAGGTGTTTCTATGAATAATCACATCGGCATTTGGCTGTTTGGTCACCGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTc                                        <  2:110626/149‑1 (MQ=255)
                                                                                                            aaTTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTc                                        >  1:838253/1‑150 (MQ=255)
                                                                                                            aaTTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTATGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTc                                        >  1:460412/1‑150 (MQ=255)
                                                                                                                gAGGGGGATGACTTGATTCTGCTATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTAATTc                                    >  1:216007/1‑150 (MQ=255)
                                                                                                                        tGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCgg                            <  2:828422/150‑1 (MQ=255)
                                                                                                                          aCTTGATTCTGCTATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCgggg                          <  2:216007/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)
                                                                                                                                  cTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCGGGGTATCTGcc                  >  1:990446/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)
                                                                                                                                          aaaaaGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCGGGGTATCTGCCGGAGCGtt          >  1:83278/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)
                                                                                                                                               gCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCGGGGTATCTGCCGGAGCGTTAACCg     >  1:717349/1‑150 (MQ=255)
                                                                                                                                                ccccGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCGGGGTATCTGCCGGAGCGTTAACCGc    >  2:764477/1‑150 (MQ=255)
                                                                                                                                                ccccGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCGGGGTATCTGCCGGAGCGTTAACCGc    >  1:1609931/1‑150 (MQ=255)
                                                                                                                                                ccccGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCGGGGTATCTGCCGGAGCGTTAACCGc    >  2:34515/1‑150 (MQ=255)
                                                                                                                                                  ccGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCGGGGTATCTGCCGGAGCGTTAACCGCaa  >  2:1269505/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
ACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCGGGGTATCTGCCGGAGCGTTAACCGCAA  >  NC_000913/2411259‑2411554

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