Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,781,050 T→A 14.2% intergenic (+33/‑165) 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,0500TA14.2% 851.0 / 23.8 331intergenic (+33/‑165)lldD/trmLL‑lactate dehydrogenase, FMN‑linked/tRNA Leu mC34,mU34 2'‑O‑methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base T (42/213);  new base A (48/2);  total (115/216)
Fisher's exact test for biased strand distribution p-value = 2.58e-28
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 335 total aligned reads displayed.

TGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCG  >  NC_000913/3780902‑3781199
                                                                                                                                                    |                                                                                                                                                     
 gACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCTGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCa                                                                                                                                                     <  2:2949933/150‑1 (MQ=255)
  aCGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAg                                                                                                                                                    <  2:853418/150‑1 (MQ=255)
  aCGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAg                                                                                                                                                    <  1:3656779/150‑1 (MQ=255)
   cGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTGAgg                                                                                                                                                   <  2:1640935/150‑1 (MQ=255)
    gCTGTCTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTCAAGAGTGGCCTGCGGCACTGGCTCCCATGGCGACCGGGAATGCGGCATAGTCGTTTGCGGCCCCTCACCCGAACCCTCTCCCGCgggg                                                                                                                                                  >  1:1298964/1‑150 (MQ=255)
    gCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGCGCCCCCTCACCCTAACCCCCACCCACAcgg                                                                                                                                                  >  2:1009369/1‑150 (MQ=255)
       gACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGgaga                                                                                                                                               <  1:1903195/150‑1 (MQ=255)
       gACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGgaga                                                                                                                                               <  2:1201224/150‑1 (MQ=255)
       gACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGgaga                                                                                                                                               <  2:3781912/150‑1 (MQ=255)
        aCTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGagag                                                                                                                                              <  1:1227587/150‑1 (MQ=255)
         cTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAgg                                                                                                                                             <  1:3411705/150‑1 (MQ=255)
         cTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAgg                                                                                                                                             <  2:1092753/150‑1 (MQ=255)
         cTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAgg                                                                                                                                             <  1:683502/150‑1 (MQ=255)
            gcgcGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGCGCCCCCCCCCCCCAAACCACACCCACAGGGGGAggggg                                                                                                                                          >  1:865426/1‑149 (MQ=255)
             cgcgAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTGAGGGAGAgggtta                                                                                                                                         <  2:3574724/150‑4 (MQ=255)
              gcgAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGAcc                                                                                                                                        <  2:1426867/150‑1 (MQ=255)
                 aaaTCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGtt                                                                                                                                     <  2:791762/150‑1 (MQ=255)
                   aTCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTct                                                                                                                                   <  2:696853/150‑2 (MQ=255)
                   aTCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCACACCCGAACCCACGCCCGCACGGAGGGGGGACCgggcg                                                                                                                                   >  1:3572189/1‑146 (MQ=255)
                    tCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCgg                                                                                                                                  <  1:4040838/150‑1 (MQ=255)
                     cGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGgc                                                                                                                                 <  2:3128033/150‑1 (MQ=255)
                      gATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCCCAACCAAACCCAAAACCCCCGGGAGAGGGGAACgagaggcg                                                                                                                                >  1:3798598/1‑143 (MQ=255)
                      gATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGCGCCCCCACACCCAAAACCACACCCGCAGGGCGAGGGGGCCCATCGgcg                                                                                                                                >  1:1694472/1‑150 (MQ=255)
                       aTCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGgcgc                                                                                                                               <  2:603017/150‑1 (MQ=255)
                        tCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCt                                                                                                                              <  2:1768553/150‑1 (MQ=255)
                        tCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCt                                                                                                                              <  2:3954388/150‑1 (MQ=255)
                        tCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCt                                                                                                                              <  1:3154400/150‑1 (MQ=255)
                        tCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCt                                                                                                                              <  1:1143227/150‑1 (MQ=255)
                         cAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCtg                                                                                                                             <  2:2149295/150‑1 (MQ=255)
                         cAGCGAAATTACGCAAGATTCGCTGGAGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCtg                                                                                                                             <  2:3867786/150‑1 (MQ=255)
                          aGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCGCACCCTAACCCACTCCCGCAGGGGGGGGGGCCCGGTCGgcgcggg                                                                                                                            >  1:3377700/1‑147 (MQ=255)
                           gCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCACACCCCAACCCACACCCACCGGGAGTGGGGACCGTTCGgcgcggga                                                                                                                           >  1:2313271/1‑146 (MQ=255)
                           gCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCACACCCCAAACCCAACCCACCCGGGGTGGGGGCCGGTCGgcgcgggt                                                                                                                           >  2:3331726/1‑146 (MQ=255)
                           gCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCACACCCAAAACCACACCCGCAGGGCGAGGGGGCCGTTCGgcgcggga                                                                                                                           >  1:4175504/1‑146 (MQ=255)
                            cgccaggaCGCACGATTCGCGGGTGCAGGGGCTGGGGAAAGAGGTGCCTGCGGCACTGGCGCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCtgtat                                                                                                                          <  1:823782/143‑1 (MQ=255)
                               aaTTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCAAACCCTCACCCGCAGGGAGAGGGGACCGTTCGGCGCtgtatgta                                                                                                                       >  1:2659416/1‑150 (MQ=255)
                                  tACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTc                                                                                                                    <  1:2727148/150‑1 (MQ=255)
                                    cGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTccc                                                                                                                  <  1:1441457/150‑1 (MQ=255)
                                    cGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTccc                                                                                                                  <  1:520608/150‑1 (MQ=255)
                                     gCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTTCGCCCCCCCACCCAAACCCACACCCAAAGGGAGAGGGGTCCGTTCGGCGCtgtttgtaaaacca                                                                                                                 >  1:572572/1‑144 (MQ=255)
                                     gCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGCGCCCCCACACCCAAACCCGCAACCACACGGAGAGGGGACCCGTACGCGCGGTATGTAATccca                                                                                                                 >  1:1037508/1‑149 (MQ=255)
                                       aaGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCa                                                                                                               <  1:1943626/150‑1 (MQ=255)
                                            tCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAATGGGAATGCGGCATAGTCGTTTGAGCCCCCGCACCCAAACCCACTCCCTACCGGGGAGGGGACCCATAGGCGCGTTATGTACACCCGCActcat                                                                                                          >  2:268224/1‑148 (MQ=255)
                                             cGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGa                                                                                                         <  1:1546089/150‑1 (MQ=255)
                                             cGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGa                                                                                                         <  1:2149338/150‑1 (MQ=255)
                                              gCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCACACCCAAAACCACTACCACAGGGAGAGGGGGACCTTAGGCGCTGTATGTAAAACCACACTCAGaa                                                                                                        >  1:2766189/1‑150 (MQ=255)
                                               cTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGaaa                                                                                                       <  1:303243/150‑1 (MQ=255)
                                                tGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAAc                                                                                                      <  2:166797/150‑1 (MQ=255)
                                                 ggTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGCGCCCCCGCACCCAAAACCAAAACCACAGGGAGAGGGGGCCCTTAGGCGCGGTATTGACTCCCACACAAAGAAACg                                                                                                     >  2:2868531/1‑150 (MQ=255)
                                                  gTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGATAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGa                                                                                                    <  2:2745614/150‑1 (MQ=255)
                                                     cAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGacac                                                                                                 <  2:3743795/150‑1 (MQ=255)
                                                       ggggCTGGGTAATGTGTTGCCTGCGGCACTGGCTCCCTTGGCGAATGGTATTGGGTAATTGGCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGTGAGGGGACCGTTCGGCGCTGTTTGTACTCCCTCACTCTGAAACGACACCg                                                                                               <  1:1611662/150‑1 (MQ=255)
                                                        ggggTGGGTAATGTGTTGCCTGCGGCACTGGCTCCCATGGCTAAAGGTAATGGGGCTTAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCATCTGGGTGTGGGGACCGTTCGGCGCTGTATTTACTCCCTCACTCTGAAACGACACCGc                                                                                              <  1:2561069/150‑1 (MQ=255)
                                                                 aaaGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCttttt                                                                                     <  1:1198690/150‑1 (MQ=255)
                                                                 aaaGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCttttt                                                                                     <  2:4185648/150‑1 (MQ=255)
                                                                 aaaGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCttttt                                                                                     <  2:2067781/150‑1 (MQ=255)
                                                                   agagTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTTTGTACTCCCTCACTCTGAAACGACACCGCACTCttttttt                                                                                   <  1:2622841/150‑1 (MQ=255)
                                                                     agTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGTGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTtct                                                                                 <  1:1836841/150‑1 (MQ=255)
                                                                       ttGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCGCCCCCAAAACCACACACGCCGGGCGGGGGGGCCCGTCGGGGGGGGTTGTTCACCCGCACACTGTAACCACACCCCCCACTTTTTTTGTcc                                                                               >  2:212573/1‑150 (MQ=255)
                                                                        tGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTccc                                                                              <  1:2949068/150‑1 (MQ=255)
                                                                        tGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTccc                                                                              <  1:3486331/150‑1 (MQ=255)
                                                                            tGCGGCACTGGCTCCCATGGCGAAAGGGACTGCGGCATAGTCGTGTGCGCCCCCCCCCCCCAACCCACGCCCACCCGGGGAGGGGCCCCGTCGGCGCGGGATGGACACCCCCACGCGGCAACCACACCCCACGCATTTTTTCTccccccc                                                                          >  2:2336671/1‑146 (MQ=255)
                                                                            tGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGc                                                                          <  1:4073195/150‑1 (MQ=255)
                                                                                  aCTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCACACCCTAAACCACTCCCTCAGGGGGAGGGGACCCTTCGGCGCTGTATGTACTCCCTCACTCTGTAACGACACCGCCCTCGTTTTTTATACCACGCCCCAcc                                                                    >  2:3192320/1‑150 (MQ=255)
                                                                                    tGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCgg                                                                  <  1:2331537/150‑1 (MQ=255)
                                                                                     ggCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCggg                                                                 <  1:2240821/150‑1 (MQ=255)
                                                                                     ggCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCggg                                                                 <  1:1586961/150‑1 (MQ=255)
                                                                                      gCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCGCACCCAAAACCACTCCCCCAGGGGGTGGGGGCCGTTCCGCCCGGGTAGTTAAACCACACTCAGTAAACACAACGCCCTCAATTTTTATCCCTCGCCCCCCCcggg                                                                >  2:563721/1‑150 (MQ=255)
                                                                                      gCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTTCGCCCCCGCACCCAAACCCACACCCGCAGGGCGAGGGGTCCCGTCGGCGCGGTAAGTACAGCCACACGAGGAAACCACACCCCCAACTTTTTTTCTGCCAAGCCCCACCcggg                                                                >  2:2534968/1‑150 (MQ=255)
                                                                                       cTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGga                                                               <  2:2991995/150‑1 (MQ=255)
                                                                                        tCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCCGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGgag                                                              <  1:3971747/150‑1 (MQ=255)
                                                                                          ccctGGCGACAGGGACTGCGGCATCGGCGTTGGCGCCCCCTCCCCCTCCCCCGCGCCCGCCGGGAGAGGGGACCGGGCGGCGCGGTATGTACGCCCTCCCTCTGAACCGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGagag                                                            <  2:43643/147‑1 (MQ=255)
                                                                                           cATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCTCCCCGGAGGGGCGAgg                                                           <  1:2410865/150‑1 (MQ=255)
                                                                                           cATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAgg                                                           <  2:2427831/150‑1 (MQ=255)
                                                                                             tGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGc                                                         <  2:3084167/150‑1 (MQ=255)
                                                                                              ggCGAACGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGcc                                                        <  1:2840860/150‑1 (MQ=255)
                                                                                               gCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCg                                                       <  2:1295234/150‑1 (MQ=255)
                                                                                                 gAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCggg                                                     <  2:4115171/150‑1 (MQ=255)
                                                                                                 gAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCGCACCCGAACCCACAACCACAGGGGGGGGGGGCCGGTCGGGGCGGTATGGTCTCCCACACACAGGAAACAAACCCCCCACTTTTTTTGTACCACGCCCCACCCGGGAGGGGGCCcgg                                                     >  2:3398921/1‑150 (MQ=255)
                                                                                                 gAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCACACCCAAAACCAAAACCCCCCGGAGAGGGGCCCCTTCCGCGCGGGATGTAAGCCCCCACACAGAAAACACACCCCACAAATTTTTTAAACCCACCCCCACCCGGGCGCGGGCCggg                                                     >  2:457144/1‑150 (MQ=255)
                                                                                                   atgggAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCACACTCTGAAACGACACCGCACTCTTTTTTTCCCCCTCGCCCCTCCGGGGAGAGGGCCGGGGt                                                   <  1:1060635/148‑1 (MQ=255)
                                                                                                     gggAATGCGGCATAGTCGTGTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGa                                                 <  2:3264665/150‑1 (MQ=255)
                                                                                                      ggAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAg                                                <  1:2281458/150‑1 (MQ=255)
                                                                                                      ggAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAg                                                <  2:500751/150‑1 (MQ=255)
                                                                                                      ggAATGCGGCATAGTCGTTTGCGCCCCCGCACCCAAAACCAATACCAAAGGGAGAGGGGCCCGTTCGGCGCTGTATGTAATCCCGCACACAGAAACCGCAACCCACACTTTTTTTCAACCACGCCCCCCCGGGGGGAGGGGCggggggtg                                                >  2:187067/1‑146 (MQ=255)
                                                                                                       gAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAgg                                               <  2:2916835/150‑1 (MQ=255)
                                                                                                       gAATGCGGCATAGTCGTGGGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAgg                                               <  1:2355622/150‑1 (MQ=255)
                                                                                                             ggCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCCGGGAGAGGGGACCGTTCGGCGCTGTATGTACGCCCTCCCTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGaaaa                                         <  2:3979909/150‑1 (MQ=255)
                                                                                                             ggCATAGTCGTTTGCGCCCCCACACCCTAACCCACTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGaaaa                                         <  2:3532999/150‑1 (MQ=255)
                                                                                                              gCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAg                                        <  2:2554930/150‑1 (MQ=255)
                                                                                                                            gcCCCCACACCCAAAACCAAACCCACAGGGAGAGGGGACCCGTCGGCGCTGTAAGTACTCCCACACACAGAAAAGACACCCCACACATTTTTTCTACCACGCCCCACCGGGGCGAGGGCCGGGGTGAGGGGAAAAGGCCGCACAGCacac                          >  1:209536/1‑150 (MQ=255)
                                                                                                                                ccTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCCCTGCTCACCAAt                      >  2:1549625/1‑150 (MQ=255)
                                                                                                                                 cTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAAtt                     >  2:4092045/1‑150 (MQ=255)
                                                                                                                                    aCCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTc                  >  1:3239050/1‑150 (MQ=255)
                                                                                                                                        tAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTAGGTACGCCCTCACTCTGAACCGACACCGCACTCTTTGGGTCTCCCTCGCCCCTCCGGGGAGCGGGCCGGGGGGAGGGGACAAGGCCGCACTGCTCACCACTTTTCTGCt              >  2:3622807/1‑150 (MQ=255)
                                                                                                                                           cccTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCTCCCCTCCGGGGAGAGGGCCGGGGTGCGGGGACAAGGCCGAACGGCTAACAAATGTTATGCGATc           >  2:2766991/1‑150 (MQ=255)
                                                                                                                                            ccTCTCCCTCGGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACGCTGAAACGACACCGCACTCTGTTTTTCTCCCTCGCCCCTCCGGGGCGCGGGCCGGGGTGAGGGGACAAGGCCGCACTGCTCACCCCTTTTCTGCTCTcc          >  2:2221833/1‑150 (MQ=255)
                                                                                                                                                tcCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGcc      >  1:4266238/1‑150 (MQ=255)
                                                                                                                                                  ccTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGcccc    >  1:726528/1‑150 (MQ=255)
                                                                                                                                                    tCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCg  >  1:3834688/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCG  >  NC_000913/3780902‑3781199

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