Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,817,774 A→T 5.2% intergenic (‑271/+10) yqaB ← / ← argQ fructose‑1‑P and 6‑phosphogluconate phosphatase/tRNA‑Arg

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,817,7740AT 5.2% 1034.6 / 24.0 302intergenic (‑271/+10)yqaB/argQfructose‑1‑P and 6‑phosphogluconate phosphatase/tRNA‑Arg
Reads supporting (aligned to +/- strand):  ref base A (104/180);  new base T (7/9);  total (113/189)
Fisher's exact test for biased strand distribution p-value = 6.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 304 total aligned reads displayed.

TTTTGAAGGAAATTTTTTGGAGAAGATGGTGCATCCGGGAGGATGACTCACTGACGTTCGCCCTTCGGGCCGTTGCTGGCGCAACGTTATCCTTCCTGGTGCCTGCGGTTAATACCGCTATCTCTTCGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAG  >  NC_000913/2817625‑2817922
                                                                                                                                                     |                                                                                                                                                    
  ttGCAGGAAATTTTTTGGAGAAGATGGTGCATCCGGGAGGATGACTCACTGACGTTCGCCCTTCGGGCCGTTGCTGGCGCTACGTTATCCTTCCTGGTGCCTGCGGTTAATACCGCTATCTCTTCGATACCTTCATTGCTGAAGGTTACg                                                                                                                                                    <  2:4025758/150‑1 (MQ=255)
      aGGATATTTTTTGGAGAAGATGGTGCATCCGGGAGGATGACTCACTGACGTTCGCCCTTCGGGCCGTTGCTGGCGCAACGTTATCCTTCCTGGTGCCTGCGGTTAATACCGCTATCTCTTCGATACCTTCATTGCTGAAGGTTACGgaga                                                                                                                                                <  1:6532183/150‑1 (MQ=255)
      aGGAAATTTTTTGGAGAAGATGGTGCATCCGGGAGGATGACTCACTGACGTTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGgaga                                                                                                                                                <  1:2843564/150‑1 (MQ=11)
        gAAATTTTTTGGAGAAGATGGTGCATCCGGGAGGATGACTCACTGACGTTCGCCCTTCGGGCCGTTGCTGGCGCAACGTTATCCTTCCTGGTGCCTGCGGTTAATACCGCTATCTCTTCGATACCTTCATTGCTGAAGGTTACGGAGAat                                                                                                                                              >  2:1499430/1‑150 (MQ=255)
        gAAATTTTTTGGAGAAGATGGTGCATCCGGGAGGATGACTCACTGACGTTCGCCCTTCGGGCCGTTGCTGGCGCAACGTTATCCTTCCTGGTGCCTGCGGGTAATACCGCTATCTCTTCGATACATTCATTGCTGAAGGTTACGGAGaaa                                                                                                                                              >  2:640039/1‑149 (MQ=255)
        gAAATTTTTTGGAGAAGATGGTGCATCCGGGAGGATGACTCACTGACGTTCGCCCTTCGGGCCGTTGCTGGCGCAACGTTATACTTACTGGTGCCTGCGGTTAATACCGCTATCTCTTAGATACCATCATTTCTGAAGGTTACGGAGAat                                                                                                                                              >  1:2948396/1‑150 (MQ=25)
        gAAATTTTTTGGAGAAGATGGTGCATCCGGGAGGATGACTCACTGACGTGCGCCCTTCGGGCCGTGGCTGGCGCCACGGTATCCTTCCTCGTGCCAGCGGGTAATACCGCTAGCTCGTCGATACCTTCCTTGCTGCAGCTTGCGGAGCat                                                                                                                                              >  2:6107046/1‑150 (MQ=21)
          aaTTTTTTGGAGAAGATGGTGCATCCGGGAGGATGACTCACTGACGTTCGCCCTTCGGGCCGTTGCTGGCGCAACGTTATCCTTCCTGGTGCCTGCGGTTAATACCGCTATCTCTTCGATACCTTCATTGCTGAAGGTTACGGAGAatat                                                                                                                                            >  2:5249866/1‑150 (MQ=255)
           aTTTTTTGGAGAAGATGGTGCATCCGGGAGGATGACTCACTGACGTTCGCCCTTCGGGCCGTTGCTGGCGCAACGTTATCCTTCCTGGTGCCTGCGGTTAATACCGCTATCTCTTCGATTCCTTCATTGCTGAAGGTTACGGAGAATATg                                                                                                                                           <  1:3011600/150‑1 (MQ=255)
           aTTTTTTGGAGAAGATGGTGCATCCGGGAGGATGACTCACTGACGTTCGCCCTTCGGGCCGTTGCTGGCGCAACGTTATCCTTCCTGGTGCCTGCGGTTAATACCGCTATCTCTTCGATACCTTCATTGCTGAAGGTTACGGAGAATATg                                                                                                                                           <  2:1647627/150‑1 (MQ=255)
           aTTTTTTGGAGAAGATGGTGCATCCGGGAGGATGACTCACTGACGTTCGCCCTTCGGGCCGTTGCTGGCGCAACGTTATCCTTCCTGGTGCCTGCGGTTAATACCGCTATCTCTTCGATACCTTCATTGCTGAAGGTTACGGAGAATATg                                                                                                                                           <  2:2372965/150‑1 (MQ=255)
           aTTTTTTGGAGAAGATGGTGCATCCGGGAGGATGACTCACTGACGTTCGCCCTTCGGGCCGTTGCTGGCGCAACGTTATCCTTCCTGGTGCCTGCGGTTAATACCGCTATCTCTTCGATACCTTCATTGCTGAAGGTTACGGAGAATATg                                                                                                                                           >  1:2596189/1‑150 (MQ=255)
           aTTTGTGGGAGCAGATGGTGCATCCGGGAGGATGACTCACTGACGTTCGCCCTTCGGGCCGTTGCTGGCGCAACGTTATCCTTCCTGGTGCCTGCGGTTAATACCGCTATCTCTTCGATACCTTCATTGCTGAAGGTTACGGAGAATATg                                                                                                                                           <  2:4520521/150‑1 (MQ=255)
            ttttttGGAGAAGATGGTGCATCCGGGAGGATGACTCACTGACGTTCGCCCTTCGGGCCGTTGCTGGCGCAACGTTATCCTTCCTGGTGCCTGCGGTTAATACCGCTATCTCTTCGATACCTTCATTGCTGAAGGTTACGGAGAATATgg                                                                                                                                          <  2:4225196/150‑1 (MQ=255)
                    agaagaTGGTGCATCCGGGAGGATGACTCACTGACGTTCGCCCTTCGGGCCGTTGCTGGCGCAACGTTATCCTTCCTGGTGCCTGCGGTTAATACCGCTATCTCTTCGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCg                                                                                                                                  <  1:2745484/150‑1 (MQ=255)
                     gaagaTGGTGCATCCGGGAGGATGACTCACTGACGTTCGCCCTTCGGGCCGTTGCTGGCGCAACGTTATCCTTCCTGGTGCCTGCGGTTAATACCGCTATCTCTTCGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCgg                                                                                                                                 <  1:5828819/150‑1 (MQ=255)
                     gaagaTGGTGCATCCGGGAGGATGACTCACTGACGTTCGCCCTTCGGGCCGTTGCTGGCGCAACGTTATCCTTCCTGGTGCCTGCGGTTAAGACCGCTATCTCTTCGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCgg                                                                                                                                 <  1:5199568/150‑1 (MQ=255)
                      aagaTGGTGCATCCGGGAGGATGACTCACTGACGTTCGCCATTCGGGCCGTTGCTGGCGCAACGTTATACTTACTTGTTCCTGCGGTTAATACCTCTATATATTCGATACCTTCATTGCGGAAGTTAACGGAGAATATGGGGAATCCggg                                                                                                                                >  2:3562975/1‑150 (MQ=255)
                       agaTGGTGCATCCGGGAGGATGACTCACTGACGTTCGCCCTTCGGGCCGTTGCTGGCGCAACGTTATCCTTCCTGGTGCCTGCGGTTAATACCGCTATCTCTTCGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGgga                                                                                                                               <  2:5730136/150‑1 (MQ=255)
                       agaTGGTGCATCCGGGAGGATGACTCACTGACGTTCGCCCTTCGGGCCGTTGCTGGCGCAACGTTATCCTTCCTGGTGCCTGCGGTTAATACCGCTATCTCTTCGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGgga                                                                                                                               <  1:3376716/150‑1 (MQ=255)
                       agaTGGTGCATCCGGGAGGATGACTCACTGACGTTCGCCCTTCGGGCCGTTGCTGGCGCAACGTTATCCTTCCTGGTGCCTGCGGTTAATACCGCTATCTCTTCGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGgga                                                                                                                               <  2:5597839/150‑1 (MQ=255)
                         aTGGTGCATCCGGGAGGATGACTCACTGACGTTCGCCCTTCGGGCCGTTGCTGGCGCAACGTTATCCTTCCTGGTGCCTGCGGTTAATACCACTATCTCTTCGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGggagg                                                                                                                             >  2:2889576/1‑150 (MQ=255)
                            gTGCATCCGGGAGGATGACTCACTGACGTTCGCCCTTCGGGCCGTTGCTGGCGCAACGTTATCCTTCCTGGTGCCTGCGGTTGATAGCGCTATCTCGTCGATTCGTTCATTGCTGAAGGTTACGGAGTATATGGTGCGTCCGGGAGGAtt                                                                                                                          >  2:4870494/1‑150 (MQ=38)
                               cATCCGGGAGGATGACTCACTGACGTTCGCCCTTCGGGCCGTTGCTGGCGCAACGTTATCCTTCCTGGTGCCTGCGGTTAATACCGCTATCTCTTCGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGa                                                                                                                       <  1:4185490/150‑1 (MQ=255)
                               cATCCGGGAGGATGACTCACTGACGTTCGCCCTTCGGGCCGTTGCTGGCGCAACGTTATCCTTCCTGGTGCCTGCGGTTAATACCGCTATCTCTTCGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGa                                                                                                                       <  2:3242484/150‑1 (MQ=255)
                                       aggaTGACTCACTGACGTTCGCCCTTCGGGCCGTTGCTGGCGCAACGTTATCCTTCCTGGTGCCTGCGGTTAATACCGCTATCTCTTCGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGa                                                                                                               <  1:2438943/150‑1 (MQ=255)
                                        ggaTGACTCACTGACGTTCGCCCTTCGGGCCGTTGCTGGCGCAACGTTATCCTTCCTGGTGCCTGCGGTTAATACCGCTATCTCTTCGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGAc                                                                                                              <  2:2593140/150‑1 (MQ=255)
                                        ggaTGACTCACTGACGTTCGCCCTTCGGGCCGTTGCTGGCGCAACGTTATCCTTCCTGGTGCCTGCGGTTAATACCGCTATCTCTTCGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGAc                                                                                                              <  2:1812140/150‑1 (MQ=255)
                                         gaTGACTCACTGACGTTCGCCCTTCGGGCCGTTGCTGGCGCAACGTTATCCTTCCTGGTGCCTGCGGTTAATACCGCTATCTCTTCGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGAcc                                                                                                             >  2:4704734/1‑150 (MQ=255)
                                             aCTCACTGACGTTCGCCCTTCGGGCCGTTGCTGGCGCAACGTTATCCTTCCTGGTGCCTGCGGTTAATACCGCTATCTCTTCGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTc                                                                                                         <  1:3808938/150‑1 (MQ=255)
                                              cTCACTGACGTTCGCCCTTCGGGCCGTTGCTGGCGCAACGTTATCCTTCCTGGTGCCTGCGGTTAATACCGCTATCTCTTCGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCg                                                                                                        <  2:1150643/150‑1 (MQ=255)
                                              cTCACTGACGTTCGCCCTTCGGGCCGTTGCTGGCGCAACGTTATCCTTCCTGGTGCCTGCGGTTAATACCGCTATCTCTTCGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCg                                                                                                        <  1:5390193/150‑1 (MQ=255)
                                              cTCACTGACGTTCGCCCTTCGGGCCGTTGCTGGCGCAACGTTATCCTTCCTGGTGCCTGCGGTTAATACCGCTATCTCTTCGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCg                                                                                                        <  2:6527486/150‑1 (MQ=255)
                                                  cTGACGTTCGCCCTTCGGGCCGTTGCTGGCGCAACGTTATCCTTCCTGGTGCCTGCGGTTAATACCGCTATCTCTTCGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTc                                                                                                    <  2:6372917/150‑1 (MQ=255)
                                                     aCGTTCGCCCTTCGGGCCGTTGCTGGCGCAACGTTATCCTTCCTGGTGCCTGCGGTTAATACCGCTATCTCTTCGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTa                                                                                                 <  2:4179954/150‑1 (MQ=255)
                                                     aCGTTCGCCCTTCGGGCCGTTGCTGGCGCAACGTTATCCTTCCTGGTGCCTGCGGTTAATACCGCTATCTCTTCGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTa                                                                                                 <  1:1610802/150‑1 (MQ=255)
                                                          cGCCCTTCGGGCCGTTGCTGGCGCAACGTTATCCTTCCTGGTGCCTGCGGTTAATACCGCTATCTCTTCGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGa                                                                                            <  1:6037122/150‑1 (MQ=255)
                                                          cGCCCTTCGGGCCGTTGCTGGCGCAACGTTATCCTTCCTGGTGCCTGCGGTTAATACCGCTATCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGcc                                                                                              >  2:5892390/1‑148 (MQ=255)
                                                             ccTTCGGGCCGTTGCTGGCGCAACGTTATCCTTCCTGGTGCCTGCGGTTAATACCGCTATCTCTTCGATACCTTCATTGCTGAAGGTTTCGGAGAATATGGTGCATCCGgga                                                                                                                               >  2:2822919/1‑112 (MQ=255)
                                                             ccTTCGGGCCGTTGCTGGCGCAACGTTATCCTTCCTGGTGCCTGCGGTTAATACCGCTATCTCTTCGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTa                                                                                         >  2:3102146/1‑150 (MQ=255)
                                                              cTTCGGGCCGTTGCTGGCGCAACGTTATCCTTCCTGGTGCCTGCGGTTAATACCGCTATCTCTTCGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTAc                                                                                        >  2:3357529/1‑150 (MQ=255)
                                                                  gggCCGTTGCTGGCGCAACGTTATCCTTCCTGGTGCCTGCGGTTAATACCGCTATCTCTTCGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTa                                                                                    <  1:5725476/150‑1 (MQ=255)
                                                                  gggCCGTTGCTGGCGCAACGTTATCCTTCCTGGTGCCTGCGGTTAATACCGCTATCTCTTCGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCAGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTa                                                                                    <  2:4692090/150‑1 (MQ=255)
                                                                  gggCCGTTGCTGGCGCAACGTTATCCTTCCTGGTGCCTGCGGTTAATACCGCTATCTCTTCGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCAGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTa                                                                                    <  2:5355302/150‑1 (MQ=255)
                                                                    gCCGTTGCTGGCGCAACGTTATCCTTCCTGGTGCCTGCGGTTAATACCGCTATCTCTTCGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATc                                                                                  <  1:4120920/150‑1 (MQ=255)
                                                                    gCCGTTGCTGGCGCAAAGTTATCCTTCCTGGTGCCTGCGGATAATACCGCTATCTCTTCGATACATTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCTGGAGGATTCGAACCTCCGACAGCTCGGTTCGTAGCCGAGTACTCTATc                                                                                  <  1:1505408/150‑1 (MQ=17)
                                                                       gTTGCTGGCGCAACGTTATCCTTCCTGGTGCCTGCGGTTAATACCGCTATCTCTTCGATACCTTCATTGCTGACGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTTGTAGCCGAGTACTCTATCCAg                                                                               >  2:587689/1‑150 (MQ=255)
                                                                       gTTGCTGGCGCAACGTTATCCTTCCTGGTGCCTGCGGTTAATACCGCTATCTCTTCGATACCTTCATTGCTGAAGTTTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTGCTCTGTCGAg                                                                               >  1:178905/1‑150 (MQ=14)
                                                                       gTTGCTGGCGCAACGTTATCCTTCCTGGTGCCTGCGGTTAATACCGCTATCTCTTCGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGGTCGTAGCCGAGTACTCTATCCAg                                                                               >  2:273441/1‑150 (MQ=255)
                                                                       gTTGCTGGCGCAACGTTATCCTTCCTGGTGCCTGCGGTTAATACCGCTATCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGa                                                                                            >  2:4340038/1‑137 (MQ=255)
                                                                           cTGGCGCAACGTTATCCTTCCTGGTGCCTGCGGTTAATACCGCTATCTCTTCGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGa                                                                           <  1:613330/150‑1 (MQ=255)
                                                                           cTGGCGCAACGTTATCCTTCCTGGTGCCTGCGGTTAATACCGCTATCTCTTCGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGa                                                                           <  1:507208/150‑1 (MQ=255)
                                                                              gcgcAACGTTATCCTTCCTGGTGCCTGCGGTTAATACCGCTATCTCTTCGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTc                                                                                                                         >  1:932559/1‑101 (MQ=255)
                                                                                gcAACGTTATCCTTCCTGGTGCCTGCGGTTAATACCGCTATCTCTTCGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTAc                                                                      >  1:517593/1‑150 (MQ=255)
                                                                                  aaCGTTATCCTTCCTGGTGCCTGCGGTTAATACCGCTATCTCTTCGATAACTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACgg                                                                    >  1:1529035/1‑150 (MQ=255)
                                                                                   aCGTTATCCTTCCTGGTGCCTGCGGTTAATACCGCTATCTCTTCGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGa                                                                   <  2:969609/150‑1 (MQ=255)
                                                                                   aCGTTATCCTTCCTGGTGCCTGCGGTTAATACCGCTATCTCTTCGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGa                                                                   <  1:6505025/150‑1 (MQ=255)
                                                                                        aTCCTTCCTGGTGCCTGCGGTTAATACCGCTATCTCTTCGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCAt                                                              <  1:5249866/150‑1 (MQ=255)
                                                                                          ccttcctGGTGCCTGCGGTTAATACCGCTATCTCTTCGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCg                                                            >  1:776662/1‑150 (MQ=255)
                                                                                           cttcctGGTGCCTGCGGTTAATACCGCTATCTCTTCGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCgg                                                           <  2:2035926/150‑1 (MQ=255)
                                                                                           cttcctGGTGCCTGCGGTTAATACCGCTATCTCTTCGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCgg                                                           <  1:6303289/150‑1 (MQ=255)
                                                                                               ctGGTGCCTGCGGTTAATACCGCTATCTCTTCGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGaaa                                                       <  2:517593/150‑1 (MQ=255)
                                                                                                tGGTGCCTGCGGTTAATACCGCTATCTCTTCGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGaaaa                                                      <  1:5134747/150‑1 (MQ=255)
                                                                                                 ggTGCCTGCGGTTAATACCGCTATCTCTTCGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAAc                                                     <  1:4546589/150‑1 (MQ=255)
                                                                                                 ggTGCCTGCGGTTAATACCGCTATCTCTTCGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAAc                                                     <  2:809710/150‑1 (MQ=255)
                                                                                                    gCCTGCGGTTAATACCGCTATCTCTTCGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTa                                                  <  2:1538932/150‑1 (MQ=255)
                                                                                                     ccTGCGGTTAATACCGCTATCTCTTCGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTAt                                                 >  1:4313568/1‑150 (MQ=255)
                                                                                                      cTGCGGTTAATACCGCTATCTCTTCGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGaaatttattt                                                >  2:5652786/1‑143 (MQ=18)
                                                                                                      cTGCGGTTAATACCGCTATCTCTTCGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGaaatttattt                                                >  2:1179029/1‑143 (MQ=18)
                                                                                                      cTGCGGTTAATACCGCTATCTCTTCGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAActtctt                                                >  1:1127852/1‑150 (MQ=255)
                                                                                                       tGCGGTTAATACCGCTATCTCTTCGATACCTTCATTGCTGAAGGTTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGaaat                                                      <  1:3483133/143‑2 (MQ=18)
                                                                                                         cGGTTAATACCGCTATCTCTTCGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGcc                                                                                              >  2:370411/1‑101 (MQ=255)
                                                                                                          ggTTAATACCGCTATCTCTTCGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGtt                                                                                                     <  2:6271748/93‑1 (MQ=255)
                                                                                                           gTTAATACCGCTATCTCTTCGATACCTTCATTGCTGAAGTTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAg                                                                               <  1:4748491/114‑1 (MQ=255)
                                                                                                           gTTAATACCGCTATCTCTTCGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAg                                                                               <  2:1011991/114‑1 (MQ=255)
                                                                                                           gTTAATACCGCTATCTCTTCGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAg                                                                               >  1:1011991/1‑114 (MQ=255)
                                                                                                           gTTAATACCGCTATCTCTTCGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTAct                                           <  1:2350106/150‑1 (MQ=255)
                                                                                                           gTTAAGACCGCTATCTCTTCGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTAct                                           <  2:2246752/150‑1 (MQ=255)
                                                                                                            ttAATACCGCTATCTCTTCGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATcc                                                                                 >  1:3208543/1‑111 (MQ=255)
                                                                                                              aaTACCGCTATCTCTTCGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTActgct                                        <  1:1384849/150‑1 (MQ=255)
                                                                                                              aaTACCGCTATCTCTTCGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTActgct                                        <  1:4002775/150‑1 (MQ=255)
                                                                                                               aTACCGCTATCTCTTCGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGa                                      >  1:2066256/1‑150 (MQ=2)
                                                                                                                            ttCGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGc                          >  1:4186052/1‑150 (MQ=255)
                                                                                                                              cGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGcta                        <  2:1935556/150‑1 (MQ=255)
                                                                                                                              cGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGcta                        <  2:6002935/150‑1 (MQ=255)
                                                                                                                               gATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTAct                                                                                       >  2:5543105/1‑86 (MQ=38)
                                                                                                                               gATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTAct                                                                                       <  1:5543105/86‑1 (MQ=38)
                                                                                                                                   ccTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTaaa                   <  1:4615163/150‑1 (MQ=255)
                                                                                                                                   ccTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTaaa                   <  2:3140414/150‑1 (MQ=255)
                                                                                                                                   ccTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTaaa                   <  2:2639527/150‑1 (MQ=255)
                                                                                                                                     ttCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTAc                                                                      >  2:1419009/1‑97 (MQ=37)
                                                                                                                                     ttCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCg                                                            >  1:2824910/1‑107 (MQ=35)
                                                                                                                                     ttCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGc                 >  2:6061712/1‑150 (MQ=255)
                                                                                                                                     ttCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGc                 >  1:6355301/1‑150 (MQ=255)
                                                                                                                                      tCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGcc                >  1:4969209/1‑150 (MQ=255)
                                                                                                                                        aTTGCTGAAGGTTACGGAGTATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCat              <  2:3429576/150‑1 (MQ=255)
                                                                                                                                             tGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGtt                                                                                                     <  2:5197828/58‑1 (MQ=255)
                                                                                                                                                 ggTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTaa     <  1:3091963/150‑1 (MQ=255)
                                                                                                                                                   ttACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAaga   <  1:4469228/150‑1 (MQ=255)
                                                                                                                                                   ttACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAaga   <  1:4526478/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TTTTGAAGGAAATTTTTTGGAGAAGATGGTGCATCCGGGAGGATGACTCACTGACGTTCGCCCTTCGGGCCGTTGCTGGCGCAACGTTATCCTTCCTGGTGCCTGCGGTTAATACCGCTATCTCTTCGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAG  >  NC_000913/2817625‑2817922

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