Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,222,684 G→C 6.0% 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 6.0% 397.2 / 3.5 146intergenic (+140/+233)yohP/dusCuncharacterized protein/tRNA‑dihydrouridine synthase C
Reads supporting (aligned to +/- strand):  ref base G (28/108);  new base C (6/4);  total (34/112)
Fisher's exact test for biased strand distribution p-value = 1.10e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 147 total aligned reads displayed.

GATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATT‑GCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACG  >  CP009273/2222535‑2222831
                                                                                                                                                       |                                                                                                                                                    
gatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATT‑GCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTg                                                                                                                                                      <  1:3306788/150‑1 (MQ=255)
   atTTTAAAATTAATTAATGTCATCAGGTCCGAACATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATT‑GCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCAc                                                                                                                                                   <  2:1356357/150‑1 (MQ=255)
     tttAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATT‑GCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCAccc                                                                                                                                                 <  1:2711285/150‑1 (MQ=255)
       tAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATT‑GCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCAccccc                                                                                                                                               <  2:4742067/150‑1 (MQ=255)
        aaaaTTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATT‑GCTTCATATAATCACTGGCGCAAGGAGCGC‑GGGGGGGGGGGGCCCGCGCCCCCGCCCCCCGCccccccc                                                                                                                                              >  1:5253482/1‑149 (MQ=255)
        aaaaTTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATT‑GCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCg                                                                                                                                              <  1:4333695/150‑1 (MQ=255)
           attaattaATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATT‑GCTTCATATAATCACTGGCGCAAGGAGCGC‑GGAGGGGGCGGGCCCACGCCCCCGCCCCCCGCACCCCCCGGc                                                                                                                                           >  2:1229431/1‑150 (MQ=255)
            ttaattaaTGTCATCTGTTCCGAAATTAACGTTAATTTTTCAGTCTCTCATCCTGTTGCGCTCCTTACATTTGCATT‑GCTTCATATAATCAGTGGCGCACGGAGCGC‑GCTGGGGGCGGCCAATCGCGGCCGCCCCCTGCACCCCCGGGct                                                                                                                                          <  2:3785755/150‑1 (MQ=21)
            ttaattaaTGTCATCAGTTCCTAAAATAACGAGAATATTTCTTTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATT‑GCTTAATATAATAACTGGCGAAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGct                                                                                                                                          <  2:4327308/150‑1 (MQ=25)
            ttaattaaTGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATT‑GCTTCATATAATCACTGGCGCAAGGAGCGCGGCGGGGGGGGGCCCA‑AGCCGCCGCCCCCCGCAcccccccccc                                                                                                                                          >  2:4088952/1‑145 (MQ=255)
            ttaattaaTGTAATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATT‑GCTTCATATAATCACTGGCGCACGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGct                                                                                                                                          <  2:1705543/150‑1 (MQ=255)
             gaattacgGTCTTCTGGTCCGAAATTAACGAGAATAGTTCTGTCTCTCTTCCTGGTGCGCTCGTGTCATGTGGATT‑GCGTCAGTTATTGAGTGGCGCATGGAGCGC‑GCTGGGGGCGGGCAAGCGCCGCGGCCCCCTGCACCCCCGGGctc                                                                                                                                         <  2:4414073/149‑1 (MQ=16)
              aattaatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATT‑GCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGctct                                                                                                                                        <  2:5313039/150‑1 (MQ=255)
               attaatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATT‑GCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTg                                                                                                                                       <  2:953470/150‑1 (MQ=255)
                  attGTCATCTGGTCCTAAAATAACTATAATATTTCTGTCTCTAATCCTGTTGCGCTCCTGTCATGTGCTTT‑GCTTCATATAATCACTGGCGCAAGGAGCGC‑GCGGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCg                                                                                                                                    <  1:518659/148‑1 (MQ=21)
                  aatGACAGCCGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATT‑GCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCg                                                                                                                                    <  2:4877159/150‑1 (MQ=255)
                   atGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATT‑GCTTCATATAATCACTGGCGCAAGGAGCGCGGCGGGGGGGGGCCAA‑CGCCGCCCCCCCCCCCCCCCCCCGCCGCAGGCG‑a                                                                                                                                  >  1:4445830/1‑150 (MQ=255)
                    tGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATT‑GCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑aa                                                                                                                                 <  2:843628/150‑1 (MQ=255)
                          ctggTCCTAAAATAACGAGAATTTTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATT‑GCTTCATATAATCACGGGCGCAAGGAGCGC‑GCCGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAAt                                                                                                                           <  1:3544988/148‑1 (MQ=255)
                          caGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATT‑GCTTCATATAATCACTGGCGCACGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAAt                                                                                                                           <  2:2565569/150‑1 (MQ=255)
                           aGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATT‑GCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATc                                                                                                                          <  2:433666/150‑1 (MQ=255)
                           aGGTCCGAAAATAACGAGAATATTTCAGTCTCTCAGCCTGTTGCGCTCCTGTCATGTGCATT‑GCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATc                                                                                                                          <  2:457234/150‑1 (MQ=255)
                               ccGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATT‑GCTTCATATAATCACTGGCGCAAGGAGCGC‑GGGGGGGGGGGGCCAGAACCCCCCCCCCCCCCACCCCCCCCCCCAGGCG‑GACaaaaaaccac                                                                                                                      >  1:3011294/1‑144 (MQ=255)
                               ccGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATT‑GCTTCATATAATCACTGGCGCAAGGAGCGC‑GGAGGGGGGGGGCAAAAACCCCCCCCCCCCCCACCCCCCCGCCCTGGCGCAACaaaaaaccc                                                                                                                       >  1:328635/1‑145 (MQ=255)
                               ccGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATT‑GCTTCATATAATCACTGGCGCAAGGAGCGC‑GGAGGGGGGGGCCAAAAGCCGCCCCCCCCCGCCCCCCCCGGCGCAGGCG‑CACaaaaaacccc                                                                                                                      >  2:3937937/1‑144 (MQ=255)
                               ccGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATT‑GCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATcgccg                                                                                                                      <  2:2423674/150‑1 (MQ=255)
                                      tAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATT‑GCTTCATATAATCACTGGCGCATGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCACCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGcgcg                                                                                                               <  1:575773/150‑1 (MQ=255)
                                         ctataaTATTTCAGTCTCTCATCCTGTTGCGATCCTGTCATGTGCATT‑GATTCATATAAACACTGGCCAACGGAGCGC‑GCAGGGGGCGGCAAAACGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCT‑AACAAAATCGCCGCTGCGCGGTg                                                                                                            <  2:1895504/146‑1 (MQ=18)
                                         cgtcattatTTCATTCTATCATCCTGTTGCGCTCCTGTCTTGTTATTT‑GATTAATATAATCCTGGGCGAATGGCGGGC‑GCCGGGGGCGGCCAATCGCCGCCGCCCCCGGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTg                                                                                                            <  1:2920768/144‑1 (MQ=17)
                                         cGAGAATATTTCTGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATT‑GCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTg                                                                                                            <  2:4267393/150‑1 (MQ=255)
                                          gagaATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCTTGTGCATT‑GCTTCATATAATCACGGGCGCACGGAGCGC‑GCTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGc                                                                                                           <  2:4067941/150‑1 (MQ=255)
                                            gaATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATT‑GCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGccc                                                                                                         <  2:1088738/150‑1 (MQ=255)
                                              ttttttCAGTCTCTCTTCCTGTTGCGCTCCTGTCTTGTGAATT‑TCTTAATATAATCACTGGCGCAAGGAGGGC‑GCTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTc                                                                                                       <  1:3460546/147‑1 (MQ=18)
                                              ttatTTCTGTATCCCATCCTGTGGCGCTCCTGTCTTGT‑CATTGGCTACATATAACCCGTGGCGCACGGTGCGC‑GCTGGGGGCGGCAAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTc                                                                                                       <  1:698175/149‑1 (MQ=14)
                                               tatTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATT‑GCTTCATATAATCCCGGGCGCACGGAGCGC‑GCCGGGGGCGGCCACCCGCCGCCGCCCCGTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCg                                                                                                      <  2:2028840/150‑1 (MQ=18)
                                                   tctgTCTCTCACCCTGTGGCGCCCCTTTATTGTGCATT‑GCTTCATAAAACCCCGGGCGCAGGGCGCGC‑GCCGGGGGCGGCCAAGCGCCCCCGCCCCCGCCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCtt                                                                                                  <  1:3201430/147‑1 (MQ=17)
                                                      gTCTCTCATCCTGTTGCGCTCCTCTCATGTGCATT‑GCTTCATATAATCACTGGCGCAAGGAGCGC‑GCCGGGGGCGGCCCAAAGGCCCCCCCCCCCCCACCCCCCGGCACAGGCG‑AACAAAAAAACCCCATCGCCGTGCCCCCCGCTTAta                                                                                               >  2:1247384/1‑149 (MQ=16)
                                                          ctcATCCTGTTGCGCTCCTGTCATGTGCATT‑GCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCtt                                                                                           <  2:2606545/150‑1 (MQ=255)
                                                           tcATCCTGTTGCGCTCGTGTATGGTTCTTT‑GCTAAATATAACCCCGGGCGCATGGAGGGC‑GCTGGGGGCGCCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTa                                                                                          <  1:1051230/150‑1 (MQ=14)
                                                           tcATCCTGTGGCGCTCCTGTCATGTGCATT‑GCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTa                                                                                          <  2:2774783/150‑1 (MQ=255)
                                                            cATCGTGTTGCGCTCCTGTCTGGTGCATG‑GCTTCATATCAGCACTGGCGCAAGGAGCGC‑GCTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCGCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTAc                                                                                         <  2:3048182/150‑1 (MQ=17)
                                                             aTCCTGTTGCGCTCCTGTCATGTGCATT‑GCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATAGCCGCCGCCCCCAGTACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATACCGTACg                                                                                        >  2:2019988/1‑150 (MQ=18)
                                                             aTCCTGTTGCGCTCCTGTCATGTGCATT‑GCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAAAAGCCGCCGCCCCCAGCACCCCCCGGCGCAGGCG‑AACAAAAACCCCCCCGCGCGGGGCCCACGGCTTATACCATAcc                                                                                        >  2:3566304/1‑149 (MQ=18)
                                                             aTCCTGTTGCGCTCCTGTCATGTGCATT‑GCTTCATATAATCACGGGCGCACGGAGCGC‑GCTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACg                                                                                        <  1:3048286/150‑1 (MQ=255)
                                                                   ttGCGCTCCTTTCATGTGCATT‑GCTTCATATAATCCCGGGCGCATGGAGCGC‑GCTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTAcc                                                                                  <  1:153133/150‑1 (MQ=17)
                                                                   ttGCGCTCCTGTCATGTGCATT‑GCTTCATATAAACACTGGCGCACGGCGCGC‑GCCGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTAcc                                                                                  <  2:2471387/150‑1 (MQ=255)
                                                                   ttGCGCTCCGGTCTTGTGCTTT‑GCTTAATATAATCACTGGCGCATGGTGCGC‑GCGGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTAcc                                                                                  <  1:153136/150‑1 (MQ=17)
                                                                   ttGCGATCCTGTCTTGGTCATT‑GATTAATATAACCCCGGGCGCAGGGCGGGC‑GCGGGGGGGGCCCAACCGCCGCCCCCCCCTCCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTAcc                                                                                  <  1:1404114/150‑1 (MQ=14)
                                                                       gcTCCTTTCTTGTGCATT‑GCTTCATATAATCACTGGCGCATGGAGCGC‑GCTGGGGGCGGCAAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGt                                                                              <  2:1884196/150‑1 (MQ=17)
                                                                        cTCCTGTCTGGTGCATT‑GCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTc                                                                             <  2:745683/150‑1 (MQ=255)
                                                                         tCCTGACTTGTGCATT‑GCTACATATAATCACTGGCGCAAGGAGCGC‑GCTGGGGGGGGCAAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCg                                                                            <  1:1677385/150‑1 (MQ=14)
                                                                                    gCATT‑GCTTCATATAAGCACTGGCGCAAGGAGCGC‑GCCGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTa                                                                 <  2:3157809/150‑1 (MQ=255)
                                                                                        t‑GCTTCATATAATCACTGGCGCATGGAGCGC‑GCTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACAt                                                             <  1:1174151/150‑1 (MQ=255)
                                                                                           cTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATcc                                                           <  2:5088796/150‑1 (MQ=40)
                                                                                             tCATATAAGCCCGGGCGCACGGCGCGC‑GCGGGGGGCGGCCAAGCGCCGCCGCCCCCTGCACCCCCGGGCTCGGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCt                                                         <  1:2123854/150‑1 (MQ=11)
                                                                                                 ataATCACTGGCGCAAGGAGCGC‑GGAGGGGGGGGGCCAACCCCGCCGCCCCCAGCACCCCCCGGCGCTGGCG‑AACAAAAACGCCCCAGCGCCGGGGCCGCGGCGTATACCATAAGGCGTCCCGGGCGGGCGCGAGGGGACAACCCAGTaa                                                     >  2:3090240/1‑150 (MQ=255)
                                                                                                 ataATCACTGGCGCAAGGAGCGC‑GCACGGGGCGGCCCAACGCCCCCCCCCCCACCCCCCCCCGGCACAGGCG‑AACCAAAACGCCCCCGCCCGGGGCCCCCGGCGTATCCCAAACCGCAACCCGGGCGGGGGCCCGGTAACAAACCCGTaa                                                     >  2:3110641/1‑150 (MQ=255)
                                                                                                        gtGGCGCAGGGTGCGC‑GCGGGGGGCGGCAAATCGCCGCCGCCCCCTCCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACgcgc                                              <  1:4061340/149‑1 (MQ=11)
                                                                                                            cgcAAGGAGCGC‑ACAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTc                                          >  2:260493/1‑150 (MQ=17)
                                                                                                             gcAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCa                                         <  1:5256545/150‑1 (MQ=32)
                                                                                                              cAAGGAGCGC‑GCAGGGGGCGGCCAAGCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAACACGCGCCCTCAg                                        <  2:2599205/150‑1 (MQ=17)
                                                                                                                  gAGCGC‑GCCGGGGGCGGCCAAGCGCCGCCGCCCCCGGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCca                                    <  1:205801/150‑1 (MQ=12)
                                                                                                                  gAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCca                                    <  2:3328529/150‑1 (MQ=32)
                                                                                                                     cgc‑ggAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACAt                                 <  2:4647913/150‑1 (MQ=17)
                                                                                                                     cgc‑gcAGGGGGCGCCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACAt                                 <  2:1362985/150‑1 (MQ=17)
                                                                                                                       c‑gcAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATcc                               <  1:2230767/150‑1 (MQ=32)
                                                                                                                         gcAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCa                              <  1:265697/150‑1 (MQ=32)
                                                                                                                          ctgggggCGGCCAAGCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAt                             <  1:1408953/148‑1 (MQ=17)
                                                                                                                          cAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAt                             <  2:3808795/150‑1 (MQ=32)
                                                                                                                          cAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAt                             <  2:3396297/150‑1 (MQ=32)
                                                                                                                             ggggCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAtgtg                          <  1:3409725/150‑1 (MQ=32)
                                                                                                                             ggggCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAtgtg                          >  1:5327231/1‑150 (MQ=32)
                                                                                                                             ggggCGGCAAATCGCCGCCGCCCCCTGCACCCCCGGTCTCTGGCG‑AACAAAATCGCCGCTGAGCGGTGCCCTAGGATTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAtgtg                          <  1:573490/150‑1 (MQ=11)
                                                                                                                              gggCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTgg                         <  1:3758801/150‑1 (MQ=32)
                                                                                                                              gggCGGCCAATCGCCGCCGCCCCCGGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTgg                         <  2:2909106/150‑1 (MQ=17)
                                                                                                                                 cggcCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCt                      <  2:4881605/150‑1 (MQ=32)
                                                                                                                                 cggcCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCt                      <  2:2337707/150‑1 (MQ=32)
                                                                                                                                   gcCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGc                    <  2:2149502/150‑1 (MQ=32)
                                                                                                                                   gcCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGc                    <  1:1091340/150‑1 (MQ=32)
                                                                                                                                     cAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGccc                  >  1:388559/1‑150 (MQ=32)
                                                                                                                                     cAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGccc                  >  2:3951661/1‑150 (MQ=32)
                                                                                                                                     cAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAGGTGGGCTGccc                  >  2:3687068/1‑150 (MQ=17)
                                                                                                                                     cAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCGCGGCGTATCCCGTACGGCTACCGGGGCGGGCGCGAGGGAACCTCCCAGTAAAACGCGCCCACAGCCCCCCACCCAGTGGGCTGccc                  >  1:4915303/1‑150 (MQ=11)
                                                                                                                                     cAATCGCCGCCGCCCCCGGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGccc                  <  1:4780415/150‑1 (MQ=17)
                                                                                                                                      acgcgccgccccCCCCGCCACCCCGGGGATCTGGCT‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGcccc                 <  1:2277008/147‑1 (MQ=11)
                                                                                                                                      aaTCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGcccc                 >  2:171402/1‑150 (MQ=32)
                                                                                                                                       ctCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCg                <  1:2863311/149‑1 (MQ=17)
                                                                                                                                       aTCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCg                >  1:4227711/1‑150 (MQ=32)
                                                                                                                                       aTCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCg                <  1:2863352/150‑1 (MQ=32)
                                                                                                                                         cgccgccgcCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGc              <  1:356095/150‑1 (MQ=32)
                                                                                                                                           ccgccgccCCCCGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCt            <  1:689146/150‑1 (MQ=17)
                                                                                                                                              ccgccCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCa         <  1:3921706/150‑1 (MQ=32)
                                                                                                                                               cgccCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAg        <  2:2604724/150‑1 (MQ=32)
                                                                                                                                                gccCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAgg       <  1:3713772/150‑1 (MQ=32)
                                                                                                                                                gccCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAgg       <  1:3713730/150‑1 (MQ=32)
                                                                                                                                                    ccTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAAc   <  1:3937937/150‑1 (MQ=32)
                                                                                                                                                    ccTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAAc   <  2:5219437/150‑1 (MQ=32)
                                                                                                                                                    ccTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCGTACGGCTACCGGGTCGGGCGCGCGGTAACATCCCTGTAAAACGCGCCCGCAGCCCCCCTCCCTGTGGGCTGCCCCGGCCGTCAGGGGAc   >  1:4159997/1‑150 (MQ=11)
                                                                                                                                                     cTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACg  <  1:2188186/150‑1 (MQ=32)
                                                                                                                                                     cTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACg  <  1:2572316/150‑1 (MQ=32)
                                                                                                                                                       |                                                                                                                                                    
GATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATT‑GCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACG  >  CP009273/2222535‑2222831

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