Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,222,684 G→C 7.5% intergenic (+140/+233) yohP → / ← dusC uncharacterized protein/tRNA‑dihydrouridine synthase C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,222,6840GC 7.5% 286.1 / 3.8 109intergenic (+140/+233)yohP/dusCuncharacterized protein/tRNA‑dihydrouridine synthase C
Reads supporting (aligned to +/- strand):  ref base G (23/77);  new base C (2/7);  total (25/84)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 120 total aligned reads displayed.

GATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGC‑CAA‑T‑CGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCC  >  CP009273/2222535‑2222833
                                                                                                                                                         |                                                                                                                                                     
gatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGC‑CAA‑T‑CGCCGCCGCCCCCTg                                                                                                                                                       <  1:660884/150‑1 (MQ=255)
gatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGC‑CAA‑T‑CGCCGCCGCCCCCTg                                                                                                                                                       <  1:660917/150‑1 (MQ=255)
   tatttAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGC‑CAA‑T‑CGCCGCCGCCCCCTGCAc                                                                                                                                                    <  2:1919390/148‑1 (MQ=255)
    ttttAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGC‑CAA‑T‑CGCCGCCGCCCCCTGCAcc                                                                                                                                                   <  2:2395765/150‑1 (MQ=255)
          aattaattAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GGGGGGGGGGGG‑CGG‑G‑CCCCCCCGccccccccccccccccc                                                                                                                                             >  1:2951259/1‑147 (MQ=255)
           attaattaATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGTGGGCGGC‑CAA‑A‑CGCCCCCCCCCCCAGCACCCCCCGGc                                                                                                                                            >  1:2654731/1‑150 (MQ=38)
            ttaattaaTGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GGGGGGGGGGGG‑GCA‑C‑CCCCCCCGccccccccccccccccccc                                                                                                                                           >  2:3806317/1‑145 (MQ=255)
            ttaattaaTGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GGGGGGGGGGGG‑CCA‑G‑AGccccccccccccccccccccccccc                                                                                                                                           >  2:1389186/1‑145 (MQ=255)
            ttaattaaTGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GGGGGGGGGGGCGCAA‑A‑CAcccccccccccccccccccccccc                                                                                                                                            >  2:1389228/1‑146 (MQ=255)
              aattaatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGGGGGGGGGGGC‑CAA‑A‑AGCC‑CCGCCCCCCCCAcccccccggcca                                                                                                                                         >  2:3733410/1‑143 (MQ=255)
               attaatGTCATCAGGTCCGAAAATAACGATAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGC‑CAA‑T‑CGCCGCCGCCCCCTGCACCCCCGGGCTCTg                                                                                                                                        <  2:375172/150‑1 (MQ=255)
               attaatGTCATCAGGTCCGAAAATAACGAGAATTTTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGC‑CAA‑T‑CGCCGCCGCCCCCTGCACCCCCGGGCTCTg                                                                                                                                        <  2:4147390/150‑1 (MQ=255)
                ttattGTCATCTGGTCCGAAAATAACGAGAATATTTCTGTCTCTCTTCCTGTTGCGCGCCTGTATTGTGCTTTGCTTCATATAATCACTGGGGCATGGTGGGC‑GCTGGGGGCGGC‑CAT‑T‑CGCGGCCGCCCCCTGCACCCCCGGGCTCTgg                                                                                                                                       <  2:1626786/150‑1 (MQ=25)
                      tcatcaGGTCCGAAAATAACGAGAATAGTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GGCCGGGGCGGG‑CCA‑C‑AACCCCCCCCCCCCGCCCCCCCCCGCTCCGGAGGACa                                                                                                                                 >  2:794848/1‑150 (MQ=255)
                          ctggTCCTAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCCGGGGGCGGC‑CAA‑T‑CGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAAt                                                                                                                             <  1:2741808/148‑1 (MQ=255)
                          caGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGC‑CAA‑T‑CGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAAt                                                                                                                             <  2:1091121/150‑1 (MQ=255)
                            ggTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGC‑CAA‑T‑CGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATcg                                                                                                                           <  2:3888458/150‑1 (MQ=255)
                             gTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGCGGGC‑CCA‑‑‑AGCCCCCCCCCCCCGCCCCCCCCGGCACAGGCCAaaaaaaaaaac                                                                                                                         >  2:1123023/1‑145 (MQ=25)
                                             aaTATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGATTCATATAATCACTGGCGCAAGGAGCGC‑GCTGGGGGCGGC‑CAA‑T‑CGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCt                                                                                                          <  2:108807/150‑1 (MQ=255)
                                              ctatTCCAGTCTCTCTTCGTGGGGCGCGCCGGTCCTGGCCTTTGCTTCATATAACCCCGGGCGCCCGGCGCGC‑GCGGGGGGCGCC‑CAA‑C‑CGCCGCCGCCCCCGCCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTc                                                                                                         <  1:3857315/149‑1 (MQ=9)
                                                 tttCAGACTCTCATCCTGTGGCGCTCGTGTCTTGTGCTTTGCTACATATAATCACGGGCGCACGGAGCGC‑GCAGGGGGCGGC‑CAA‑A‑CGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGc                                                                                                      <  2:3689980/150‑1 (MQ=17)
                                                 ttaaaGTCTCTCATCGTGTTGCGCTCCTGTCATGTGATTTGCTTAATATAAACCCGGGCGCATGGCGCGC‑GCCGGGGGCGGC‑CAA‑T‑GGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGc                                                                                                      <  1:379305/146‑1 (MQ=17)
                                                     aGTCTCTCATCGTGTTGCGCTCCTGTCTTGTGATTTGATTAATATAAACACGGGCGCATGGAGCGC‑GCGGGGGGCGGC‑AAA‑T‑CGCCGCCGCCCCCTCCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTAt                                                                                                  <  1:3398474/150‑1 (MQ=17)
                                                       tctctcATCCTGTTGGGCTCCTGTCTGGTGATTTGCTTAATATAATCACTGGCGCAAGGAGCGC‑GCTGGGGGCGGC‑CAA‑T‑CGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATcc                                                                                                <  1:1471656/150‑1 (MQ=17)
                                                          ctcATCGTGTGGCGCGCCTGTCATGTGCATTGCTTCATATAAGCACGGGCGCACGGAGCGC‑GCAGGGGGCGGC‑CAA‑T‑CGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCtt                                                                                             <  2:4419686/150‑1 (MQ=17)
                                                           tcATCGTGTGGCGCTCGTGTCTTGTGATTTGCTTCATATAACCCCGGGCGCACGGCGCGC‑GC‑GGGGGGGGC‑CAA‑TGCCCCCCCCCCCCCTCCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTa                                                                                            <  1:515654/150‑1 (MQ=9)
                                                           taaTCTTGTTGGGCTCGTGTCTTTTGATTTGTTTCATATAATCACTGGGGCACGGAGCGC‑GCTGGGGGCGGC‑CAA‑T‑CGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTa                                                                                            <  2:3180914/148‑1 (MQ=14)
                                                             agccggTTGCGCTCCTGTCATGTGCATTGCTTCATATATTCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGC‑CAC‑T‑CGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACg                                                                                          <  2:2601983/145‑1 (MQ=255)
                                                             aTCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGC‑CAA‑T‑CGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACg                                                                                          <  2:2601435/150‑1 (MQ=255)
                                                              tCCTGTTGCGCTCCTGTCTTGTGAATTGCTTAATATAATCCCGGGCGCATGGTGCGC‑GCTGGGGGCGGC‑AAA‑T‑CGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACgg                                                                                         <  2:2017361/150‑1 (MQ=17)
                                                                 tgtgGCGCTCCTGTCTTGTGATTTGCTTAATATAAACACGGGCGCAGGGCGGGC‑GCTGGGGGGGGC‑CAA‑T‑CGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTa                                                                                      <  1:1518530/150‑1 (MQ=14)
                                                                     gcgcTCCTTTCTTGTGCATTGCTTAATATAATCACTGGCGAAAGGAGCGC‑GCTGGGGGCGGC‑CAA‑T‑CGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCgg                                                                                  <  2:2610724/150‑1 (MQ=17)
                                                                      cgcTCCGGTCCTGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCTGGGGGCGGC‑CAA‑T‑CGCCGCCGCCCCCGGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCggg                                                                                 <  2:4489809/150‑1 (MQ=255)
                                                                           ttttCATGTGATTTGCTTAATATAATCACTGCCGCCAGGAGCGC‑GCTGGGGGCGGC‑CAA‑G‑CGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCggg                                                                            <  1:2206086/147‑1 (MQ=14)
                                                                           gtGTCTTGTGCTTTGGTAAATATAATCCGGGGCGCAGGGGGGGG‑GCGGGGGGGGCC‑AAA‑T‑CGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCggg                                                                            <  1:2206555/149‑1 (MQ=11)
                                                                             gTCATGTGCATTGCTTCATATAATCCCTGGCGCAAGGAGCGC‑GCAGGGGGCGGC‑CAA‑T‑CGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGgcg                                                                          <  2:4105333/150‑1 (MQ=255)
                                                                              tCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGC‑CAA‑T‑CGCCGCCGCCCCCGGCACCCCCGGGCACTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGgcgc                                                                         <  2:2277790/150‑1 (MQ=255)
                                                                                aggtgCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGC‑CAA‑T‑CGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGa                                                                       <  2:3745421/148‑1 (MQ=255)
                                                                                     cATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGC‑CAA‑T‑CGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTaa                                                                  <  2:1071707/150‑1 (MQ=255)
                                                                                           ttCATATAATCACTGGCGCATGGAGCGC‑GCAGGGGGCGGC‑CAA‑T‑CGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATccc                                                            <  2:179311/150‑1 (MQ=255)
                                                                                           ttCATATAATCACTGGCGCACGGAGCGC‑GCAGGGGGCGGC‑CAA‑T‑CGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATccc                                                            <  2:179000/150‑1 (MQ=255)
                                                                                                      aCTGGCGCAAGGAGCGC‑GCAGGGGGCGGC‑CAA‑T‑CGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAcgcg                                                 <  2:4590367/150‑1 (MQ=32)
                                                                                                          gcgcAAGGAGCGC‑GCAGGGGGCGGC‑CAA‑T‑CGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCt                                             <  2:1846335/150‑1 (MQ=32)
                                                                                                            tcatgttgcgc‑gcTGGGGG‑GGC‑CAATT‑CGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCTAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCa                                           <  1:1014928/143‑1 (MQ=11)
                                                                                                             cAAGGAGCGC‑GCTGGGGGCGGC‑CAA‑T‑CGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAg                                          <  2:1472631/150‑1 (MQ=17)
                                                                                                             cAAGGAGCGC‑GCAGGGGGCGGC‑CAA‑T‑CGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAg                                          <  2:4709866/150‑1 (MQ=32)
                                                                                                             aatgttgcgc‑gcGGGGGGCGCC‑CAA‑G‑CGCGGCCGCCCCCGCCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAg                                          <  1:1019237/144‑1 (MQ=11)
                                                                                                              aaGGAGCGC‑GCAGGGGGCGGC‑CAA‑T‑CGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGc                                         <  2:3313908/150‑1 (MQ=32)
                                                                                                               aGGAGCGC‑GCAGGGGGCGGC‑CAA‑T‑CGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGcc                                        <  2:2667430/150‑1 (MQ=32)
                                                                                                                 gtgcgc‑gcTGGGGGCGGC‑AAA‑T‑CGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCca                                      <  1:1130398/148‑1 (MQ=12)
                                                                                                                 gAGCGC‑GCAGGGGGCGCC‑CAA‑T‑CGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCca                                      <  1:1907931/150‑1 (MQ=17)
                                                                                                                        gcAGGGGGCGGC‑CAA‑T‑CGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCa                                <  1:4295242/150‑1 (MQ=32)
                                                                                                                        gcAGGGGGCGGC‑CAA‑T‑CGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCa                                <  1:4295209/150‑1 (MQ=32)
                                                                                                                         ctgggggCGGC‑CAC‑T‑CGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAt                               <  1:3427214/148‑1 (MQ=17)
                                                                                                                         ctgggggCGGC‑CAA‑G‑CGCCGCCGCCCCGGGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAt                               <  1:3426968/148‑1 (MQ=12)
                                                                                                                         ccgggggCGGC‑CAA‑T‑CGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAt                               <  2:3866484/148‑1 (MQ=17)
                                                                                                                         cAGGGGGCGGC‑CAA‑T‑CGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCag                               >  2:3067185/1‑149 (MQ=17)
                                                                                                                          aGGGGGCGGC‑CAA‑T‑CGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCGTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATACCAGTAAAACGCGCCCTCAGCCCCCAAcccag                              >  1:512218/1‑147 (MQ=11)
                                                                                                                           gggggCGGC‑CAA‑T‑CGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAtgt                             >  2:668237/1‑150 (MQ=32)
                                                                                                                           gggggCGGC‑CAA‑T‑CGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAtgt                             <  2:2599353/150‑1 (MQ=32)
                                                                                                                           gggggCGGC‑CAA‑T‑CGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAtgt                             <  2:2423577/150‑1 (MQ=32)
                                                                                                                           gggggCGGC‑AAT‑T‑CGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAtgt                             <  2:2423629/150‑1 (MQ=17)
                                                                                                                            ggggCGGC‑CAA‑T‑CGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAtgtg                            >  1:1550729/1‑150 (MQ=32)
                                                                                                                            ggggCGGC‑CAA‑T‑CGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAtgtg                            >  1:1551498/1‑150 (MQ=32)
                                                                                                                            ggggCGGC‑CAA‑T‑CGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCCGCTTATAACTTACGGCTACCCGGTCGGGCGCGAGGTAACCTCCCTGTAAAAAGCCCCCACCGCCCCCAACCAagtg                            >  1:3322607/1‑150 (MQ=11)
                                                                                                                            ggggCGGC‑CAA‑T‑CGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCGCGGCTTATCCCTTACGGCTACCCGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCACAGCCCCCATCCAtgtg                            >  1:202921/1‑150 (MQ=255)
                                                                                                                             gggCGGC‑CAA‑T‑CGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTgg                           <  1:2614906/150‑1 (MQ=32)
                                                                                                                                cggc‑CAA‑T‑CGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCt                        <  1:1495467/150‑1 (MQ=32)
                                                                                                                                cggc‑CAA‑T‑CGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCGCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCAGTAAAACGCGCCCACAGCCCACAACCATGTGGGCt                        >  1:2121384/1‑150 (MQ=255)
                                                                                                                                 ggc‑CAA‑T‑CGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTg                       <  2:3499682/150‑1 (MQ=32)
                                                                                                                                 ggc‑CAA‑T‑CGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTg                       >  2:77948/1‑150 (MQ=32)
                                                                                                                                  gggcAA‑G‑CGCCGCGGCGCCCTGATCCCCCGCGCGCTGGCGAGAAAAATATCCGCGGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTg                       <  2:4645055/147‑1 (MQ=11)
                                                                                                                                   c‑cAA‑G‑CGCCGCCCCCCCCTCCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGcc                     <  1:915493/150‑1 (MQ=12)
                                                                                                                                     cAA‑T‑CGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGccc                    >  1:3078280/1‑150 (MQ=32)
                                                                                                                                     cAA‑T‑CGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGccc                    >  1:3078285/1‑150 (MQ=32)
                                                                                                                                     cAA‑T‑CGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGccc                    >  1:3078473/1‑150 (MQ=32)
                                                                                                                                     cAA‑T‑CGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGccc                    <  2:1841716/150‑1 (MQ=32)
                                                                                                                                      aa‑T‑CGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGcccc                   >  2:603955/1‑150 (MQ=32)
                                                                                                                                            gccgccgccCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGcc               <  2:889267/150‑1 (MQ=32)
                                                                                                                                             ccgccgccCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCt              <  2:3848825/150‑1 (MQ=32)
                                                                                                                                               gccgccCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTc            <  1:619163/150‑1 (MQ=32)
                                                                                                                                               gccgccCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCGTATACCGTACGGCTACCGGGTCGGGCGCGAGGTGACATCCCTGTAAAACGCGCCCACAGCCCACCACCATGTGGGGTGCCCCGGCCGTc            >  1:3294327/1‑150 (MQ=11)
                                                                                                                                               gcccccCCGTCCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTc            <  1:3706447/150‑1 (MQ=12)
                                                                                                                                                cgccccccGGCACCCCCGGGCTCTGTCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCa           <  1:926744/147‑1 (MQ=12)
                                                                                                                                                 cgccCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACCTCCATGTGGGCTGCCCCGGCCTTCAg          >  2:1671430/1‑150 (MQ=17)
                                                                                                                                                 cgccCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCGTCAg          >  1:4208786/1‑150 (MQ=17)
                                                                                                                                                   ccccccGCACCCCCGGGCTCTGGCGAACAAACTCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAggg        <  1:28208/150‑1 (MQ=17)
                                                                                                                                                   ccccccGCACCCCCGGGCTCTGGCGAACAAACTCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAggg        <  1:3609622/150‑1 (MQ=17)
                                                                                                                                                   ccccccGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAggg        <  1:27911/150‑1 (MQ=17)
                                                                                                                                                   cccccTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCGTATACCGTAAGGCGACCCGGTCGGGCGCGCGGGAACAAACCAGTAAAAACCCCCCACAGCCCACCACCAAGTGGGGTGGCCCCGCCGTCAcgg        >  1:2794695/1‑150 (MQ=2)
                                                                                                                                                   cccccGGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAggg        <  1:1216451/150‑1 (MQ=17)
                                                                                                                                                   cccccGCCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAggg        <  1:1216206/150‑1 (MQ=17)
                                                                                                                                                    ccccTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGa       <  1:1796728/150‑1 (MQ=32)
                                                                                                                                                    ccccTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGa       <  1:2355273/150‑1 (MQ=32)
                                                                                                                                                       cTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACg    <  2:115825/150‑1 (MQ=32)
                                                                                                                                                       cTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACg    <  1:1912644/150‑1 (MQ=32)
                                                                                                                                                         gCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACCTCCATGTGGGCTGCCCCGGCCTTCAGGGAACGcc  <  2:3567498/150‑1 (MQ=17)
                                                                                                                                                         gCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGcc  >  1:3039755/1‑150 (MQ=32)
                                                                                                                                                         gCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGcc  >  1:3039707/1‑150 (MQ=32)
                                                                                                                                                         gCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGcc  <  1:3821242/150‑1 (MQ=32)
                                                                                                                                                         |                                                                                                                                                     
GATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGC‑CAA‑T‑CGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCC  >  CP009273/2222535‑2222833

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: ‑