Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,639,349 G→A 6.0% intergenic (+208/+36) pitA → / ← uspB metal phosphate:H(+) symporter PitA/putative universal stress (ethanol tolerance) protein B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,639,3490GA 6.0% 518.0 / 12.3 163intergenic (+208/+36)pitA/uspBmetal phosphate:H(+) symporter PitA/putative universal stress (ethanol tolerance) protein B
Reads supporting (aligned to +/- strand):  ref base G (64/88);  new base A (3/7);  total (68/95)
Fisher's exact test for biased strand distribution p-value = 5.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 168 total aligned reads displayed.

TTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCACCCGCTCACAGCGGCGAATAAATTCATCATCGT  >  NC_000913/3639200‑3639497
                                                                                                                                                     |                                                                                                                                                     
tttGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCATTCTCGAGCGGGTCAGTTTACTgg                                                                                                                                                       <  2:405899/150‑1 (MQ=255)
tttGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTgg                                                                                                                                                       <  2:2266779/150‑1 (MQ=255)
tttGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTgg                                                                                                                                                       <  2:1960554/150‑1 (MQ=255)
 ttGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGc                                                                                                                                                      <  1:929292/150‑1 (MQ=255)
 ttGTTCATGCCGGATGCGGCGTGAACGCCTTATCAGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATAAGTCTCGAGCGGGTCAGTTTACTGGc                                                                                                                                                      <  1:747965/150‑1 (MQ=255)
  tGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGCCAGTAAACTGAcc                                                                                                                                                     <  1:2355891/150‑1 (MQ=255)
   gTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGccc                                                                                                                                                    >  2:700747/1‑150 (MQ=255)
   gTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGccc                                                                                                                                                    <  2:2427523/150‑1 (MQ=255)
   gTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGccc                                                                                                                                                    <  1:1666427/150‑1 (MQ=255)
     tCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGc                                                                                                                                                  >  2:2121799/1‑150 (MQ=255)
      cATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCt                                                                                                                                                 >  2:2450815/1‑150 (MQ=255)
      cATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGCCAGTAAACTGACCCGCt                                                                                                                                                 <  1:320644/150‑1 (MQ=255)
       aTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGCCAGTAAACTGACCCGCtc                                                                                                                                                <  1:716813/150‑2 (MQ=255)
       aTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGCCAGTAAACTGACCCGCtc                                                                                                                                                <  2:2461450/150‑2 (MQ=255)
        tGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCttt                                                                                                                                               >  1:2424783/1‑150 (MQ=255)
         gCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCtttt                                                                                                                                              <  1:1394611/150‑1 (MQ=255)
         gCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGCCAGTAAACTGACCCGCtcga                                                                                                                                              <  2:1358984/150‑4 (MQ=18)
         gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCtttt                                                                                                                                              >  2:2096747/1‑150 (MQ=21)
           cGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCtttttt                                                                                                                                            <  1:734272/150‑1 (MQ=255)
           cGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGCCAGTAAACTGACCCGCtcgaga                                                                                                                                            <  2:484649/150‑6 (MQ=18)
             gataaggcgTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTa                                                                                                                                          >  1:2448400/6‑150 (MQ=25)
                gcggcgTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGc                                                                                                                                       <  2:2187356/150‑1 (MQ=255)
                gcggcgTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGc                                                                                                                                       <  1:2371182/150‑1 (MQ=255)
                 cggcggGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCg                                                                                                                                      <  1:814035/150‑1 (MQ=255)
                 cggcgTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCg                                                                                                                                      >  1:2194040/1‑150 (MQ=255)
                 cggcgTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCg                                                                                                                                      >  1:2214666/1‑150 (MQ=255)
                  ggcgTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGt                                                                                                                                     >  2:288780/1‑150 (MQ=255)
                   gcgTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGtt                                                                                                                                    >  1:1430546/1‑150 (MQ=255)
                    cgTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGtta                                                                                                                                   >  2:1934816/1‑150 (MQ=255)
                     gTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGttat                                                                                                                                  <  2:1879226/150‑1 (MQ=255)
                       gAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGttattt                                                                                                                                >  2:2038945/1‑149 (MQ=255)
                         aCGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTattatt                                                                                                                              <  2:1562020/150‑1 (MQ=255)
                         aCGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTattatt                                                                                                                              <  1:1203077/150‑1 (MQ=255)
                           gCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGc                                                                                                                            <  2:7058/150‑1 (MQ=255)
                             cTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTg                                                                                                                          >  1:1415396/1‑150 (MQ=255)
                             cTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTACGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTg                                                                                                                          >  1:1226852/1‑150 (MQ=255)
                                aTCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTTCTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAg                                                                                                                       <  2:1674645/150‑1 (MQ=255)
                                aTCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAg                                                                                                                       <  1:1339773/150‑1 (MQ=255)
                                aTCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAg                                                                                                                       <  2:1666498/150‑1 (MQ=255)
                                aTCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAg                                                                                                                       <  2:366996/150‑1 (MQ=255)
                                    ggCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCt                                                                                                                   >  2:882021/1‑150 (MQ=255)
                                     gCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCtt                                                                                                                  <  1:655143/150‑1 (MQ=255)
                                      ccTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTa                                                                                                                 <  1:1130045/150‑1 (MQ=255)
                                          cAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGc                                                                                                             <  1:1180054/150‑1 (MQ=255)
                                            aaaaCATGGAAATTCAATAGATTGCAGAGATTATTTAGGCCTGATTATCGTAGCGCATCAGGCAATTTTGCGTTTGTAATGATTCTCGATCGTGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCa                                                                                                           <  2:2071205/150‑1 (MQ=37)
                                               aCATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAat                                                                                                        <  2:2214666/150‑1 (MQ=255)
                                               aCATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAat                                                                                                        <  2:121084/150‑1 (MQ=255)
                                                 aTGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAatca                                                                                                      <  2:79088/150‑1 (MQ=255)
                                                  tGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAatcat                                                                                                     <  2:665068/150‑1 (MQ=255)
                                                  tGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAatcat                                                                                                     >  1:2376173/1‑150 (MQ=255)
                                                   gCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAatcatc                                                                                                    <  2:2292052/150‑1 (MQ=255)
                                                    cAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAAtcatca                                                                                                   >  1:255593/1‑150 (MQ=255)
                                                        ttCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGc                                                                                               >  1:2116837/1‑150 (MQ=255)
                                                         tCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCa                                                                                              <  2:1119109/150‑1 (MQ=255)
                                                               gATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAgg                                                                                       >  1:734108/1‑150 (MQ=255)
                                                                aTTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGc                                                                                      <  1:2120269/150‑1 (MQ=255)
                                                                 ttGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGTGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCt                                                                                     <  1:1840179/150‑1 (MQ=255)
                                                                   gCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGa                                                                                   <  2:2194040/150‑1 (MQ=255)
                                                                   gCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGa                                                                                   >  1:2044540/1‑150 (MQ=255)
                                                                    cAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGac                                                                                  <  1:466095/150‑1 (MQ=255)
                                                                      gagaTTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGacca                                                                                <  2:1756590/150‑1 (MQ=255)
                                                                         aTTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGaccacca                                                                             >  2:63180/1‑150 (MQ=255)
                                                                              gTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGa                                                                        >  1:631921/1‑150 (MQ=255)
                                                                                aGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGAcc                                                                      <  1:700747/150‑1 (MQ=255)
                                                                                 ggCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACca                                                                     <  2:2446749/150‑1 (MQ=255)
                                                                                 ggCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGGCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACca                                                                     <  1:1309780/150‑1 (MQ=255)
                                                                                  gCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACcac                                                                    >  1:2173810/1‑150 (MQ=255)
                                                                                  gCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACcac                                                                    <  2:582814/150‑1 (MQ=255)
                                                                                   ccTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACcaca                                                                   <  1:2281168/150‑1 (MQ=255)
                                                                                   ccTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACcaca                                                                   >  2:2269716/1‑150 (MQ=255)
                                                                                    cTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGTGCGTGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGGTGACCACCACCAGTCcacac                                                                  <  2:2317781/150‑1 (MQ=255)
                                                                                           gCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGTGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATgcgc                                                           >  1:1594190/1‑150 (MQ=255)
                                                                                              tAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTgg                                                        <  2:1560158/150‑1 (MQ=255)
                                                                                               aGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGt                                                       <  1:1934816/150‑1 (MQ=255)
                                                                                                gcgcATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAgg                                                                                       <  1:1572997/117‑1 (MQ=255)
                                                                                                  gcATCAGGCAATTTTGCGTTTGTTATAAGTCTCGAGCTGGTCAGTTTACTTGCGCGCTTTTTTTATGCGTTATTATTGCTTAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCTCACAATGCGCTGGTCAg                                                    <  2:1853374/150‑1 (MQ=255)
                                                                                                   cATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGa                                                   >  1:1688582/1‑150 (MQ=255)
                                                                                                      cAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATa                                                <  2:1366792/150‑1 (MQ=255)
                                                                                                      cAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATa                                                <  1:282364/150‑1 (MQ=255)
                                                                                                          cAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACt                                            >  2:1590343/1‑150 (MQ=255)
                                                                                                                tGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGAc                                      <  1:1993925/150‑1 (MQ=255)
                                                                                                                 gCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACg                                     <  1:349889/150‑1 (MQ=255)
                                                                                                                  cGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGACAAATCATCAATGCAATCAGGCTGCCCACCACCAGACCCCACAATGCGCTGGTCAGAATAAACTGCCGACGc                                    >  2:1254681/1‑150 (MQ=255)
                                                                                                                    tttGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCAc                                  <  1:941729/150‑1 (MQ=255)
                                                                                                                     ttGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCAcc                                 <  1:672288/150‑1 (MQ=255)
                                                                                                                        tAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCACCCGc                              <  1:1210060/150‑1 (MQ=255)
                                                                                                                         aaTCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCACCCGCt                             <  2:1681590/150‑1 (MQ=255)
                                                                                                                          aTCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCACCCGCTc                            <  1:1090385/150‑1 (MQ=255)
                                                                                                                           tCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCACCCGCTca                           >  1:1495784/1‑150 (MQ=255)
                                                                                                                            cAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCACCCGCTcac                          <  1:1711327/150‑1 (MQ=255)
                                                                                                                              gTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCACCCGCTCACAg                        >  1:2317398/1‑150 (MQ=255)
                                                                                                                                ctcGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACAACACAATGCGCTGGTCAGAATAAACTGCCGACGCACCCGCTCACAgcg                      >  2:337005/1‑150 (MQ=255)
                                                                                                                                  cGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAATTTTTTTTTTAATGATTGGTGTACAACAGATTTCTACACAATGCGCTGGTCAGAATAAACTGCCGACGCACCCGCTCACAgcggc                    <  2:320033/150‑1 (MQ=255)
                                                                                                                                     gCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCACCCGCTCACAGCGGCGaa                 <  1:2198711/150‑1 (MQ=255)
                                                                                                                                        ggTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCACCCGCTCACAGCGGCGAATaa              >  1:213484/1‑150 (MQ=255)
                                                                                                                                             gTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCACCCGCTCACAGCGGCGAATAAATtca         <  2:1430546/150‑1 (MQ=255)
                                                                                                                                             gTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCACCCGCTCACAGCAGCGAATAAATtca         >  2:2146796/1‑150 (MQ=255)
                                                                                                                                                  cTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACCCACCCGCTCACAGCGGCGAATAAATTCatcatc    >  2:2353779/1‑150 (MQ=255)
                                                                                                                                                    ggCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCACCCGCTCACAGCGGCGAATAAATTCATCATCGt  >  2:1161179/1‑150 (MQ=255)
                                                                                                                                                    ggCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCACCCGCTCACAGCGGCGAATAAATTCATCATCGt  >  1:615758/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCACCCGCTCACAGCGGCGAATAAATTCATCATCGT  >  NC_000913/3639200‑3639497

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: