Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,781,038 T→A 15.3% intergenic (+21/‑177) lldD → / → trmL L‑lactate dehydrogenase, FMN‑linked/tRNA Leu mC34,mU34 2'‑O‑methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,781,0380TA15.3% 1318.5 / 47.0 501intergenic (+21/‑177)lldD/trmLL‑lactate dehydrogenase, FMN‑linked/tRNA Leu mC34,mU34 2'‑O‑methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base T (47/336);  new base A (63/2);  total (160/341)
Fisher's exact test for biased strand distribution p-value = 8.97e-44
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 503 total aligned reads displayed.

GAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTC‑TCCCTCAG‑GGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCT  >  NC_000913/3780891‑3781187
                                                                                                                                                   |                                                                                                                                                       
gAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTaa                                                                                                                                                       <  1:5015803/150‑1 (MQ=255)
gAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTaa                                                                                                                                                       <  1:2395724/150‑1 (MQ=255)
     tGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCctc                                                                                                                                                  <  1:752307/150‑1 (MQ=255)
     tGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGCGCCCCCACACCCGAAaccac                                                                                                                                                  >  1:4480437/1‑148 (MQ=255)
      ggCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCctc‑t                                                                                                                                                <  2:4257196/150‑1 (MQ=255)
      ggCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCACACCCAAAaccaca                                                                                                                                                 >  1:2632508/1‑147 (MQ=255)
       gCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGCGCCCCCAAACCCAAACCCac‑tc                                                                                                                                               >  1:5820493/1‑150 (MQ=255)
        cGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTC‑Tcc                                                                                                                                              <  2:4162233/150‑1 (MQ=255)
        cGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGAGTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTC‑Tcc                                                                                                                                              <  2:556556/150‑1 (MQ=255)
         gATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCCCCCCCAAACCCCC‑Tccc                                                                                                                                             >  2:4686695/1‑150 (MQ=255)
          aTGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGCGCCCCCACACCCCAAaccacaacca                                                                                                                                             >  2:3985/1‑143 (MQ=255)
            gACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTC‑TCCCTCa                                                                                                                                          <  1:3063099/150‑1 (MQ=255)
            gACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCGCCCCCGAACCCAC‑Tcccgcc                                                                                                                                          >  2:837596/1‑147 (MQ=255)
             aCGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTC‑TCCCTCAg                                                                                                                                         <  1:2534126/150‑1 (MQ=255)
             aCGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTC‑TCCCTCAg                                                                                                                                         <  1:2000004/150‑1 (MQ=255)
             aCGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTC‑TCCCTCAg                                                                                                                                         <  2:5260128/150‑1 (MQ=255)
             aCGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTC‑TCCCTCAg                                                                                                                                         <  1:4653719/150‑1 (MQ=255)
             aCGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTC‑TCCCTCAg                                                                                                                                         <  1:4197820/150‑1 (MQ=255)
               gCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGACTGCGGCATAGTCGTGTGCGCCCCCCCCCCCCAACCCAC‑GCCCCCCg‑gg                                                                                                                                      >  1:1279322/1‑150 (MQ=255)
                 tGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCCCACACTAACCCAC‑GCCCACAg‑gggg                                                                                                                                    >  1:2946050/1‑148 (MQ=255)
                   aCTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTC‑TCCCTGAG‑GGagag                                                                                                                                  <  2:427532/150‑1 (MQ=255)
                    cTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTC‑TCCCTGAG‑GGAGAgg                                                                                                                                 <  1:4760798/150‑1 (MQ=255)
                    cTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTC‑TCCCTGAG‑GGAGAgg                                                                                                                                 <  1:4372512/150‑1 (MQ=255)
                    cTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTC‑TCCCTCAG‑GGAGAgg                                                                                                                                 <  1:1593865/150‑1 (MQ=255)
                     tGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTC‑TCCCTCAG‑GGAGAggg                                                                                                                                <  2:2100526/150‑1 (MQ=255)
                     tGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTC‑TCCCTCAG‑GGAGAggg                                                                                                                                <  1:4510444/150‑1 (MQ=255)
                       gcgcGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTC‑TCCCTCAG‑GGAGAGGGGa                                                                                                                              <  2:5102899/150‑1 (MQ=255)
                       gcgcGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGCGCCCCCTCACCCTAACCCTC‑TCCCTCAG‑GGAGAGGGGa                                                                                                                              <  1:2583980/150‑1 (MQ=255)
                       gcgcGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGCGCCCCCCCCCCCCAACCCAC‑ACCCCCAC‑Gggggggggg                                                                                                                              >  1:11996/1‑149 (MQ=255)
                       gcgcGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGCGCCCCCACACCCAAACCCAC‑GCCCGCAG‑GGTGAggggc                                                                                                                              >  1:1855065/1‑149 (MQ=255)
                        cgcgAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCCTGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTC‑TCCCTGAG‑GGAGAgggtta                                                                                                                             <  2:5699367/150‑4 (MQ=255)
                            aaaTCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCAC‑TCCCACAG‑GGAGAGGGGGCCggt                                                                                                                         >  1:5604553/1‑148 (MQ=255)
                             aaTCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTC‑TCCCTCAG‑GGAGAGGGGACCGTTc                                                                                                                        <  2:4982650/150‑1 (MQ=255)
                                 gATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTC‑TCCCTCAG‑GGAGAGGGGACCGTTCGgcg                                                                                                                    <  1:3941500/150‑1 (MQ=255)
                                 gATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTC‑TCCCTCAG‑GGAGAGGGGACCGTTCGgcg                                                                                                                    <  1:3007204/150‑1 (MQ=255)
                                 gATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCAAAACCAC‑AACCACAG‑GGAGAGGGGACCGGTCGgcg                                                                                                                    >  1:3203379/1‑150 (MQ=255)
                                  aTCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGCGCCCCCCCACCCAAACCCAC‑ACCCGCCG‑GGCGAGGGGACCCGTCGgcgg                                                                                                                   >  2:3575628/1‑149 (MQ=255)
                                   tCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTC‑TCCCTCAG‑GGAGAGGGGACCGTTCGGCGCt                                                                                                                  <  1:5170068/150‑1 (MQ=255)
                                   tCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTC‑TCCCTCAG‑GGAGAGGGGACCGTTCGGCGCt                                                                                                                  <  1:1363960/150‑1 (MQ=255)
                                      gCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGCGCCCCCTCAACCGAACCCTC‑TGCCTCAG‑GGGGAGGGGCCCGTTGGGGGCggta                                                                                                               >  2:3388137/1‑150 (MQ=255)
                                      gCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGCGCCCCCACACCCAAAACCAC‑AACCACAG‑GGAGTGGGGAACCTTCGGCGCtgtt                                                                                                               >  2:2863869/1‑149 (MQ=255)
                                       cGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTC‑TCCCTCAG‑GGAGAGGGGACCGTTCGGCGCtgtat                                                                                                              <  2:1681164/150‑1 (MQ=255)
                                       cGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGCGCCCCCTCACCCGAAACCAC‑TCCCACAG‑GGAGAGGGGACCGGTCGGCGCtgtat                                                                                                              >  1:5342335/1‑150 (MQ=255)
                                       cGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGCGCCCCCCCACCCAAAACCAC‑AACCACCG‑GGAGAGGGGACCCGTAGgcgcgggat                                                                                                              >  2:1337293/1‑145 (MQ=255)
                                         aaaTTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTTCGCCCCCACACCCAAAACCAC‑ACCCTCAG‑GGAGAGGGGGCCCTTGGGCGCggtatgt                                                                                                            >  2:1877728/1‑150 (MQ=255)
                                         aaaTTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCACACCCAAACCCAC‑AACCGCAG‑GGAGAGGGGACCCTTAGGCGCgggatgt                                                                                                            >  1:3301171/1‑150 (MQ=255)
                                          aaTTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTC‑TCCCTCAG‑GGAGAGGGGACCGTTCGGCGCtgtatgta                                                                                                           <  2:2276257/150‑1 (MQ=255)
                                          aaTTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTC‑TCCCTCAG‑GGAGAGGGGACCGTTCGGCGCtgtatgta                                                                                                           <  2:1141349/150‑1 (MQ=255)
                                          aaTTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTC‑TCCCTCAG‑GGAGAGGGGACCGTTCGGCGCtgtatgta                                                                                                           <  1:2833466/150‑1 (MQ=255)
                                          aaTTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTC‑TCCCTCAG‑GGAGAGGGGACCGTTCGGCGCtgtatgta                                                                                                           <  2:5788825/150‑1 (MQ=255)
                                           aTTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGCGCCCCCGCCCCCAAAACCAC‑ACCCTCAT‑GGGGAGGGGGCCGGTCGGCGCGGTATGTAc                                                                                                          >  2:3874957/1‑150 (MQ=255)
                                              aCGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCACACCCAAAAACCA‑TCCCGCAA‑GGGGAGGGGGCCCGTCGGCCCtgtatttaaacc                                                                                                       >  1:5463793/1‑146 (MQ=255)
                                               cGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTC‑TCCCTCAG‑GGAGAGGGGACCGTTCGGCGCTGTATGTACTccc                                                                                                      <  1:2768873/150‑1 (MQ=255)
                                               cGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTC‑TCCCTCAG‑GGAGAGGGGACCGTTCGGCGCTGTATGTACTccc                                                                                                      <  1:3162901/150‑1 (MQ=255)
                                                gCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCCCCCCCAAAACCAC‑ACCCAACG‑GGCGAGGGGGACGTTAGGCGCggtatgtaaaacca                                                                                                     >  2:320956/1‑144 (MQ=255)
                                                gCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCCCCCCCAAAAACCC‑CCCCCCAC‑GGAGAGGGGCCCGGGCGGCGCGGGATTTACAcccc                                                                                                     >  1:2311496/1‑149 (MQ=255)
                                                 cAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGGGCCCCCACACACTAAACCAA‑AAACGCAG‑GGGGTGGGGACCGTTAGGCGCTGTATGTTATcccac                                                                                                    >  1:3601950/1‑148 (MQ=255)
                                                 cAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCACACCCAAACCCTC‑AACCACAG‑GGAGGGGGGCCCGGTCGGCGCGGGATGTTCTaccac                                                                                                    >  1:3408393/1‑148 (MQ=255)
                                                  aaGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGTTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCCC‑TCCCTCAG‑GGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCa                                                                                                   <  1:2227094/150‑1 (MQ=255)
                                                    gATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTC‑TCCCTCAG‑GGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCAct                                                                                                 <  2:1979773/150‑1 (MQ=255)
                                                     aTTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTC‑TCCCTCAG‑GGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCActc                                                                                                <  1:911638/150‑1 (MQ=255)
                                                         gCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTC‑TCCCTCAG‑GGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGaa                                                                                            <  1:675215/150‑1 (MQ=255)
                                                          cTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTC‑TCCCTCAG‑GGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGaaa                                                                                           <  1:995822/150‑1 (MQ=255)
                                                          cTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTC‑TCCCTCAG‑GGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGaaa                                                                                           <  1:507306/150‑1 (MQ=255)
                                                                       tGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTC‑TCCCTCAG‑GGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCActc                                                                              <  2:5401303/150‑1 (MQ=255)
                                                                       tGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCACACCCAAAACCAC‑AACCACAC‑GGAGAGGGGACCGGTAGGCGCAGTATGTTATACCACACAATGAAACCACACCCCAcac                                                                              >  1:3505013/1‑148 (MQ=255)
                                                                       ggggTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTC‑TCCCTCAG‑GGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCCCTCTGAAACGACCCCGCActc                                                                              <  2:1523856/149‑1 (MQ=255)
                                                                            aaaGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTC‑TCCCTCAG‑GGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCttttt                                                                         <  2:1180473/150‑1 (MQ=255)
                                                                            aaaGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTC‑TCCCTCAG‑GGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCttttt                                                                         <  2:5446620/150‑1 (MQ=255)
                                                                            aaaGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGGTTGCGCCCCCTCACCCTACCCCTC‑TCCCTCCG‑GGAGAGGGGACCGTTCGGCGCTGTATGTCCTCCCTCACTCTGAAACGACACCGCACTCttttt                                                                         <  2:4792761/150‑1 (MQ=255)
                                                                              agagTTGCCTGCGGCACTGGCTCCCATGGCGAAAGTGAATGAGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTC‑TCCCTCAG‑GGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAAAGACACCGCACTCttttttt                                                                       <  2:3720812/150‑1 (MQ=255)
                                                                              agagGTGCCTGCGGCACTGGCTCCCAGGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTC‑TCCCTCAG‑GGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCttttttt                                                                       <  2:4083014/150‑1 (MQ=255)
                                                                                  ttGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTC‑TCCCTCAG‑GGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTcc                                                                   <  2:4383031/150‑1 (MQ=255)
                                                                                  ttGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTC‑TCCCTCAG‑GGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTcc                                                                   <  2:352525/150‑1 (MQ=255)
                                                                                    gCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTC‑TCCCTCAG‑GGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCt                                                                 <  2:5213388/150‑1 (MQ=255)
                                                                                    gCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTC‑TCCCTCAG‑GGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCt                                                                 <  2:3332392/150‑1 (MQ=255)
                                                                                     ccTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCCTAGTCGTTTGCGCCCCCTCCCCCGAACCCTC‑TCCCTCCG‑GGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTc                                                                <  2:2146113/150‑1 (MQ=255)
                                                                                     ccTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTC‑TCCCTCAG‑GGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTc                                                                <  1:100272/150‑1 (MQ=255)
                                                                                     ccTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGCGCCCCCGCACCCCAACCCAC‑AACCACCC‑GGAGAGGGGGCCCGTCGGCGCGGGGTGTACTCCCGCCCACAGAAAAGACCCCCCCCACgttttttttaccgc                                                                >  1:4237815/1‑143 (MQ=255)
                                                                                         cGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTC‑TCCCTCAG‑GGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGccc                                                            <  1:5720687/150‑1 (MQ=255)
                                                                                         cGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGCGCCCCCCCAACCAAAACCAA‑TCCC‑CAGCGGGGAGGGGCCCGGTCGGCGCGGGAGGTACTCCCGCACACGGAAACCACACCCCGCACATTTTTTCTCCCGCGccc                                                            >  2:3033766/1‑150 (MQ=255)
                                                                                             aCTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCGCACCCGAACCCAA‑AACCACCG‑GGGGGGGGGACCGGTCGGCGCGGTATGTTCACCCGCAATCAGGAACGACACCGCCCACTTTTTTTCTCCCACCCCCCAcc                                                        >  2:2381504/1‑150 (MQ=255)
                                                                                                ggCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTC‑TCCCTCAG‑GGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCggg                                                     <  1:410313/150‑1 (MQ=255)
                                                                                                  cTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTC‑TCCCTCAG‑GGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGga                                                   <  2:1116283/150‑1 (MQ=255)
                                                                                                       aTGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTC‑TCCCTCAG‑GGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCTCCCCGGAGGGGCGAggg                                              <  1:3707608/150‑1 (MQ=255)
                                                                                                            gAAAGGGAATGCGGCATAGTCGTTTTCGCCCCCCCACCCGAACCCAC‑ACCCGCAG‑GGAGAGGGGGCCCGTCGGCGCGGTTTGTTAACCCGCACGCAGAAACCACACCCCACACTTTTTTTATCCCCCGCCCCACCCGGGGGAGGGGCggg                                         >  2:2737414/1‑150 (MQ=255)
                                                                                                            gAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTC‑TCCCTCAG‑GGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCggg                                         <  2:4496903/150‑1 (MQ=255)
                                                                                                               aGGGAATGCGGCATAGTCGTTTGCGCCCCCGCAACCGAACCCAC‑CCCCCCCG‑GGGGGGGGGACCGGTCGGGGCAGTCTGGACTCCCGCACACGGTAACGACACCGCCCGCGTTTTTTCGACCGCCCCCCACCCGGGCGGGGGCCCGGGTg                                      >  2:4830248/1‑150 (MQ=255)
                                                                                                                gggAATGCGGCATAGTCGTTTGCGCCCCCCCCCCCGAACCCAC‑ACCCAAAG‑GGCGAGGGGTCCCGTAGGCGCGGTAATTAATCCCACACAAAGAAAAGACACCCCACAAATTTTTTCTCACACCCCCCACCGGGGGGCGGGGCGGGGTGa                                     >  2:148027/1‑150 (MQ=255)
                                                                                                                 ggAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTC‑TCCCTCAG‑GGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAg                                    <  2:2228173/150‑1 (MQ=255)
                                                                                                                  gAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTC‑TCCCTCAG‑GGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAgg                                   <  2:4607779/150‑1 (MQ=255)
                                                                                                                  gAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTC‑TCCCTCAG‑GGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAgg                                   <  1:5234053/150‑1 (MQ=255)
                                                                                                                  gAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTC‑TCCCTCAG‑GGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAgg                                   <  2:4296036/150‑1 (MQ=255)
                                                                                                                   aaTGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTC‑TCCCTCAG‑GGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAggg                                  <  2:2618427/150‑1 (MQ=255)
                                                                                                                        ggCATAGTCGTTTGCGCCCCCCCACCCAAAACCAC‑AACCAAAC‑GGGGAGGGGGCCCTTAGGCCCTGTATGTAAACCCACACACAGAAAACACACCCCACAATATTTTTAAACCAAGCCCCACCCGGGAGAGGGCCCGGGGGGGGGGgaaa                             >  2:3160438/1‑150 (MQ=255)
                                                                                                                          cATTGTCGTTTGCGCCCCCTCACCCT‑ACCCTCACCCCCCCG‑GGCTAGGTGACCGTTGGGCGTTTTATGTACCCCCGCACTCTGACACGACACCTCACTATTTTTTTCACCCTCGCCGCTCGGGGGTTAGGGCGGGGTTGAGGGGAAAAgg                           <  2:3915958/150‑1 (MQ=255)
                                                                                                                                         ccccTCACCCTAACCCTC‑TCCCTCAG‑GGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGACAAGGCCGCTCTGCTCACCa            >  2:731604/1‑150 (MQ=255)
                                                                                                                                                 ccTAACCCTC‑TCCCTCAG‑GGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTg    >  1:1711306/1‑150 (MQ=255)
                                                                                                                                                 ccTAACCCTC‑TCCCTCAG‑GGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTg    >  1:1484388/1‑150 (MQ=255)
                                                                                                                                                 ccTAACCCTC‑TCCCTCAG‑GGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTg    >  2:1259540/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                       
GAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTC‑TCCCTCAG‑GGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCT  >  NC_000913/3780891‑3781187

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