Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,781,022 T→G 7.4% intergenic (+5/‑193) 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,0220TG 7.4% 624.2 / 4.0 212intergenic (+5/‑193)lldD/trmLL‑lactate dehydrogenase, FMN‑linked/tRNA Leu mC34,mU34 2'‑O‑methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base T (46/144);  new base G (19/3);  total (65/147)
Fisher's exact test for biased strand distribution p-value = 1.54e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 216 total aligned reads displayed.

TCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTG  >  NC_000913/3780873‑3781170
                                                                                                                                                     |                                                                                                                                                     
  tGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTg                                                                                                                                                     >  1:834642/1‑150 (MQ=255)
   gATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTgc                                                                                                                                                    <  2:2461627/150‑1 (MQ=255)
    aTCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTgcg                                                                                                                                                   >  2:448973/1‑150 (MQ=255)
    aTCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTgcg                                                                                                                                                   >  1:3226856/1‑150 (MQ=255)
    aTCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTgcg                                                                                                                                                   >  1:1791138/1‑150 (MQ=255)
     tCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTgcgc                                                                                                                                                  <  2:1167410/150‑1 (MQ=255)
       gAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGccc                                                                                                                                                >  1:2135742/1‑150 (MQ=255)
        aaaaaGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGcccc                                                                                                                                               <  2:281864/150‑1 (MQ=255)
          aaaGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCt                                                                                                                                             >  2:695716/1‑150 (MQ=255)
          aaaGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGCGcccccc                                                                                                                                             >  2:2592018/1‑149 (MQ=255)
              agaTGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCAcc                                                                                                                                         >  1:4064998/1‑150 (MQ=255)
                aTGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCGCAccca                                                                                                                                       >  2:2141113/1‑149 (MQ=255)
                aTGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCACAcccc                                                                                                                                       >  2:2141220/1‑149 (MQ=255)
                aTGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAAGTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGccccctccccct                                                                                                                                       <  2:769665/150‑1 (MQ=255)
                    aaGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAAcc                                                                                                                                   <  1:4127851/150‑1 (MQ=255)
                        ggCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCctct                                                                                                                               <  1:2822695/150‑1 (MQ=255)
                         gCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGACTGCGGCATAGTCGTGTGCGCCCCCACACCCGAACCCactc                                                                                                                              >  1:2553918/1‑150 (MQ=255)
                         gCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCGCACCCAAAcccacac                                                                                                                              >  1:4042063/1‑146 (MQ=255)
                           gATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCGCACCCAAACCCACAccc                                                                                                                            >  1:2607090/1‑150 (MQ=255)
                           gATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGCGCCCCCACAACCCAAACCACTacc                                                                                                                            >  1:1485249/1‑150 (MQ=255)
                                cGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTGAgg                                                                                                                       <  1:595775/150‑1 (MQ=255)
                                cGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAgg                                                                                                                       <  1:4015255/150‑1 (MQ=255)
                                 gCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCGCCCCCGAACCCACTCCCACCggg                                                                                                                      >  2:611691/1‑150 (MQ=255)
                                 gCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGCGCCCCCACACCCAAACCCACTCCCGCAcgg                                                                                                                      >  2:1378683/1‑150 (MQ=255)
                                 gCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGCGCCCCCACACCCAAAACCACTACCACAggg                                                                                                                      >  2:1728960/1‑150 (MQ=255)
                                    gACTGGCGCGAAATCGATCAGCGAAAGTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGgaga                                                                                                                   <  2:1187391/150‑1 (MQ=255)
                                     aCTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGagag                                                                                                                  <  2:1541532/150‑1 (MQ=255)
                                        ggCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTGGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAgggg                                                                                                               <  1:2983924/150‑1 (MQ=255)
                                         gcgcGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGa                                                                                                              <  1:3603257/150‑1 (MQ=255)
                                         gcgcGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGa                                                                                                              <  2:521970/150‑1 (MQ=255)
                                            cgAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTGAGGGAGAgggttagg                                                                                                           <  2:2011012/150‑6 (MQ=255)
                                                 tCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCgg                                                                                                      <  2:2550370/150‑1 (MQ=255)
                                                 tCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCgg                                                                                                      <  1:2868348/150‑1 (MQ=255)
                                                 tCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCGGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCgg                                                                                                      <  1:1465170/150‑1 (MQ=255)
                                                  cGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGgc                                                                                                     <  1:283176/150‑1 (MQ=255)
                                                  cGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGgc                                                                                                     <  2:2787601/150‑1 (MQ=255)
                                                  cGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGAAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGgc                                                                                                     <  1:2664215/150‑1 (MQ=255)
                                                      cAGCGAAAGTCCGCAAGATGCGCCGGGGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGCATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGgcgccg                                                                                                 <  1:1925333/150‑3 (MQ=255)
                                                        gCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCtgta                                                                                               <  1:83699/150‑1 (MQ=255)
                                                        gCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCtgta                                                                                               <  2:930792/150‑1 (MQ=255)
                                                          gAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCtgtatg                                                                                             <  2:3699369/150‑1 (MQ=255)
                                                             aTTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTAc                                                                                          <  1:750292/150‑1 (MQ=255)
                                                               tACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTACCCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTc                                                                                        <  2:1791138/150‑1 (MQ=255)
                                                                 cGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTccc                                                                                      >  1:3036320/1‑150 (MQ=255)
                                                                 cGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCCCACCCAAAACCACACCCACCCGGCGGGGGGGCCCTTCGGCCCAGTGTGTTCTccc                                                                                      >  1:1345003/1‑150 (MQ=255)
                                                                 cGCAAGAGTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCCGCGGCACTGGCGCCCATGGCGAAAGGCAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTccc                                                                                      <  2:3198625/150‑1 (MQ=255)
                                                                   cAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCCCACCCAAAACCACACCCACAGGGGGTGGGGACCCATAGGCGCCGGATGTAcaaccac                                                                                    >  2:3247233/1‑144 (MQ=255)
                                                                    aaGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGCGCCCCCGCACCCGAACCCACTCCCGCCGGGGGGGGGGGCCCGTCGGCGCGGTATGTACTCCcaca                                                                                   >  2:3342672/1‑150 (MQ=255)
                                                                       aTTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCActc                                                                                <  1:877801/150‑1 (MQ=255)
                                                                         tCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTg                                                                             <  2:1400407/150‑1 (MQ=255)
                                                                         tCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTg                                                                             <  1:3192630/150‑1 (MQ=255)
                                                                          cGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGa                                                                            <  1:3162611/150‑1 (MQ=255)
                                                                           gCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGaa                                                                           <  1:2598548/150‑1 (MQ=255)
                                                                           gCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGaa                                                                           <  2:2242478/150‑1 (MQ=255)
                                                                            cTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGaaa                                                                          <  1:947679/150‑1 (MQ=255)
                                                                               gTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGa                                                                       <  2:3783883/150‑1 (MQ=255)
                                                                               gTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCGCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGa                                                                       <  1:3560938/150‑1 (MQ=255)
                                                                               gTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGa                                                                       <  2:3019989/150‑1 (MQ=255)
                                                                                   aGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCGCACCCAAACCCACACCCGCACGGTGAGGGGGCCGGTCGGGGGGGTTTGTTCTCCCGCACTCTGGAACGACAcc                                                                   >  2:3664911/1‑150 (MQ=255)
                                                                                   aGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGCCCCCCCCCACCCCAAACCCCACCCACACGGAGAGGGGACCCGGCGCCGCAGGAAGTACACCCCCACACAGAAACCACAcc                                                                   >  2:247954/1‑150 (MQ=255)
                                                                                           ggTAAAGTGTTGCGTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATTGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCTGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCtt                                                           <  1:383370/150‑1 (MQ=255)
                                                                                           ggTAAAGTGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTTAAACGACACCGCACTCtt                                                           <  1:936001/150‑1 (MQ=255)
                                                                                             tAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCtttt                                                         <  1:2370236/150‑1 (MQ=255)
                                                                                             tAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCGCACCCCAAACCACTCCCGCAGGGGGAGGGGGCCGGTCGGCGCGGGTTGTTCACCCGCACACAGAAACCACAACCCACGCtttt                                                         >  1:2247502/1‑150 (MQ=255)
                                                                                                  agTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTtct                                                    <  1:520015/150‑1 (MQ=255)
                                                                                                       ccTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCACACCCAAACCCACACCCCCCGCGCGAGGGGCCCCTTCGGGGCGGGTTGGACTCCCCCACGCTGAAACCTCACCCCACTCATTTTTtcgaccgc                                               >  2:681948/1‑144 (MQ=255)
                                                                                                        cTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCg                                              <  1:2999490/150‑1 (MQ=255)
                                                                                                         tGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGc                                             <  2:405592/150‑1 (MQ=255)
                                                                                                         tGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGc                                             <  1:2467345/150‑1 (MQ=255)
                                                                                                         tGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGc                                             <  2:796264/150‑1 (MQ=255)
                                                                                                         tGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGc                                             <  1:2079217/150‑1 (MQ=255)
                                                                                                         tGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGc                                             <  1:623649/150‑1 (MQ=255)
                                                                                                         tGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGc                                             <  2:3811126/150‑1 (MQ=255)
                                                                                                           cGGCACTGGCTCCCATGGCGAAAGGGAAGGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGccc                                           <  2:905233/150‑1 (MQ=255)
                                                                                                            ggCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTAGGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGcccc                                          <  2:3862456/150‑1 (MQ=255)
                                                                                                             gCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCACACCCCAAACCACGCCCACCGGGGGAGGGGGCCCGTCGGCGCAGTATGTAAAACCACACACGGAAAACACACCCCCCAAGTTTTTTCTCCCACGccccc                                         >  2:3914509/1‑149 (MQ=255)
                                                                                                               aCTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTTCGCCCCCCAACCCAAAACCTATACCACAGGGGGGGGGGGCCCTTCGGCGCAGTATGTACACCCACACTCGTAAACGACACCCCAAACTTTTTTTATCCCTCCCCCCAcc                                       >  2:2045299/1‑150 (MQ=255)
                                                                                                               aCTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCGCACCCGAACCCCCACCCCCCGGGGGGGGGGACCCGTCGGCGCGGGATGTTCTCCCGCACACGGCAACCACCCCCCCCACGTTTTTTGACCCACGCCCCGcc                                       >  2:213559/1‑150 (MQ=255)
                                                                                                                  ggCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCggg                                    <  2:1101018/150‑1 (MQ=255)
                                                                                                                  ggCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCggg                                    <  1:1630004/150‑1 (MQ=255)
                                                                                                                  ggCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCggg                                    <  1:3792747/150‑1 (MQ=255)
                                                                                                                  ggCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCggg                                    <  1:1062814/150‑1 (MQ=255)
                                                                                                                  ggCTCCCATGGCGAAAGGGAATGCGGCATAGGCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCggg                                    <  2:1591329/150‑1 (MQ=255)
                                                                                                                   gcgcccATGGCGAAAGGGAATGCGGCATAGTCGTGTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACCCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCgggg                                   <  2:2951262/147‑1 (MQ=255)
                                                                                                                   gCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCCCACCCCAACCCACGCCCACAGGGGGAGGGGCCCCGTGGGGGCGGGTTGTGCTCCCGCACACAGAAAACACCCCCCCCTCATTTTTTTTACCCCGCCCCCCCcggg                                   >  1:3732964/1‑150 (MQ=255)
                                                                                                                    cTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCCCACCCGAACCCACACCCACAGGGGGTGGGGACCCGTGGGGGCTGTATGTAATCCCACACACAGAAACGACACCCCACTCTTTTTTTATCCCACGCCCCGCCggggg                                  >  2:2225769/1‑149 (MQ=255)
                                                                                                                       ccATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCAAACCCACTCCCACAGGGAGAGGGGTCCGTTCGGGGCTGTATGTTATACCGCACTCTGGAAACGCACCCCACACTTTTTTTATCCCGCGCCCCACCgggggggg                               >  2:2420411/1‑146 (MQ=255)
                                                                                                                       ccATGGCGAAAGGGAATGCGGCATAGTCGTGTGCGCCCCCACACCCGAAACCACACCCTCAGGGGGAGGGGACCCTTCGGCGCGGTATTTACACCCACACACAGAAACGACAACGCACACATTTTTTATACCACGCCCCACCCGGGagag                               >  2:1016910/1‑150 (MQ=255)
                                                                                                                        cATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAgg                              <  1:1498332/150‑1 (MQ=255)
                                                                                                                         aTGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAggg                             <  2:34773/150‑1 (MQ=255)
                                                                                                                         aTGGCGAAAGGGAATGCGGCATAGTCGTGTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGGTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTGTCTCCCTCGCCCCTCCGGGGAGAggg                             <  2:2721424/150‑1 (MQ=255)
                                                                                                                          tGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCGCACCCGAAACCACTCCCACAGGGCGAGGGGACCCGTCGGCGCGGTATGTACACCCGCACTCTGAAACCACACCCCACTCTTTTTTTCTCCCGCGCCCCACCCGGGCGAGGGc                            >  2:4054117/1‑150 (MQ=255)
                                                                                                                             cGAAAGGGAATGCGGCATAGTCGTGTGCGCCCCCACAACCAAAACCTCAACCACAGGGAGAGGGGACCCTTCGGGCCTGTTTGTAAGCCCGCAAGATGAAACGACACCCCACAATTTTTTTATCCCACCCCCCCACGGGGGGTGGGGCgg                         >  2:957495/1‑150 (MQ=255)
                                                                                                                               aaaGGGAATGCGGCATAGTCGTTTTCGCCCCCACACCCAAACCCAAACCCACAGGGGGTGGGGACCCGTCGGCGCGGTATGTACACCCGCCCACGGAAAACACACCGCACACGTTTTTTATCCCACCCCCCGCCGGGGCGAGGGCCgggg                       >  2:1347375/1‑150 (MQ=255)
                                                                                                                                   ggAATGCGGCATAGTCGTTTGCGCCCCCTCACCCGAACCCTCTCCCTCCGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAg                   <  2:1477077/150‑1 (MQ=255)
                                                                                                                                    gAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCCCTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAgg                  <  2:1130010/150‑1 (MQ=255)
                                                                                                                                          ggCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCCCTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGaaaa            <  2:3129770/150‑1 (MQ=255)
                                                                                                                                          ggCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGaaaa            <  2:2344421/150‑1 (MQ=255)
                                                                                                                                          ggCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCGCCGGGGAGTGGGCCGGGGTGAGGGGaaaa            <  1:905547/150‑1 (MQ=255)
                                                                                                                                                    tttGCGCCCCCTCAACCAAACCCACAACCACAGGGAGAGGGGACCGGTCCGCGCAGTATTTAAACCCTCACACTGAAACCACACCCCACACATTTTTTCTCCCACGCCCCACCGGGGCGAGGGCCGGGGTGAGGGGAAAAGGCCCCAcag  >  2:3161437/1‑148 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTG  >  NC_000913/3780873‑3781170

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: