Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,811,223 A→G 7.0% intergenic (‑71/+27) 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,2230AG 7.0% 508.9 / 21.0 172intergenic (‑71/+27)mutM/rpmGformamidopyrimidine/5‑formyluracil/ 5‑hydroxymethyluracil DNA glycosylase/50S ribosomal subunit protein L33
Reads supporting (aligned to +/- strand):  ref base A (74/86);  new base G (3/9);  total (77/95)
Fisher's exact test for biased strand distribution p-value = 2.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.81e-01

Only 100 of 175 total aligned reads displayed.

CTGCATGAAGAATGGTTGCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACC  >  NC_000913/3811074‑3811372
                                                                                                                                                     |                                                                                                                                                     
  gCATGAAGAATGGGTGCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAgg                                                                                                                                                     <  1:1008458/150‑1 (MQ=255)
       aagaaTTGTTGCACCAACGAGATGCTGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGttttt                                                                                                                                                <  1:1089574/150‑1 (MQ=255)
        agaaTGGTTGCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCTCGCTCCGGCGGGGtttttt                                                                                                                                               <  2:1143085/150‑1 (MQ=255)
        agaaTGGTTGCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGtttttt                                                                                                                                               <  2:1407615/150‑1 (MQ=255)
        agaaTGGTTGCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGtttttt                                                                                                                                               <  2:226743/150‑1 (MQ=255)
         gaaTGGTTGCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCTCGCTCCGGCGGGGTTTTTTg                                                                                                                                              <  1:1258131/150‑1 (MQ=255)
         gaaTGGTTGCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTg                                                                                                                                              <  2:1313760/150‑1 (MQ=255)
         gaaTGGTTGCACCAACGAGATGCGGGTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCTCGCTCCGGCGGGGTTTTTTg                                                                                                                                              <  1:322746/150‑1 (MQ=255)
             ggTTGCACCAACGAGCTGCTGTTCTATGCCGCGGCTGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGTAATGAACAGATGCAGTCAATATGGGGGCAAGCAGAAAACAAAAAACCTCGCTCCGGCGGGGTTTTTTGTTat                                                                                                                                          <  2:1694485/150‑2 (MQ=38)
             ggTTGCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCCGAGCGAGGTTTTTTGTTAc                                                                                                                                          >  1:890763/1‑150 (MQ=255)
              gTTGCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCTCGCTCCGGCGGGGTTTTTTGTTatc                                                                                                                                         <  2:538755/150‑3 (MQ=39)
                tGCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAAGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCTCGCTCCGGCGGGGTTTTTTGTTatctg                                                                                                                                       <  2:241280/150‑5 (MQ=38)
                 gCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCa                                                                                                                                      >  1:627789/1‑150 (MQ=255)
                 gCACCAACGAGATGCGGTGCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCa                                                                                                                                      <  2:1550803/150‑1 (MQ=255)
                  cACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCaa                                                                                                                                     >  1:743433/1‑150 (MQ=255)
                        cGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCgag                                                                                                                               >  2:1588120/1‑150 (MQ=255)
                          agaTGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGaa                                                                                                                             <  1:638738/150‑1 (MQ=255)
                             tGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAatta                                                                                                                          <  1:978427/150‑1 (MQ=255)
                             tGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAatta                                                                                                                          >  1:561090/1‑150 (MQ=255)
                             tGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAatta                                                                                                                          >  2:1373442/1‑150 (MQ=255)
                               cGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAattatt                                                                                                                        <  1:957472/150‑1 (MQ=255)
                                 gTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTg                                                                                                                      >  1:988380/1‑150 (MQ=255)
                                  ttCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGa                                                                                                                     >  2:727520/1‑150 (MQ=255)
                                  ttCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGa                                                                                                                     <  1:1223372/150‑1 (MQ=255)
                                      aTGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCAGCCGTAGCGAGGTTTTTTGTTACATCAAAGCGACAATTATTTGAtttt                                                                                                                 >  2:93641/1‑150 (MQ=255)
                                         ccGCGGCGGCTGGTTTCAACTTCGTGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGc                                                                                                              <  1:1580657/150‑1 (MQ=255)
                                               cggcTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCttt                                                                                                        >  2:1505431/1‑150 (MQ=255)
                                                ggcTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTg                                                                                                       <  2:466802/150‑1 (MQ=255)
                                                ggcTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTg                                                                                                       <  2:729017/150‑1 (MQ=255)
                                                ggcTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTg                                                                                                       <  2:1211505/150‑1 (MQ=255)
                                                  cTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTa                                                                                                     <  2:709235/150‑1 (MQ=255)
                                                  cTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAAAATCAAATAATTCTCGctttgatgta                                                                                                     <  1:202638/150‑8 (MQ=255)
                                                  cTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAAATTGGGGGCAAGCAGATAACAAAAAACCCCGCCGGATCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTa                                                                                                     <  1:475578/150‑1 (MQ=255)
                                                  cTGGTTTCAACTTCGGGTAATTAAGGCATAGCATCTACAGGAATGAACATATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTa                                                                                                     <  1:315860/150‑1 (MQ=255)
                                                    ggTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGa                                                                                                   <  1:469656/150‑1 (MQ=255)
                                                    ggTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGa                                                                                                   >  1:1650118/1‑150 (MQ=255)
                                                    ggTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGa                                                                                                   <  1:1380373/150‑1 (MQ=255)
                                                       ttCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCa                                                                                                <  1:1186653/150‑1 (MQ=255)
                                                           aCTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGcttctt                                                                                                         <  2:652350/137‑1 (MQ=255)
                                                              tCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTg                                                                                         <  2:837379/150‑1 (MQ=255)
                                                               cGGGTAATTCAGGCATAGCATATCCAGGAATGACCAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTgg                                                                                        <  1:93423/150‑1 (MQ=255)
                                                                 ggTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCg                                                                                      >  1:1418018/1‑150 (MQ=255)
                                                                 ggTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCg                                                                                      >  2:559057/1‑150 (MQ=255)
                                                                  gTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGa                                                                                     >  2:1132416/1‑150 (MQ=255)
                                                                       tCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGaacaac                                                                                >  1:1588695/1‑150 (MQ=255)
                                                                         aGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTg                                                                              >  1:1145422/1‑150 (MQ=255)
                                                                          ggCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTgg                                                                             <  2:743433/150‑1 (MQ=255)
                                                                           gCATAGCATCTCCACGAATGAACAGATGCAGTCAATATCGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATATTCGCATCTTTGTAGATCACGTGCTGGCGAACAACTGGa                                                                            >  1:1612933/1‑150 (MQ=255)
                                                                             aTAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATc                                                                          <  2:1422200/150‑1 (MQ=255)
                                                                              tAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCg                                                                         >  1:480948/1‑150 (MQ=255)
                                                                               aGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGa                                                                                                   <  1:1145261/123‑1 (MQ=255)
                                                                               aGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGATCAACTGGATCGa                                                                        >  2:4904/1‑150 (MQ=255)
                                                                                gCATCTCCAGGAATGAATAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACAGGATCGaa                                                                       >  1:450569/1‑150 (MQ=255)
                                                                                    ctcCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGCATCGAAtttt                                                                   >  2:446317/1‑150 (MQ=255)
                                                                                      ccAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAAtttttt                                                                 >  1:1107310/1‑150 (MQ=255)
                                                                                      ccAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAAtttttt                                                                 <  2:1135544/150‑1 (MQ=255)
                                                                                       cAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATATCAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTc                                                                <  2:1280973/150‑1 (MQ=255)
                                                                                         ggAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAg                                                              >  2:1606668/1‑150 (MQ=255)
                                                                                              gAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCa                                                         >  1:1037369/1‑150 (MQ=255)
                                                                                              gAACAGATGCAGTAAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCa                                                         >  1:276938/1‑150 (MQ=255)
                                                                                               aaCAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCat                                                        >  2:674546/1‑149 (MQ=255)
                                                                                               aaCAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCAAGAATTATTTGATTTTCGCTTCTTTGTTCATCACGTGCAGGCGAACAACTGGATCGAATTTTTTCAGTTACAg                                                        >  2:386667/1‑150 (MQ=255)
                                                                                                aCAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGt                                                       >  2:422548/1‑150 (MQ=255)
                                                                                                   gATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGtttt                                                    >  1:553860/1‑150 (MQ=255)
                                                                                                    aTGCAGTCAATATGGTGGCACGCAGAAAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACAGCTAAGCGAGAATTATTTGATTTTCGCTGCTTTGTAGATCACGTGCTGGCGAACAACTGGCTCGAATGTTTTCAGTTCCAGttttt                                                   <  1:289833/150‑1 (MQ=255)
                                                                                                    aTGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGttttt                                                   <  1:391518/150‑1 (MQ=255)
                                                                                                       cAGGCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCg                                                <  2:114228/150‑1 (MQ=255)
                                                                                                            aaTATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTa                                           >  1:1434050/1‑150 (MQ=255)
                                                                                                             atatGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAg                                          <  2:1418018/150‑1 (MQ=255)
                                                                                                                 ggggtCAAGCAGATAACAAAAAACGCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTTGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACg                                      <  1:1259095/150‑1 (MQ=255)
                                                                                                                   gggCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGtt                                    <  2:627789/150‑1 (MQ=255)
                                                                                                                   gggCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGtt                                    <  2:1620647/150‑1 (MQ=255)
                                                                                                                    ggCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGttt                                   <  1:1641195/150‑1 (MQ=255)
                                                                                                                      cAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTa                                        >  1:341901/1‑143 (MQ=255)
                                                                                                                      cAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTa                                        <  2:341901/143‑1 (MQ=255)
                                                                                                                       aaGCAGATAACAAAAAACCCCGCCGGAGCGAGGTGTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGCTCACGTGCTGGCGAACAACTGGATCGAAGTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGtt                                >  1:788764/1‑150 (MQ=255)
                                                                                                                          cAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGttctt                             >  2:337661/1‑150 (MQ=255)
                                                                                                                          cAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGttctt                             <  2:878946/150‑1 (MQ=255)
                                                                                                                           aGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGttcttc                            <  2:173866/150‑1 (MQ=255)
                                                                                                                             aTAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGt                          >  2:366572/1‑150 (MQ=255)
                                                                                                                             aTAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGt                          <  2:1320429/150‑1 (MQ=255)
                                                                                                                              tAACAAAACACCCCGCCGGTGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTa                         <  2:71192/150‑1 (MQ=255)
                                                                                                                              tAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTa                         >  2:1024102/1‑150 (MQ=255)
                                                                                                                               aaCAAAAAACCTCGCTCCGGCGGGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAAttttt                                                                  >  2:1435902/1‑108 (MQ=37)
                                                                                                                               aaCAAAAAACCTCGCTCCGGCGGGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAAttttt                                                                  <  1:1435902/108‑1 (MQ=37)
                                                                                                                               aaCAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAg                        <  2:1637701/150‑1 (MQ=255)
                                                                                                                                aCAAAAAACCGCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGt                       <  2:301614/150‑1 (MQ=255)
                                                                                                                                aCAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGt                       >  2:1161879/1‑150 (MQ=255)
                                                                                                                                 cAAAAAACCTCGCTCCGGCGGGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTg                      >  1:173517/1‑150 (MQ=255)
                                                                                                                                 cAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTg                      >  2:789254/1‑150 (MQ=255)
                                                                                                                                 cAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTg                      >  2:1021713/1‑150 (MQ=255)
                                                                                                                                  aaaaaaCCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTgg                     <  2:289573/150‑1 (MQ=255)
                                                                                                                                    aaaaCCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGGta                   >  1:400414/1‑150 (MQ=255)
                                                                                                                                       aCCCCGCCGGCGCGAGGTTTTTTGTTACATCAAAGCGAGAATTAATTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGGTATAg                <  2:820442/150‑1 (MQ=255)
                                                                                                                                          ccgccgGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGGTATAGAAg             <  2:638503/150‑1 (MQ=255)
                                                                                                                                            gccgGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGGTATAGAAGTg           <  1:619464/150‑1 (MQ=255)
                                                                                                                                               ggAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGGTATAGAAGTGAcc        <  1:1505431/150‑1 (MQ=255)
                                                                                                                                                   ggAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTa    <  1:674546/149‑1 (MQ=255)
                                                                                                                                                     aGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTAcc  <  1:349468/150‑1 (MQ=255)
                                                                                                                                                     aGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTAcc  <  1:844768/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTGCATGAAGAATGGTTGCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACC  >  NC_000913/3811074‑3811372

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