Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,781,067 T→G 10.5% intergenic (+50/‑148) 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,0670TG10.5% 1299.9 / 24.8 431intergenic (+50/‑148)lldD/trmLL‑lactate dehydrogenase, FMN‑linked/tRNA Leu mC34,mU34 2'‑O‑methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base T (137/234);  new base G (56/2);  total (195/236)
Fisher's exact test for biased strand distribution p-value = 1.51e-19
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 432 total aligned reads displayed.

AATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTAT  >  NC_000913/3780920‑3781216
                                                                                                                                                   |                                                                                                                                                     
aaTCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTc                                                                                                                                                     <  2:4982650/150‑1 (MQ=255)
    gATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGgcg                                                                                                                                                 <  1:4117209/150‑1 (MQ=255)
    gATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGgcg                                                                                                                                                 <  1:3007204/150‑1 (MQ=255)
    gATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGgcg                                                                                                                                                 <  1:503187/150‑1 (MQ=255)
     aTCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGCGCCCCCCCACCCAAACCCACACCCGCCGGGCGAGGGGACCCGTCGgcgg                                                                                                                                                >  2:3575628/1‑149 (MQ=255)
     aTCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCAGGGCGAAAGGGAATGCGGCAGAGTCGTTTTCGCCCCCGCGCCCACACCCGCTCCCGCAGGGCGGGGGGGCCGGTCggggg                                                                                                                                                >  2:4798889/1‑147 (MQ=255)
      tCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCt                                                                                                                                               <  2:3366863/150‑1 (MQ=255)
      tCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCt                                                                                                                                               <  2:2668890/150‑1 (MQ=255)
        aGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCGCACCCGAACCCACACCCACAGGGAGAGGGGACCCTTCGGCGCtgt                                                                                                                                             >  1:2411000/1‑150 (MQ=255)
         gCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCtgta                                                                                                                                            <  1:307878/150‑1 (MQ=255)
         gCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCCCACCCAAAACCACACCCACAGGGGGCGGGGGCCGGTCGGCGCagta                                                                                                                                            >  1:4759909/1‑150 (MQ=255)
          cGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGCGCCCCCGCACCCTAAACCACTCCCCCACGGAGGGGGGGCCGTTCGGCGCggttt                                                                                                                                           >  1:662099/1‑148 (MQ=255)
            aaaTTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCAAAACCAAACCCACACGGAGAGGGGACCCTTCGGCGCgggatgt                                                                                                                                         >  2:1877732/1‑150 (MQ=255)
             aaTTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCtgtatgta                                                                                                                                        <  2:5788825/150‑1 (MQ=255)
             aaTTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCtgtatgta                                                                                                                                        <  1:2833658/150‑1 (MQ=255)
             aaTTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCtgtatgta                                                                                                                                        <  2:1587875/150‑1 (MQ=255)
              aTTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCACACCCAAACCCACTCCCACAGGGCGAGGGGACCGTTCGGCGCTGTATGTAc                                                                                                                                       >  2:3875096/1‑150 (MQ=255)
                tACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTc                                                                                                                                     <  2:1487639/150‑1 (MQ=255)
                tACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGCGCCCCCCCACCCAAAACCACACCCACCCGGGGAGGGGGCCCGTCGGCGCGGTGTGTACta                                                                                                                                     >  1:4563455/1‑149 (MQ=255)
                  cGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTccc                                                                                                                                   <  1:3163186/150‑1 (MQ=255)
                  cGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCCCCCCCAAACCCACACCCCCAAGGAGAGGGGGCCCGTCGGCGGggtgtgtaaaacc                                                                                                                                   >  1:4393340/1‑145 (MQ=255)
                   gCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCGCACCCAAAACCACAACCACAGGGAGTGGGGGCCCGTAGGCGCGGTTTTTAATccca                                                                                                                                  >  1:929742/1‑149 (MQ=255)
                   gCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCCCCCCCAAAAACCCCCCCCCACGGAGAGGGGCCCGGGCGGCGCGGGATTTACAcccc                                                                                                                                  >  1:2311496/1‑149 (MQ=255)
                    cAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCCCCCCCAAAACCACACCCGCAGGGAGAGGGGGCCCGTCGGCGCGGTGTGTAAACCCTc                                                                                                                                 >  1:2605943/1‑150 (MQ=255)
                    cAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGCGCCCCCCCCCCCAAAACCACTACCACACGGGGAGGGGGCCCGTGGGCGCGGTATGTAcaaccac                                                                                                                                 >  1:4911609/1‑144 (MQ=255)
                     aaGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCa                                                                                                                                <  2:1323268/150‑1 (MQ=255)
                     aaGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCa                                                                                                                                <  1:5454090/150‑1 (MQ=255)
                       gATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCAct                                                                                                                              <  2:4090614/150‑1 (MQ=255)
                        aTTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCActc                                                                                                                             <  1:4268759/150‑1 (MQ=255)
                           cGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGa                                                                                                                          <  1:5250273/150‑1 (MQ=255)
                             cTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGaaa                                                                                                                        <  2:3292576/150‑1 (MQ=255)
                             cTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGaaa                                                                                                                        <  2:3687226/150‑1 (MQ=255)
                             cTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCAAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGaaa                                                                                                                        <  2:2030724/150‑1 (MQ=255)
                              tGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAAc                                                                                                                       <  1:5396022/150‑1 (MQ=255)
                                      gggCTGGTTAAAGTGTTGCCTGCGGCACTGGCTCCCTTGGCGAATGGGAATGCGGCATTGTGGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCCGGGAGTGGGGACCGTTGGGCGCTGTTTGTACTCCCTCACTCTGAAACGACACCGc                                                                                                               <  1:4482045/150‑1 (MQ=255)
                                         cTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGCGCCCCCACACCCCAACCCCCGACCGCACGGAGAGGGGGCCGGTCGGCGCGGTTTTTTCACCCCCACACTGAAAAGCCAcccccca                                                                                                            >  1:2722873/1‑145 (MQ=255)
                                             gTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGCCGCCCCACACCCAAACCCACGCCCGCAGGGGTAGGGGTCCCTTCGGCGCGGTATTTAATCCCTCACACAGAAACCACACCCCACTCttt                                                                                                        >  2:1257607/1‑150 (MQ=255)
                                               aaaGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCttttt                                                                                                      <  1:4038270/150‑1 (MQ=255)
                                                  gagTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTtc                                                                                                   <  2:525595/150‑1 (MQ=255)
                                                  gagTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCGCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTtc                                                                                                   <  2:2984610/150‑1 (MQ=255)
                                                   agTTGCCTGCGGCACGGGCTCCCATGGCGAAAGGGACTGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTtct                                                                                                  <  1:2819240/150‑1 (MQ=255)
                                                     ttGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTcc                                                                                                <  2:4383031/150‑1 (MQ=255)
                                                     ttGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCCCACCCTAAACCACTCCCTCAGGGCGAGGGGCCCGGTCGGCGCTGTATGTACTCCCCCACTCTGAAACGACACCGCACGCGTTTTTTATcc                                                                                                >  2:4287951/1‑150 (MQ=255)
                                                       gCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCt                                                                                              <  2:5213388/150‑1 (MQ=255)
                                                        ccTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCCTAGTCGTTTGCGCCCCCTCCCCCGAACCCTCTCCCTCCGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTc                                                                                             <  2:2146113/150‑1 (MQ=255)
                                                        ccTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTc                                                                                             <  1:100272/150‑1 (MQ=255)
                                                          tGCGGCACTGGCGCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGc                                                                                           <  2:4015795/150‑1 (MQ=255)
                                                           gCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGcc                                                                                          <  1:5149664/150‑1 (MQ=255)
                                                           gCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGcc                                                                                          <  1:3797699/150‑1 (MQ=255)
                                                           gCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGcc                                                                                          <  2:460346/150‑1 (MQ=255)
                                                            cGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGccc                                                                                         <  1:5720687/150‑1 (MQ=255)
                                                                   ggCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTTTGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCggg                                                                                  <  1:1736245/150‑1 (MQ=255)
                                                                   ggCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCggg                                                                                  <  2:5420105/150‑1 (MQ=255)
                                                                    gCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCgggg                                                                                 <  1:2649448/150‑1 (MQ=255)
                                                                     cTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGga                                                                                <  2:1116283/150‑1 (MQ=255)
                                                                       cccATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGgaga                                                                              <  1:1589074/150‑1 (MQ=255)
                                                                         cATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAgg                                                                            <  1:2105251/150‑1 (MQ=255)
                                                                          aTGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAggg                                                                           <  2:2756246/150‑1 (MQ=255)
                                                                          aTGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCCCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCTCCCCGGAGGGGCGAggg                                                                           <  2:1774425/150‑1 (MQ=255)
                                                                                 aaGGGAATGCGGCATAGTCGTTTTCGCCCCCGCACCCTAAACCACTCCCGCAGGGAGTGGGGGCCCTTGGGCGCGGTTTTTAATGCCGAACACGGAAAAGAAACCCCACAAATTTTTTATACCTCGCCCCACCCGGGGGAGGGGCggggg                                                                    >  2:4713763/1‑149 (MQ=255)
                                                                                  aGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGa                                                                  <  1:3705615/150‑1 (MQ=255)
                                                                                  aGGGAATGCGGCATAGTCGTTTGCGCCCCCGAACCCAAAACCGCTCCCACCGGGGGAGGGGACCGTTCGGCGGTGTATGTACGACCACCCTCAGTAACGACACCGCACAGGTTTTTTCCCCCACGCCCCACCGGGGCGCGGGCCgggggg                                                                   >  2:4830286/1‑148 (MQ=255)
                                                                                   gggAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTTTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGa                                                                  <  2:2100127/150‑1 (MQ=255)
                                                                                   gggAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGa                                                                  <  1:3251521/150‑1 (MQ=255)
                                                                                   gggAATGCGGCATAGTCGTTTGCGCCCCCACACCCAAAACCACACCCACAGGGAGCGGGGACCCGTCGGCGCGGTATTTACTCCCACACTCGGAAACCACACCCCACTCTTTTTTTCTCCCGCGCCCCACCCGGGCGCGGGCCGGGGGGa                                                                  >  2:3731098/1‑150 (MQ=255)
                                                                                   gggAATGCGGCATAGTCGTTAGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGa                                                                  <  2:3470667/150‑1 (MQ=255)
                                                                                    ggAATGCGGCATAGTCGTTGGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAg                                                                 <  1:1360920/150‑1 (MQ=255)
                                                                                      aaTGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAggg                                                               <  2:2585471/150‑1 (MQ=255)
                                                                                           ggaaTAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTTTGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGTGGGCCGGGGTGAGGGGaaaa                                                          <  1:1180384/147‑1 (MQ=255)
                                                                                           ggCATAGTCGTTTGCGCCCCCTCACCCTACCCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGaaaa                                                          <  2:423182/150‑1 (MQ=255)
                                                                                            gCATAGTCGTTTGCGCCCCCCCACCCAAACCCACACCCGCACGGGGAGGGGCCCCGTAGGCGCGGTATGTACTCCCTCACACAGAAACCCAACCCCCCACTTTTTTTCTACCGCCCCCCCCCGGGGCGCGGGCCGGGGGGGCGGGGAAAg                                                         >  2:3175788/1‑150 (MQ=255)
                                                                                                          gcCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCAc                                           >  2:2305225/1‑150 (MQ=255)
                                                                                                            ccccTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGACAAGGCCGCTCTGCTCACCa                                         >  2:731604/1‑150 (MQ=255)
                                                                                                            ccccTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCa                                         >  2:86077/1‑150 (MQ=255)
                                                                                                            ccccTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCa                                         >  2:1475680/1‑150 (MQ=255)
                                                                                                            ccccTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCa                                         >  2:3041069/1‑150 (MQ=255)
                                                                                                              ccTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAAt                                       >  2:2236831/1‑150 (MQ=255)
                                                                                                                  aCCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTc                                   >  1:1890531/1‑150 (MQ=255)
                                                                                                                  aCCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTc                                   >  1:73599/1‑150 (MQ=255)
                                                                                                                  aCCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTc                                   >  1:3599491/1‑150 (MQ=255)
                                                                                                                    ccTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTg                                 >  1:1484391/1‑150 (MQ=255)
                                                                                                                       tccccTCTCCCTGAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTa                              >  2:1078482/3‑150 (MQ=255)
                                                                                                                       aaCCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTa                              >  1:1264289/1‑150 (MQ=255)
                                                                                                                          ccTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATcc                           >  1:1634194/1‑150 (MQ=255)
                                                                                                                           ctctCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCt                          >  2:5336917/1‑150 (MQ=255)
                                                                                                                            tctcCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGTAACGGCCGCACTGCTCACCAATTTTCTGCTATCCTg                         >  1:852994/1‑150 (MQ=255)
                                                                                                                             ctcCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGc                        >  1:5198753/1‑150 (MQ=255)
                                                                                                                             ctcCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGc                        >  1:2930885/1‑150 (MQ=255)
                                                                                                                             ctcCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGc                        >  1:3454669/1‑150 (MQ=255)
                                                                                                                             ctcCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGc                        >  2:4338449/1‑150 (MQ=255)
                                                                                                                                  tgaGGGAGAGGGGACCGTTCGGCGATGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTCTCCTGCCCCCg                   >  1:2185191/3‑150 (MQ=255)
                                                                                                                                  tCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACCCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCg                   >  2:1057518/1‑150 (MQ=255)
                                                                                                                                      ggAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACt               >  2:4491120/1‑150 (MQ=255)
                                                                                                                                        tgcggggACCTTTCGGCTGTTTTTTTACTCCCTCACTCTGAACCTACACCGAACTCTTTTTTTCTCCCTCGCCCGTCGGGGGTGTGGGCGGGGTTGTGGGTAAAAGGTCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTaa             <  2:11996/147‑1 (MQ=255)
                                                                                                                                         gagGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAg            >  1:5580380/1‑150 (MQ=255)
                                                                                                                                          agGGGACTGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAACAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACGAAgg           >  2:5643367/1‑150 (MQ=255)
                                                                                                                                            gggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTCAgggg         >  1:3321658/1‑150 (MQ=255)
                                                                                                                                            gggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAgggg         >  1:423454/1‑150 (MQ=255)
                                                                                                                                            gggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAgggg         >  2:5485571/1‑150 (MQ=255)
                                                                                                                                                   ttCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTAt  >  1:4687333/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
AATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTAT  >  NC_000913/3780920‑3781216

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: