Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,411,408 A→T 5.8% 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 5.8% 364.2 / 9.7 103intergenic (‑56/+31)yfbS/yfbTputative transporter/sugar phosphatas
Reads supporting (aligned to +/- strand):  ref base A (50/47);  new base T (2/4);  total (52/51)
Fisher's exact test for biased strand distribution p-value = 4.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01

Only 100 of 108 total aligned reads displayed.

CAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCGGGGTATCTGCCGGAGCGTTAACCGCAA  >  NC_000913/2411260‑2411554
                                                                                                                                                    |                                                                                                                                                  
caaaGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCcga                                                                                                                                                    <  1:70369/149‑1 (MQ=255)
 aaaGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCcgac                                                                                                                                                   <  3:634505/149‑1 (MQ=255)
 aaaGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCcgac                                                                                                                                                   >  3:111943/1‑149 (MQ=255)
       aCCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACg                                                                                                                                             >  4:125324/1‑149 (MQ=255)
         cGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACggg                                                                                                                                           <  2:1090529/149‑1 (MQ=255)
             tCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTa                                                                                                                                       >  3:617675/1‑149 (MQ=255)
              ccATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGTCTCGTCGGGGCttt                                                                                                                                      <  2:599091/149‑3 (MQ=255)
               cATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGTCTCGTCGGGGCtttt                                                                                                                                     <  4:579418/149‑4 (MQ=255)
                aTACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAAt                                                                                                                                    >  1:1073450/1‑149 (MQ=255)
                 tACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATc                                                                                                                                   >  3:26269/1‑149 (MQ=255)
                      aCTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGAt                                                                                                                              <  1:926695/149‑1 (MQ=255)
                        tctGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCa                                                                                                                            <  3:81853/149‑1 (MQ=255)
                        tctGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCa                                                                                                                            <  2:820520/149‑1 (MQ=255)
                              cGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGtt                                                                                                                      <  3:20425/149‑1 (MQ=255)
                              cGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGtt                                                                                                                      <  1:410798/149‑1 (MQ=255)
                               gTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGttt                                                                                                                     <  3:723943/149‑1 (MQ=255)
                               gTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGttt                                                                                                                     <  4:26269/149‑1 (MQ=255)
                                    aaaCAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTg                                                                                                                >  1:790639/1‑149 (MQ=255)
                                      aCAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGgtgt                                                                                                                       <  1:995475/140‑1 (MQ=255)
                                      aCAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGgtgt                                                                                                                       >  2:995475/1‑140 (MQ=255)
                                       cAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGaat                                                                                                             >  4:76329/1‑149 (MQ=255)
                                          gACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGaataat                                                                                                          <  4:429466/149‑1 (MQ=255)
                                            cGATGGCAACGGCTAGTAATGAAAGAACCCAAATCATTTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAAtaataa                                                                                                        <  3:107982/149‑1 (MQ=255)
                                            cGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAAtaataa                                                                                                        <  3:418413/149‑1 (MQ=255)
                                                ggCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATc                                                                                                    <  3:100769/149‑1 (MQ=255)
                                                  cAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAg                                                                                                                           <  3:177633/124‑1 (MQ=255)
                                                  cAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAg                                                                                                                           >  4:177633/1‑124 (MQ=255)
                                                   aaCGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACggg                                                                                                                                           <  3:631524/107‑1 (MQ=255)
                                                   aaCGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACggg                                                                                                                                           >  4:631524/1‑107 (MQ=255)
                                                       gCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCAttt                                                                                             <  1:946013/149‑1 (MQ=255)
                                                        ccAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTg                                                                                            >  2:50465/1‑149 (MQ=255)
                                                           gTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGTCTCGTCGGGGCtttt                                                                                                                                     <  3:702043/105‑4 (MQ=255)
                                                           gTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGTCTCGTCGGGGCtttt                                                                                                                                     >  4:702043/1‑102 (MQ=255)
                                                             aaTGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGt                                                                                       >  1:419282/1‑149 (MQ=255)
                                                               tGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGttt                                                                                     >  1:3506/1‑149 (MQ=255)
                                                               tGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGttt                                                                                     >  2:735824/1‑149 (MQ=255)
                                                               tGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGttt                                                                                     <  2:3506/149‑1 (MQ=255)
                                                                gAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTg                                                                                    >  2:972277/1‑149 (MQ=255)
                                                                gAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTg                                                                                    >  4:867668/1‑149 (MQ=255)
                                                                gAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTg                                                                                    >  4:145361/1‑149 (MQ=255)
                                                                 aaaGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTgg                                                                                   >  3:845069/1‑149 (MQ=255)
                                                                  aaGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGt                                                                                  >  2:6248/1‑149 (MQ=255)
                                                                    gAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTca                                                                                >  4:211476/1‑149 (MQ=255)
                                                                    gAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTca                                                                                >  1:353709/1‑149 (MQ=255)
                                                                        ccAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACgggg                                                                                                                                          >  1:188371/1‑87 (MQ=255)
                                                                        ccAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACgggg                                                                                                                                          <  2:188371/87‑1 (MQ=255)
                                                                        ccAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGt                                                                            >  4:1126737/1‑149 (MQ=255)
                                                                          aaaTCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTaa                                                                          >  2:513209/1‑149 (MQ=255)
                                                                          aaaTCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTaa                                                                          <  1:50465/149‑1 (MQ=255)
                                                                          aaaTCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTaa                                                                          >  1:604065/1‑149 (MQ=255)
                                                                          aaaTCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTaa                                                                          >  3:844557/1‑149 (MQ=255)
                                                                           aaTCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAAt                                                                         <  1:6248/149‑1 (MQ=255)
                                                                           aaTCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAAt                                                                         <  3:125324/149‑1 (MQ=255)
                                                                            aTCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAAtt                                                                        <  1:513209/149‑1 (MQ=255)
                                                                            aTCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAAtt                                                                        >  1:764704/1‑149 (MQ=255)
                                                                                aTTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCt                                                                    <  2:366222/149‑1 (MQ=255)
                                                                                 ttCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTc                                                                   >  3:955339/1‑149 (MQ=255)
                                                                                 ttCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTc                                                                   <  3:76866/149‑1 (MQ=255)
                                                                                    aCCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAg                                                                >  2:942308/1‑149 (MQ=255)
                                                                                        ttCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACtg                                                            >  3:223957/1‑149 (MQ=255)
                                                                                           acaACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTgg                                                         <  3:76329/149‑1 (MQ=255)
                                                                                            caACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTgga                                                        <  4:845069/149‑1 (MQ=255)
                                                                                            caACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTgga                                                        <  2:790639/149‑1 (MQ=255)
                                                                                                    ccTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCaaa                                                <  4:955339/149‑1 (MQ=255)
                                                                                                      ttGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATc                                              >  2:138290/1‑149 (MQ=255)
                                                                                                        gTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGaa                                                                                                              <  4:309640/83‑1 (MQ=255)
                                                                                                        gTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGaa                                                                                                              >  3:309640/1‑83 (MQ=255)
                                                                                                        gTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGa                                            <  2:353709/149‑1 (MQ=255)
                                                                                                        gTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGa                                            >  2:892112/1‑149 (MQ=255)
                                                                                                         ttAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGAc                                           <  3:867668/149‑1 (MQ=255)
                                                                                                            aTTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTc                                                                                                                    >  3:698342/1‑73 (MQ=255)
                                                                                                            aTTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTc                                                                                                                    <  4:698342/73‑1 (MQ=255)
                                                                                                            aTTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTc                                        >  2:951412/1‑149 (MQ=255)
                                                                                                              tGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCAt                                      <  4:111943/149‑1 (MQ=255)
                                                                                                                 gggggATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTggag                                                       >  1:232291/1‑129 (MQ=255)
                                                                                                                 gggggATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTggag                                                       <  2:232291/129‑1 (MQ=255)
                                                                                                                 gggggATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCa                                   <  3:211476/149‑1 (MQ=255)
                                                                                                                  ggggATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAg                                  >  2:1062545/1‑149 (MQ=255)
                                                                                                                     gATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGgcg                               >  3:591753/1‑149 (MQ=255)
                                                                                                                        gACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATc                                                                                                                                   <  1:184886/46‑1 (MQ=255)
                                                                                                                        gACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATc                                                                                                                                   >  2:184886/1‑46 (MQ=255)
                                                                                                                          cTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGgcgc                              <  3:132988/145‑1 (MQ=255)
                                                                                                                          cTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGgcgc                              >  4:132988/1‑145 (MQ=255)
                                                                                                                           ttGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGAc                                           <  3:967240/131‑1 (MQ=255)
                                                                                                                           ttGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCGGGGt                         >  4:855837/1‑149 (MQ=255)
                                                                                                                               ttCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCGGGGTATCt                     <  2:653451/149‑1 (MQ=255)
                                                                                                                                    ccATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCGGGGTATCTGCCgg                >  3:661198/1‑149 (MQ=255)
                                                                                                                                       tAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCgg                            >  2:1134552/1‑134 (MQ=255)
                                                                                                                                       tAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCgg                            <  1:1134552/134‑1 (MQ=255)
                                                                                                                                       tAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCGGGGTATCTGCCGGAGc             <  4:617675/149‑1 (MQ=255)
                                                                                                                                       tAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCGGGGTATCTGCCGGAGc             <  2:1073450/149‑1 (MQ=255)
                                                                                                                                       tAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCGGGGTATCTGCCGGAGc             <  4:223957/149‑1 (MQ=255)
                                                                                                                                        aaaaaaGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCGGGGTATCTGCCGGAGCg            >  2:437114/1‑149 (MQ=255)
                                                                                                                                          tttaGCCCCGTCTCGTCGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCGGGGTATCTGCCGGAGCGtt          >  3:927180/4‑149 (MQ=255)
                                                                                                                                                cccGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCGGGGTATCTGCCGGAGCGTTAACCGc    <  4:661198/149‑1 (MQ=255)
                                                                                                                                                 ccGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCa                                                  >  4:426607/1‑102 (MQ=255)
                                                                                                                                                 ccGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCa                                                  <  3:426607/102‑1 (MQ=255)
                                                                                                                                                 ccGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCGGGGTATCTGCCGGAGCGTTAACCGCa   >  4:183920/1‑149 (MQ=255)
                                                                                                                                                  cgacgaGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCGGGGTATCTGCCGGAGCGTTAACCGCaa  <  4:66266/149‑1 (MQ=255)
                                                                                                                                                  cgacgaGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCGGGGTATCTGCCGGAGCGTTAACCGCaa  <  1:892112/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCGGGGTATCTGCCGGAGCGTTAACCGCAA  >  NC_000913/2411260‑2411554

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: