Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,811,212 C→T 7.1% intergenic (‑60/+38) 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,2120CT 7.1% 508.1 / 18.2 171intergenic (‑60/+38)mutM/rpmGformamidopyrimidine/5‑formyluracil/ 5‑hydroxymethyluracil DNA glycosylase/50S ribosomal subunit protein L33
Reads supporting (aligned to +/- strand):  ref base C (77/80);  new base T (3/9);  total (80/91)
Fisher's exact test for biased strand distribution p-value = 1.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.21e-01

Only 100 of 174 total aligned reads displayed.

TTGCGCACCACTGCATGAAGAATGGTTGCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGGTATAGAAG  >  NC_000913/3811064‑3811361
                                                                                                                                                    |                                                                                                                                                     
    gcACCACTGCATGAAGAATGGTTGCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGATTGAACAGATGCAGTCAATATGGGTGCAAGCAGATAACAAAAAACCccgccg                                                                                                                                                  >  1:117889/1‑150 (MQ=255)
     caccacTGCATGAAGAATGGTTGCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCgg                                                                                                                                                 >  2:1374132/1‑150 (MQ=255)
     caccacTGCATGAAGAATGGTTGCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCgg                                                                                                                                                 >  1:301614/1‑150 (MQ=255)
         acTGCAACAAGAATCGTACCACCAACCAGATGCCGTTCTATGCCCCGGACGCTGCAATAAACTTAGGCAAATTCAGGCATAGCAACACCATGAAAGAACAGATGCAGTCAATATCGGGGCAAGCAGATAACAAAAAACCCCGCCGGACCg                                                                                                                                             >  2:1259095/1‑150 (MQ=255)
            gCATGAAGAATGGGTGCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAgg                                                                                                                                          <  1:1008458/150‑1 (MQ=255)
                 aagaaTTGTTGCACCAACGAGATGCTGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGttttt                                                                                                                                     <  1:1089574/150‑1 (MQ=255)
                  agaaTGGTTGCACCAACGAGATGCGGTTCTATGCCTCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCTAGCAGATAACAAAAAACCCCGCCGGAGCGAGGtttttt                                                                                                                                    <  1:1393706/150‑1 (MQ=255)
                  agaaTGGTTGCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCTCGCTCCGGCGGGGtttttt                                                                                                                                    <  2:1143085/150‑1 (MQ=255)
                  agaaTGGTTGCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGtttttt                                                                                                                                    <  2:1407615/150‑1 (MQ=255)
                   gaaTGGTTGCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCTCGCTCCGGCGGGGTTTTTTg                                                                                                                                   <  1:1199541/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)
                        gTTGCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCTCGCTCCGGCGGGGTTTTTTGTTatc                                                                                                                              <  2:538755/150‑3 (MQ=39)
                          tGCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAAGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCTCGCTCCGGCGGGGTTTTTTGTTatctg                                                                                                                            <  2:241280/150‑5 (MQ=38)
                           gCACCAACGAGATGCGGTGCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCa                                                                                                                           <  2:1550803/150‑1 (MQ=255)
                              ccAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCATGCATAGCATCTCCAGGAATGAACATATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACACCTCCGGAGCGAGGTTTTTTTTTACAtataag                                                                                                                        >  2:315860/1‑145 (MQ=255)
                              ccAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAg                                                                                                                        <  1:1628005/150‑1 (MQ=255)
                                    agaTGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGaa                                                                                                                  >  2:1596622/1‑150 (MQ=255)
                                      aTGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAatt                                                                                                                <  1:1267408/150‑1 (MQ=255)
                                       tGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAatta                                                                                                               <  1:978427/150‑1 (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)
                                              cTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGAtt                                                                                                        >  1:628798/1‑150 (MQ=255)
                                               tATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGTCATAGCATCTCCAGGAATGATCAGATGCAGTCAATATGGGGGCAAGCAGTTAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGAttt                                                                                                       <  1:864034/150‑1 (MQ=255)
                                                aTGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCAGCCGTAGCGAGGTTTTTTGTTACATCAAAGCGACAATTATTTGAtttt                                                                                                      >  2:93641/1‑150 (MQ=255)
                                                 tGCCGCGGCGGCTGGTTTCAACTTCTGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGTAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTc                                                                                                     <  1:1644472/150‑1 (MQ=255)
                                                   ccGCGGCGGCTGGTTTCAACTTCGTGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGc                                                                                                   <  1:1580657/150‑1 (MQ=255)
                                                       ggcggcTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAAAATCAAATAATTCTCGctttg                                                                                               <  2:352670/150‑3 (MQ=255)
                                                          ggcTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTg                                                                                            <  2:1211505/150‑1 (MQ=255)
                                                          ggcTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTg                                                                                            <  2:729017/150‑1 (MQ=255)
                                                          ggcTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTg                                                                                            <  2:460484/150‑1 (MQ=255)
                                                            cTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTa                                                                                          <  2:709235/150‑1 (MQ=255)
                                                              ggTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGa                                                                                        <  1:469656/150‑1 (MQ=255)
                                                                tttCAACTTCGGGTAATTCAGGCAAAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATc                                                                                      <  2:294430/150‑1 (MQ=255)
                                                                 ttCAACTTCGGGTTATTCAGGCATAGCATCTCCAGGAATGATCAGATGCAGTCAATTTGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCa                                                                                     <  1:203958/150‑1 (MQ=255)
                                                                  tCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCAc                                                                                    <  1:471585/150‑1 (MQ=255)
                                                                     aCTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGcttctt                                                                                              >  1:652350/1‑137 (MQ=255)
                                                                     aCTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGcttctt                                                                                              <  2:652350/137‑1 (MQ=255)
                                                                        tCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTg                                                                              <  2:837379/150‑1 (MQ=255)
                                                                         cGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTgg                                                                             <  1:926875/150‑1 (MQ=255)
                                                                         cGGGTAATTCAGGCATAGCATATCCAGGAATGACCAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTgg                                                                             <  1:93423/150‑1 (MQ=255)
                                                                          gggTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGc                                                                            <  1:461148/150‑1 (MQ=255)
                                                                           ggTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCg                                                                           >  1:287128/1‑150 (MQ=255)
                                                                           ggTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCg                                                                           >  1:1418018/1‑150 (MQ=255)
                                                                            gTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGa                                                                          >  2:1132416/1‑150 (MQ=255)
                                                                                 tCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGaacaac                                                                     >  1:1588695/1‑150 (MQ=255)
                                                                                   aGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTg                                                                   >  1:347416/1‑150 (MQ=255)
                                                                                   aGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTg                                                                   >  1:1145422/1‑150 (MQ=255)
                                                                                    ggCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTgg                                                                  <  2:743433/150‑1 (MQ=255)
                                                                                      cATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGAt                                                                >  2:1368450/1‑150 (MQ=255)
                                                                                       aTAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATc                                                               <  2:526867/150‑1 (MQ=255)
                                                                                       aTAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATc                                                               <  2:1422200/150‑1 (MQ=255)
                                                                                        tAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCg                                                              >  1:480948/1‑150 (MQ=255)
                                                                                                ccAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAAtttttt                                                      <  2:1135544/150‑1 (MQ=255)
                                                                                                ccAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAAtttttt                                                      >  1:1107310/1‑150 (MQ=255)
                                                                                                 cAGGACTGAACAGATGCAGTCAATATGGGTGCAAGCAGATAACAAAACACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTc                                                     <  2:1610452/150‑1 (MQ=255)
                                                                                                 cAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTc                                                     >  1:690466/1‑150 (MQ=255)
                                                                                                        gAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCa                                              >  1:1037369/1‑150 (MQ=255)
                                                                                                        gAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCa                                              >  2:421812/1‑150 (MQ=255)
                                                                                                         aaCAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCat                                             >  2:674546/1‑149 (MQ=255)
                                                                                                         aaCAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAg                                             >  2:1140963/1‑150 (MQ=255)
                                                                                                          aCAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGt                                            >  2:422548/1‑150 (MQ=255)
                                                                                                             gATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGtttt                                         >  1:553860/1‑150 (MQ=255)
                                                                                                              aTGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGttttt                                        <  1:391518/150‑1 (MQ=255)
                                                                                                               tGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTc                                       <  2:370816/150‑1 (MQ=255)
                                                                                                                gCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTcc                                      <  1:407653/150‑1 (MQ=255)
                                                                                                                gCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTcc                                      <  2:1647595/150‑1 (MQ=255)
                                                                                                                     cAATATGGGGGCAAACAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCtt                                 >  2:1418303/1‑150 (MQ=255)
                                                                                                                      aaTATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTa                                >  1:1434050/1‑150 (MQ=255)
                                                                                                                        tatGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGt                              >  1:610373/1‑150 (MQ=255)
                                                                                                                            ggggCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTTGTACGt                          >  1:1473857/1‑150 (MQ=255)
                                                                                                                            ggggCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGt                          >  2:1699737/1‑150 (MQ=255)
                                                                                                                             gggCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGtt                         <  2:276938/150‑1 (MQ=255)
                                                                                                                                cAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTa                             >  1:341901/1‑143 (MQ=255)
                                                                                                                                 aaGCAGATAACAAAAAACCCCGCCGGAGCGAGGTGTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGCTCACGTGCTGGCGAACAACTGGATCGAAGTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGtt                     >  1:788764/1‑150 (MQ=255)
                                                                                                                                  aGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGttc                    >  2:22482/1‑150 (MQ=255)
                                                                                                                                   gCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGttct                   <  2:569094/150‑1 (MQ=255)
                                                                                                                                    cAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGttctt                  <  2:878946/150‑1 (MQ=255)
                                                                                                                                    cAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGttctt                  >  2:337661/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)
                                                                                                                                         aaCAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAg             <  2:1637701/150‑1 (MQ=255)
                                                                                                                                          aCAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGt            <  1:1132297/150‑1 (MQ=255)
                                                                                                                                          aCAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGt            >  2:1161879/1‑150 (MQ=255)
                                                                                                                                           cAAAAAACCTCGCTCCGGCGGGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTg           >  1:173517/1‑150 (MQ=255)
                                                                                                                                           cAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTg           >  2:1017219/1‑150 (MQ=255)
                                                                                                                                           cAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTg           >  2:789254/1‑150 (MQ=255)
                                                                                                                                            aaaaaaCCTCGCTCCGGCGGGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGATTAGTACGTTTGTTCTTCGTAGTgg          >  1:1664196/1‑150 (MQ=255)
                                                                                                                                            aaaaaaCCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTgg          <  2:471008/150‑1 (MQ=255)
                                                                                                                                            aaaaaaCCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTgg          <  1:556977/150‑1 (MQ=255)
                                                                                                                                             aaaaaCCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGGt         >  2:1579707/1‑150 (MQ=255)
                                                                                                                                              aaaaCCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGGta        >  1:2931/1‑150 (MQ=255)
                                                                                                                                                 aCCCCGCCGGCGCGAGGTTTTTTGTTACATCAAAGCGAGAATTAATTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGGTATAg     <  2:820442/150‑1 (MQ=255)
                                                                                                                                                 aCCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGGTATAg     <  2:1461080/150‑1 (MQ=255)
                                                                                                                                                    ccgccgGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGGTATAGAAg  <  2:638503/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TTGCGCACCACTGCATGAAGAATGGTTGCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGGTATAGAAG  >  NC_000913/3811064‑3811361

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: