Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,279,861 G→T 5.5% intergenic (+33/‑119) ghxP → / → yjcE guanine/hypoxanthine permease, high affinity; guanine/hypoxanthine:H+ symporter/putative cation/proton antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,279,8610GT 5.5% 351.2 / 7.0 118intergenic (+33/‑119)ghxP/yjcEguanine/hypoxanthine permease, high affinity; guanine/hypoxanthine:H+ symporter/putative cation/proton antiporter
Reads supporting (aligned to +/- strand):  ref base G (57/53);  new base T (2/5);  total (60/58)
Fisher's exact test for biased strand distribution p-value = 2.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.10e-01

Only 100 of 121 total aligned reads displayed.

TTCGCGACTCTGGTGATTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTTCACCATACTGATAAT  >  NC_000913/4279712‑4280008
                                                                                                                                                     |                                                                                                                                                   
ttCGCGACTCTGGTGATTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTGTCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATg                                                                                                                                                     <  2:1111075/150‑1 (MQ=255)
ttCGCGACTCTGGTGATTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATg                                                                                                                                                     <  2:571162/150‑1 (MQ=255)
ttCGCGACTCTGGTGATTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATg                                                                                                                                                     <  1:1398936/150‑1 (MQ=255)
       ctctGGTGATTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGCCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCAACGGATTCTGAAAAAGTGAGGAAATGGCTGccc                                                                                                                                              >  2:407452/1‑150 (MQ=255)
         ctGGTGATTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGt                                                                                                                                            >  2:1436854/1‑150 (MQ=255)
            gTGATTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGtttt                                                                                                                                         >  1:1218379/1‑150 (MQ=255)
             tGATTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGttttt                                                                                                                                        >  2:356440/1‑150 (MQ=255)
             tGATTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGttttt                                                                                                                                        <  2:720245/150‑1 (MQ=255)
              gATTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGCAGCCATTGCCACCCGttttca                                                                                                                                       <  1:1433252/150‑3 (MQ=255)
               aTTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGTATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTAt                                                                                                                                      <  1:1509461/150‑1 (MQ=255)
                ttGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTAtt                                                                                                                                     <  1:822146/150‑1 (MQ=255)
                 ttgtcgtcTGGTTTCCGGCGAATGGCGCAAGTTTAAAATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGGGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTAttt                                                                                                                                    <  2:1499026/148‑1 (MQ=255)
                 tGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGCAGCCATTGCCACCCGttttcagaa                                                                                                                                    <  2:608985/150‑6 (MQ=255)
                   gtcgtcTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTtc                                                                                                                                  >  1:518261/1‑150 (MQ=255)
                   gtcgtcTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTtc                                                                                                                                  <  2:1561367/150‑1 (MQ=255)
                    tcgtcTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTtct                                                                                                                                 <  1:1338065/150‑1 (MQ=255)
                    tcgtcTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTtct                                                                                                                                 >  1:640980/1‑150 (MQ=255)
                     cgtcTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGTCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTtctc                                                                                                                                <  2:1159359/150‑1 (MQ=255)
                      gtcTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTcc                                                                                                                               <  2:1698725/150‑1 (MQ=255)
                       tcTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCg                                                                                                                              >  1:1366690/1‑150 (MQ=255)
                       tcTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCg                                                                                                                              <  1:1687214/150‑1 (MQ=255)
                       tcTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCg                                                                                                                              >  1:1611782/1‑150 (MQ=255)
                          ggTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTGTTTTCTCCGcac                                                                                                                           >  2:1023944/1‑150 (MQ=255)
                            tttACGGCGAATTGAGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGATTTTGACCTTCTATGCGGGTTTCTGGGCTATCTAATCTTTCATCGTATTCTGAAAACGTGTGGCAATGTCTGCCCGTTTTTATTTTCTCCGCACAt                                                                                                                         <  1:521792/150‑1 (MQ=255)
                               ccGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAAGCTTTCAGCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATtgt                                                                                                                      >  1:280640/1‑150 (MQ=255)
                                cGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCGGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTGTATTTTCTCCGcacaatctt                                                                                                                     >  1:107156/1‑145 (MQ=255)
                                  gCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATtgtgtg                                                                                                                   <  2:1223680/150‑1 (MQ=255)
                                  gCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATtgtgtg                                                                                                                   <  2:570525/150‑1 (MQ=255)
                                   cGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGa                                                                                                                  <  2:348753/150‑1 (MQ=255)
                                     aaTGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGc                                                                                                                >  1:608410/1‑150 (MQ=255)
                                     aaTGGCGCAAGTTGAACATCGGTACGGGCGTTATCGCCGTGGCGCTGGGGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGc                                                                                                                <  2:445029/150‑1 (MQ=255)
                                      aTGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTACGCACATTGTGTGAGCt                                                                                                               >  1:1539713/1‑150 (MQ=255)
                                        ggCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCttt                                                                                                             <  1:504152/150‑1 (MQ=255)
                                          cgcAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTg                                                                                                           >  2:1001699/1‑150 (MQ=255)
                                           gcAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGc                                                                                                          <  1:1111716/150‑1 (MQ=255)
                                           gcAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGc                                                                                                          <  1:356440/150‑1 (MQ=255)
                                           gcAAGTTGAACAGCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGc                                                                                                          <  2:76271/150‑1 (MQ=255)
                                           gcAAGTTGAACAGCGGTACGGGCGTTATCGCCGTGGCGCTGGGGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGc                                                                                                          <  1:1379967/150‑1 (MQ=255)
                                            cAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCg                                                                                                         <  2:1346929/150‑1 (MQ=255)
                                            cAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCg                                                                                                         >  1:1028428/1‑150 (MQ=255)
                                                ttGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGatat                                                                                                     >  1:1025746/1‑150 (MQ=255)
                                                    aCATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTg                                                                                                 >  1:632678/1‑150 (MQ=255)
                                                    aCATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTg                                                                                                 <  1:186947/150‑1 (MQ=255)
                                                    aCATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGGTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTg                                                                                                 >  2:867323/1‑150 (MQ=255)
                                                     cATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGa                                                                                                >  1:423179/1‑150 (MQ=255)
                                                      aacGGTACGGTGGTTATCGCCGTGGCGGTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAACCGGGTGGCAATGGCGGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGaa                                                                                               <  2:280640/148‑1 (MQ=255)
                                                             cGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATg                                                                                        <  2:1040747/150‑1 (MQ=255)
                                                               gTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATgag                                                                                      >  1:201701/1‑150 (MQ=255)
                                                                    tATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTc                                                                                 <  1:1001699/150‑1 (MQ=255)
                                                                     aTCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCa                                                                                <  2:1303542/150‑1 (MQ=255)
                                                                       cGCCGTGGCGCTGGTGTCCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAgg                                                                              <  2:283171/150‑1 (MQ=255)
                                                                       cGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTACTCTTTCATCGGATTCTGAACACGGGTGGCTCTGGCTGCCCGTGTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGATAAAAATGAGAATTCAgg                                                                              >  1:252789/1‑150 (MQ=255)
                                                                        gCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGc                                                                             <  2:538018/150‑1 (MQ=255)
                                                                        gCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGc                                                                             >  1:1389276/1‑150 (MQ=255)
                                                                        gCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTACGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGc                                                                             >  2:100187/1‑150 (MQ=255)
                                                                             ggCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAAc                                                                        <  1:1148608/150‑1 (MQ=255)
                                                                               cgcTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGt                                                                      >  1:1375491/1‑150 (MQ=255)
                                                                                gcTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATAGCTGCCCGTTTTTATTTTCTCAGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTc                                                                     >  2:861099/1‑150 (MQ=255)
                                                                                   ggTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCttt                                                                  <  2:1366690/150‑1 (MQ=255)
                                                                                    gtgtCCTTCTTTTCGGGTGGATGTGCTATCTAATATTTCATCGGTTTCTGAAAACTGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGAAGATTGTGTGAGCTTTTGCGATATTCTGTAAAAAATGAGAATTCAGGCATAACGTCTTtc                                                                 <  1:1053255/150‑1 (MQ=255)
                                                                                    gTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTtc                                                                 <  2:255269/150‑1 (MQ=255)
                                                                                     tGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTtct                                                                <  1:1277172/150‑1 (MQ=255)
                                                                                       aCCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTCTGAGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAAAGTCTTTCTCa                                                              >  1:29612/1‑150 (MQ=255)
                                                                                        ccTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAg                                                             >  1:377563/1‑150 (MQ=255)
                                                                                        ccTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGa                                                            <  1:1118101/150‑1 (MQ=255)
                                                                                        ccTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCAGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAg                                                             >  1:935454/1‑150 (MQ=255)
                                                                                        ccTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGAATCTGAAAACGGGTGGCAATGGCTGCCGGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGAGATATTAAGAATAAAATGAGAATTCAGGCATAACGTATTTCTCAg                                                             >  1:112697/1‑150 (MQ=255)
                                                                                         cTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGa                                                            >  2:1412535/1‑150 (MQ=255)
                                                                                          ttCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGAt                                                           <  1:192649/150‑1 (MQ=255)
                                                                                              aTGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATa                                                       >  1:1235928/1‑150 (MQ=255)
                                                                                                 cGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATaaaa                                                    >  1:1004090/1‑150 (MQ=255)
                                                                                                       gCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCTGATCATAAAAAACATc                                              <  2:535080/150‑1 (MQ=255)
                                                                                                       gCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATc                                              >  1:597595/1‑150 (MQ=255)
                                                                                                       gCTGGGATATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATc                                              <  1:100083/150‑1 (MQ=255)
                                                                                                        cTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCg                                             >  1:909975/1‑150 (MQ=255)
                                                                                                         tGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGc                                            >  2:280203/1‑150 (MQ=255)
                                                                                                          tggCTATCTAATCTTTCATAGGATTCTGAAAACGGGTGGCAATTTCTGCCAGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAATACATCGCa                                           <  1:407452/149‑1 (MQ=255)
                                                                                                              tATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACa                                       >  2:863294/1‑150 (MQ=255)
                                                                                                              tATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACa                                       <  2:589182/150‑1 (MQ=255)
                                                                                                              tATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACa                                       >  1:201811/1‑150 (MQ=255)
                                                                                                                  tAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGa                                   <  1:780017/150‑1 (MQ=255)
                                                                                                                  tAATCTATCATCGGATTCTGAAAACGGGAGGCAATGGCTGCCCGTTTTTATTTTATCCTCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAAGGAGAATACAGGCATCACGTCGTTCTCAGATCATAAACAACATCGCAAACAGCGa                                   >  1:263821/1‑150 (MQ=255)
                                                                                                                   aaTCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGaa                                  >  2:631440/1‑150 (MQ=255)
                                                                                                                      cTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGc                               >  2:1058913/1‑150 (MQ=255)
                                                                                                                       tttCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCa                              >  1:864638/1‑150 (MQ=255)
                                                                                                                         tCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATg                            <  2:632678/150‑1 (MQ=255)
                                                                                                                         tCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATg                            >  1:770202/1‑150 (MQ=255)
                                                                                                                         tCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATg                            <  2:423179/150‑1 (MQ=255)
                                                                                                                         tCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATg                            <  1:1436854/150‑1 (MQ=255)
                                                                                                                             cGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGaaa                        <  2:201811/150‑1 (MQ=255)
                                                                                                                                    tGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATcttctt                 >  1:1665822/1‑150 (MQ=255)
                                                                                                                                     aaaaaCGGGCAGCCATTGCCACCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGTGAACGCATGTAAATCttcttc                >  1:1188436/2‑150 (MQ=255)
                                                                                                                                      aaaaCGGGTGGCAATGGCGGCCCGTTTTTATTTTCACCGCACATTGTGTGAGCTTATGCGATATTCTGAAAAAACTGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTTCa               <  1:161385/150‑1 (MQ=255)
                                                                                                                                       aaaCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAAATCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTTCAc              <  2:1130523/150‑1 (MQ=255)
                                                                                                                                        aaCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTTCAcc             >  2:1231453/1‑150 (MQ=255)
                                                                                                                                          cGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTTCACCAt           >  2:154158/1‑150 (MQ=255)
                                                                                                                                             gTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTTCACCATACt        <  2:1388218/150‑1 (MQ=255)
                                                                                                                                                gCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTTCACCATACTGAt     <  1:491293/150‑1 (MQ=255)
                                                                                                                                                 cAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGCACGCATGGAAATCTTCTTCACCATACTGATa    >  2:939582/1‑150 (MQ=255)
                                                                                                                                                   aTGGCTGCCCGTTTTTATTTTCTCCGCACATTGAGTGAGCTTTTGCGATATTCGGAAAAAAATGAGACATCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACTCATGGAAATCTTCTTCACCATACTGATAAt  >  1:35559/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TTCGCGACTCTGGTGATTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTTCACCATACTGATAAT  >  NC_000913/4279712‑4280008

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: