Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,411,410 G→T 5.8% intergenic (‑58/+29) 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,4100GT 5.8% 374.6 / 8.3 103intergenic (‑58/+29)yfbS/yfbTputative transporter/sugar phosphatas
Reads supporting (aligned to +/- strand):  ref base G (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 = 1.00e+00

Only 100 of 109 total aligned reads displayed.

ACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCGGGGTATCTGCCGGAGCGTTAACCGCAATGAC  >  NC_000913/2411267‑2411558
                                                                                                                                               |                                                                                                                                                    
  cGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACggg                                                                                                                                               <  2:1090529/149‑1 (MQ=255)
      tCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATTAAAGAACCCAAATCAGTTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGTCTCGTCGGGGCtt                                                                                                                                           <  1:639091/149‑2 (MQ=255)
      tCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTa                                                                                                                                           >  3:617675/1‑149 (MQ=255)
         aTACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAAt                                                                                                                                        >  1:1073450/1‑149 (MQ=255)
          tACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATc                                                                                                                                       >  3:26269/1‑149 (MQ=255)
              cACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGa                                                                                                                                   >  4:100769/1‑149 (MQ=255)
               aCTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGAt                                                                                                                                  <  1:926695/149‑1 (MQ=255)
                 tctGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCa                                                                                                                                <  2:820520/149‑1 (MQ=255)
                 tctGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCa                                                                                                                                <  3:81853/149‑1 (MQ=255)
                       cGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGtt                                                                                                                          <  1:410798/149‑1 (MQ=255)
                       cGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGtt                                                                                                                          <  3:20425/149‑1 (MQ=255)
                        gTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGttt                                                                                                                         <  4:26269/149‑1 (MQ=255)
                        gTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGttt                                                                                                                         <  3:723943/149‑1 (MQ=255)
                             aaaCAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTg                                                                                                                    >  1:790639/1‑149 (MQ=255)
                               aCAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGgtgt                                                                                                                           <  1:995475/140‑1 (MQ=255)
                                cAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGaat                                                                                                                 >  4:76329/1‑149 (MQ=255)
                                 aaGACGCTGGCAACGGACAGTAATGAAAGAACACAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAGGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGaata                                                                                                                >  2:177137/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                                                                                                                                               >  4:631524/1‑107 (MQ=255)
                                            aaCGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACggg                                                                                                                                               <  3:631524/107‑1 (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                                                                                         <  2:3506/149‑1 (MQ=255)
                                                        tGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGttt                                                                                         >  1:3506/1‑149 (MQ=255)
                                                        tGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGttt                                                                                         >  2:735824/1‑149 (MQ=255)
                                                         gAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTg                                                                                        >  4:145361/1‑149 (MQ=255)
                                                         gAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTg                                                                                        >  2:972277/1‑149 (MQ=255)
                                                         gAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTg                                                                                        >  4:867668/1‑149 (MQ=255)
                                                          aaaGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTgg                                                                                       >  3:845069/1‑149 (MQ=255)
                                                           aaGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGt                                                                                      >  2:6248/1‑149 (MQ=255)
                                                             gAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTca                                                                                    >  1:353709/1‑149 (MQ=255)
                                                             gAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTca                                                                                    >  4:211476/1‑149 (MQ=255)
                                                                 ccAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACgggg                                                                                                                                              <  2:188371/87‑1 (MQ=255)
                                                                 ccAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACgggg                                                                                                                                              >  1:188371/1‑87 (MQ=255)
                                                                 ccAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGt                                                                                >  4:1126737/1‑149 (MQ=255)
                                                                   aaaTCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTaa                                                                              >  2:513209/1‑149 (MQ=255)
                                                                   aaaTCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTaa                                                                              >  3:844557/1‑149 (MQ=255)
                                                                   aaaTCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTaa                                                                              <  1:50465/149‑1 (MQ=255)
                                                                    aaTCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAAt                                                                             <  3:125324/149‑1 (MQ=255)
                                                                    aaTCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAAt                                                                             <  1:6248/149‑1 (MQ=255)
                                                                     aTCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAAtt                                                                            >  1:764704/1‑149 (MQ=255)
                                                                     aTCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAAtt                                                                            <  1:513209/149‑1 (MQ=255)
                                                                         aTTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCt                                                                        <  2:366222/149‑1 (MQ=255)
                                                                          ttCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTc                                                                       >  3:955339/1‑149 (MQ=255)
                                                                          ttCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTc                                                                       <  3:76866/149‑1 (MQ=255)
                                                                          ttCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTc                                                                       >  4:559286/1‑149 (MQ=255)
                                                                             aCCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAg                                                                    >  2:942308/1‑149 (MQ=255)
                                                                                 ttCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACtg                                                                >  3:223957/1‑149 (MQ=255)
                                                                                     caACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTgga                                                            <  4:845069/149‑1 (MQ=255)
                                                                                     caACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTgga                                                            <  2:790639/149‑1 (MQ=255)
                                                                                      aaCTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTggag                                                           <  2:764704/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)
                                                                                                 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                                                           <  2:232291/129‑1 (MQ=255)
                                                                                                          gggggATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTggag                                                           >  1:232291/1‑129 (MQ=255)
                                                                                                          gggggATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCa                                       <  3:211476/149‑1 (MQ=255)
                                                                                                           ggggATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAg                                      >  3:66266/1‑149 (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)
                                                                                                                 gACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCgg                                <  3:145361/149‑1 (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                                <  1:1134552/134‑1 (MQ=255)
                                                                                                                                tAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCGGGGTATCTGCCGGAGc                 <  4:223957/149‑1 (MQ=255)
                                                                                                                                tAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCGGGGTATCTGCCGGAGc                 <  2:1073450/149‑1 (MQ=255)
                                                                                                                                tAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCGGGGTATCTGCCGGAGc                 <  4:617675/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)
                                                                                                                                              cgaGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCGGGGTATCTGCCGGAGCGTTAACCGCAatga   <  2:604065/149‑1 (MQ=255)
                                                                                                                                              cgaGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCGGGGTATCTGCCGGAGCGTTAACCGCAatga   <  3:183920/149‑1 (MQ=255)
                                                                                                                                               gagaCGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCGGGGTATCTGCCGGAGCGTTAACCGCAatgac  >  1:265356/1‑149 (MQ=255)
                                                                                                                                               gagaCGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCGGGGTATCTGCCGGAGCGTTAACCGCAatgac  >  3:494419/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
ACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCGGGGTATCTGCCGGAGCGTTAACCGCAATGAC  >  NC_000913/2411267‑2411558

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