Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,811,219 A→G 6.1% intergenic (‑67/+31) mutM ← / ← rpmG formamidopyrimidine/5‑formyluracil/ 5‑hydroxymethyluracil DNA glycosylase/50S ribosomal subunit protein L33

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,811,2190AG 6.1% 342.4 / 8.8 114intergenic (‑67/+31)mutM/rpmGformamidopyrimidine/5‑formyluracil/ 5‑hydroxymethyluracil DNA glycosylase/50S ribosomal subunit protein L33
Reads supporting (aligned to +/- strand):  ref base A (47/60);  new base G (5/2);  total (52/62)
Fisher's exact test for biased strand distribution p-value = 2.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.82e-01

Only 100 of 114 total aligned reads displayed.

GCATGAAGAATGGTTGCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGGTATAGAAGTGAC  >  NC_000913/3811076‑3811365
                                                                                                                                               |                                                                                                                                                  
gCATGAAGAATGGTTGCACCAACGAGCTGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGAAGCAGACAATATGCGGGCAAGCAGATAACAAAAAACCCCGCCGGACCGAgg                                                                                                                                              >  2:1164704/1‑150 (MQ=255)
gCATGAAGAATGGTTGCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCTGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAgg                                                                                                                                              <  1:678939/150‑1 (MQ=255)
gCATGAACACTGGTTGCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATCTGGGGGCAAGCAGATAACAAAAACCCCCGCCGGAGCGAgg                                                                                                                                              <  2:653503/150‑1 (MQ=255)
     aagaaTGGTTGCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAAATTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATTGGGGCAAGCAGATAACAAAAAACACCTCCGGAGCGAGGttttt                                                                                                                                         >  2:392510/1‑150 (MQ=255)
      agaaTGGTTGCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGtttttt                                                                                                                                        >  1:1321745/1‑150 (MQ=255)
       gaaTGGTTGCACCAACGAGATGCGGTGCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCTCGCTCCGGCGGGGTTTTTTg                                                                                                                                       <  2:1207740/150‑1 (MQ=255)
        aaTGGTTGCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGt                                                                                                                                      <  1:984273/150‑1 (MQ=255)
          tGGTTGCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTa                                                                                                                                    >  2:947346/1‑150 (MQ=255)
           ggTTGCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCATGCAGATAAGAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTAc                                                                                                                                   <  2:706958/150‑1 (MQ=255)
           ggTTGCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTAc                                                                                                                                   <  2:1097771/150‑1 (MQ=255)
            gTTGCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCTCGCTCCGGCGGGGTTTTTTGTTatc                                                                                                                                  <  1:1286445/150‑3 (MQ=39)
              tGCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGCTGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATc                                                                                                                                >  1:1066981/1‑150 (MQ=255)
               gCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCa                                                                                                                               <  2:562007/150‑1 (MQ=255)
                  ccATCGAGATGCGGTTCTATGCCGCGGCGGCTGGATTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATCATCAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAg                                                                                                                            <  1:976332/150‑1 (MQ=255)
                  ccAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAg                                                                                                                            <  1:386979/150‑1 (MQ=255)
                   cAACGAGATGCAGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGc                                                                                                                           <  2:1256878/150‑1 (MQ=255)
                    aaCGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATCGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCACAGCg                                                                                                                          >  2:825989/1‑150 (MQ=255)
                     aCGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCga                                                                                                                         >  2:1197707/1‑150 (MQ=255)
                     aCGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGCCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCga                                                                                                                         <  1:253117/150‑1 (MQ=255)
                      cGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTATTGTTACATCAAAGCgag                                                                                                                        >  2:1323236/1‑150 (MQ=255)
                          aTGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAatt                                                                                                                    <  2:967588/150‑1 (MQ=255)
                              ggTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTAttt                                                                                                                >  2:362454/1‑150 (MQ=255)
                              ggTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTAttt                                                                                                                <  2:632346/150‑1 (MQ=255)
                              ggTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTAttt                                                                                                                >  2:340431/1‑150 (MQ=255)
                              ggTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAAATCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTAttt                                                                                                                <  1:110455/150‑1 (MQ=255)
                               gTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTg                                                                                                               <  2:1014991/150‑1 (MQ=255)
                               gTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTg                                                                                                               <  1:514634/150‑1 (MQ=255)
                                          cggcggCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTGCAGGAATGAACAGATGCAGTGAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGcttc                                                                                                    <  1:308742/150‑1 (MQ=255)
                                            gcggcTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGcttctt                                                                                                  <  1:547689/150‑1 (MQ=255)
                                               gcTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGt                                                                                               <  1:273763/150‑1 (MQ=255)
                                                cTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGTGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTAAATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTa                                                                                              >  1:660148/1‑150 (MQ=255)
                                                 tGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAg                                                                                             <  2:303589/150‑1 (MQ=255)
                                                  ggTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGa                                                                                            >  2:846117/1‑150 (MQ=255)
                                                  ggTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGa                                                                                            >  2:745612/1‑150 (MQ=255)
                                                    tttCAACTTCTGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATc                                                                                          <  1:771323/150‑1 (MQ=255)
                                                    tttCAACTTCGGGTATTTCAGTCATAGCATCTCCATGAATGAACTGATGCAGTCTATATGTGGGCAAGCAGATAACAAAATTCCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGATAATTATTTGATTTTCGCTTCTTTGTAGATc                                                                                          <  1:37741/150‑1 (MQ=255)
                                                    tttCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCATTCAATATTGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTATCGCTTCTTTGTAGATc                                                                                          >  2:689898/1‑150 (MQ=255)
                                                    tttCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACACCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATc                                                                                          >  1:582050/1‑150 (MQ=255)
                                                       cAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACg                                                                                       <  2:662621/150‑1 (MQ=255)
                                                          cTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGATGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGc                                                                                    >  1:892895/1‑150 (MQ=255)
                                                           ttCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGTTGCAGTCAATATGGGGGCAAGCCGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCt                                                                                   <  2:403806/150‑1 (MQ=255)
                                                           ttCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCt                                                                                   <  1:477595/150‑1 (MQ=255)
                                                           ttCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCt                                                                                   <  1:899117/150‑1 (MQ=255)
                                                            tCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTg                                                                                  <  2:406338/150‑1 (MQ=255)
                                                             cGGGTAATTCAGTCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACTAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTgg                                                                                 <  1:871172/150‑1 (MQ=255)
                                                             cGGGTAATTCAGGCATAGCATCTCCAGGAATCAACAGATGCAGTCAATATGGGGGAAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTgg                                                                                 <  1:1039878/150‑1 (MQ=255)
                                                                  aaTTCAGGCATAGCATCTCCAGGAATGATCAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCCCGTGCTGGCGaac                                                                            <  1:494198/150‑1 (MQ=255)
                                                                      cAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCTAGAATTATTTGATTTTCGCTTCTTTATAGATCACGTGCTGGCGAACAACt                                                                        >  2:883009/1‑150 (MQ=255)
                                                                          cATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGAt                                                                    <  2:732594/150‑1 (MQ=255)
                                                                            tAGCATCTCCAGGAATGAACATATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGATCTGGCGAACAACTTGATCg                                                                  >  2:570674/1‑150 (MQ=255)
                                                                             aGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGa                                                                 <  1:667625/150‑1 (MQ=255)
                                                                               cATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAAt                                                               <  1:271569/150‑1 (MQ=255)
                                                                                aTCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCTAACAACTTCATCGAAtt                                                              >  2:1013172/1‑150 (MQ=255)
                                                                                  ctcCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAAtttt                                                            <  1:1324478/150‑1 (MQ=255)
                                                                                  ctcCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAAtttt                                                            <  1:24758/150‑1 (MQ=255)
                                                                                   tcCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAAttttt                                                           >  1:240834/1‑150 (MQ=255)
                                                                                     cAGGAATGAACAGATGCAGTCAATATGGTGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTc                                                         <  1:1340738/150‑1 (MQ=255)
                                                                                     cAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTc                                                         <  1:969945/150‑1 (MQ=255)
                                                                                     cAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTc                                                         <  2:1049510/150‑1 (MQ=255)
                                                                                       ggAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAACACCCCGCCGGAGCGAGTTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAg                                                       <  1:1164704/150‑1 (MQ=255)
                                                                                          aTGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTc                                                    >  1:423636/1‑150 (MQ=255)
                                                                                           tGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTcc                                                   <  2:165808/150‑1 (MQ=255)
                                                                                                 gATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGtttt                                             >  2:1066573/1‑150 (MQ=255)
                                                                                                 gATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATAACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGtttt                                             >  1:1162136/1‑150 (MQ=255)
                                                                                                  aTGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGttttt                                            >  1:1294884/1‑150 (MQ=255)
                                                                                                      aGTCAATATGGGGGCAAGCAGATAACAAAACACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCgg                                        <  2:721631/150‑1 (MQ=255)
                                                                                                      aGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCgg                                        <  2:150817/150‑1 (MQ=255)
                                                                                                        tCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCt                                      <  1:589972/150‑1 (MQ=255)
                                                                                                         cAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCtt                                     >  2:233598/1‑150 (MQ=255)
                                                                                                         cAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCtt                                     >  1:109174/1‑150 (MQ=255)
                                                                                                          aaTATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTa                                    >  2:951919/1‑150 (MQ=255)
                                                                                                            tatGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGt                                  <  1:550759/150‑1 (MQ=255)
                                                                                                                 gggCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTGGCTTGTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGGTTTTGCGGCTGAGTACGtt                             >  2:117716/1‑150 (MQ=255)
                                                                                                                 gggCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGtt                             <  1:947346/150‑1 (MQ=255)
                                                                                                                    cAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGt                          <  2:982446/150‑1 (MQ=255)
                                                                                                                      aGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGttc                        >  2:1141027/1‑150 (MQ=255)
                                                                                                                      aGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTGCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCGGGCTTAGTACGTTTGtta                        >  1:843601/1‑149 (MQ=255)
                                                                                                                       tgatgtAACAAAAAACCTCGCTCCGGCGGGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGttct                       >  1:113795/6‑150 (MQ=38)
                                                                                                                       gCAGATAACAAAAAACCCCGCCGGAGCGAGGTGTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTGGTAGATCACGTGCTGGCGAACTACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGttct                       >  1:33462/1‑150 (MQ=255)
                                                                                                                        cAGTTAAGAAAATACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGATTTTTTTGATTTTCTCTTCTTTGTAGATCACGTGCTGGCGAACTACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGttctt                      <  1:145728/150‑1 (MQ=255)
                                                                                                                        cAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGttctt                      <  1:1075044/150‑1 (MQ=255)
                                                                                                                         aGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGttcttc                     <  2:1321745/150‑1 (MQ=255)
                                                                                                                           gtAACAATAAACCTCGCTCCGGCGGGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGt                   >  1:1246022/2‑150 (MQ=255)
                                                                                                                           gtAACAAAAAACCTCGCTCCGGCGGGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGt                   >  2:265405/2‑150 (MQ=39)
                                                                                                                           aTAACAAAAAACCCCGCCGGAGCGATGTTTTTTGTTACATCAAAGCGAGTATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTTGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGt                   <  1:669368/150‑1 (MQ=255)
                                                                                                                           aTAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGATTTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGt                   >  1:498103/1‑150 (MQ=255)
                                                                                                                            tAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTa                  >  2:473304/1‑150 (MQ=255)
                                                                                                                              aCAAAAAACCTCGCTCCGGCGGGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGt                >  1:934429/1‑150 (MQ=255)
                                                                                                                               cAAAAAACCTCGCTCCGGCGGGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATATTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTACAGTACCAGATTTTCCGGCTTAGTACGTTTGTTCTTCGGAGTg               >  1:536014/1‑150 (MQ=255)
                                                                                                                               cAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGtt               >  2:879007/1‑149 (MQ=255)
                                                                                                                               cAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTg               >  1:319333/1‑150 (MQ=255)
                                                                                                                                aaaaaaCCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTAGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTgg              <  2:157763/150‑1 (MQ=255)
                                                                                                                                 aaaaaCCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGGt             <  1:269432/150‑1 (MQ=255)
                                                                                                                                      ccccGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGGTATAGa        >  1:886238/1‑150 (MQ=255)
                                                                                                                                       cccGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACACCTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGGTTGTTCTTCGTAGTGGTATAGaa       >  1:287201/1‑150 (MQ=255)
                                                                                                                                       cccGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGGTATAGaa       >  1:823693/1‑150 (MQ=255)
                                                                                                                                       cccGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGGTATAGaa       >  2:800943/1‑150 (MQ=255)
                                                                                                                                       cccGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTAGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGGTATAGat       >  2:478406/1‑149 (MQ=255)
                                                                                                                                            cgGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGGTATAGAAGTGAc  <  1:392510/150‑1 (MQ=255)
                                                                                                                                            cgGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGGTATAGAAGTGAc  <  2:1199919/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
GCATGAAGAATGGTTGCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGGTATAGAAGTGAC  >  NC_000913/3811076‑3811365

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