Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,405,730 G→T 100% intergenic (+13/+104) aidB → / ← yjfN DNA alkylation damage repair protein; flavin‑containing DNA binding protein, weak isovaleryl CoA dehydrogenase/DUF1471 family periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,405,7300GT98.5% 1077.1 / 4.9 331intergenic (+13/+104)aidB/yjfNDNA alkylation damage repair protein; flavin‑containing DNA binding protein, weak isovaleryl CoA dehydrogenase/DUF1471 family periplasmic protein
Reads supporting (aligned to +/- strand):  ref base G (3/2);  new base T (138/188);  total (141/190)
Fisher's exact test for biased strand distribution p-value = 6.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 331 total aligned reads displayed.

ATGTTGAAATATGCTTCTCCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCGCCTGTACGCGCCAGTTATCAGGCTG  >  CP009273/4405583‑4405875
                                                                                                                                                   |                                                                                                                                                 
aTGTTGAAATATGCTTCTCCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTAt                                                                                                                                                 <  1:3011131/150‑1 (MQ=255)
aTGTTGAAATATGCTTCTCCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTAt                                                                                                                                                 <  2:1102844/150‑1 (MQ=255)
  gTTGAAATATGCTTCTCCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGc                                                                                                                                               <  2:806664/150‑1 (MQ=255)
  gTTGAAATATGCTTCTCCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGc                                                                                                                                               >  1:983702/1‑150 (MQ=255)
  gTTGAAATATGCTTCTCCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGc                                                                                                                                               >  1:1228847/1‑150 (MQ=255)
         tatGCTTCTCCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTg                                                                                                                                        <  2:3313815/150‑1 (MQ=255)
         tatGCTTCTCCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTg                                                                                                                                        >  1:815120/1‑150 (MQ=255)
           tGCTTCTCCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGt                                                                                                                                      <  1:3051692/150‑1 (MQ=255)
            gCTTCTCCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGtt                                                                                                                                     <  1:2771200/150‑1 (MQ=255)
             cgtctcCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGttt                                                                                                                                    <  1:4717043/148‑1 (MQ=255)
              ttCTCCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTc                                                                                                                                   >  1:4559525/1‑150 (MQ=255)
              ttCTCCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTc                                                                                                                                   <  1:295514/150‑1 (MQ=255)
               tctcCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCg                                                                                                                                  <  2:617803/150‑1 (MQ=255)
                  ccGCCTCTGGCGCTGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGAtt                                                                                                                               <  2:4408206/150‑1 (MQ=255)
                   cGCCACTGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTa                                                                                                                              <  2:2644251/150‑1 (MQ=255)
                     ccAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAAc                                                                                                                            <  2:1757654/150‑1 (MQ=255)
                      cAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACt                                                                                                                           >  2:246132/1‑150 (MQ=255)
                      cAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACt                                                                                                                           >  2:1651151/1‑150 (MQ=255)
                       aaTGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTa                                                                                                                          >  1:4134856/1‑150 (MQ=255)
                       aaTGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTa                                                                                                                          >  1:2987548/1‑150 (MQ=255)
                       aaTGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTa                                                                                                                          >  2:1032501/1‑150 (MQ=255)
                       aaTGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTa                                                                                                                          <  2:3897216/150‑1 (MQ=255)
                        aTGGCGCAGGCGTGGTGTCAGGTGATGTGAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTaa                                                                                                                         <  2:192652/150‑1 (MQ=255)
                          ggCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAAt                                                                                                                       >  1:1078941/1‑150 (MQ=255)
                            cgcAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGa                                                                                                                     >  2:5103276/1‑150 (MQ=255)
                                          cAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAAtt                                                                                                       <  1:5062589/150‑1 (MQ=255)
                                             gTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTa                                                                                                    >  2:2120580/1‑150 (MQ=255)
                                               gATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAgg                                                                                                  >  2:2353610/1‑150 (MQ=255)
                                                  gTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTGATGCTACGCTGGTTTCGATTAACTAAAGGAAATATGTGAAAATTGTTGGCCg                                                                                               >  2:5157698/1‑150 (MQ=255)
                                                    tAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGa                                                                                             >  1:2032754/1‑150 (MQ=255)
                                                     aGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGAc                                                                                            >  2:1760758/1‑150 (MQ=255)
                                                     aGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGAc                                                                                            >  1:88548/1‑150 (MQ=255)
                                                     aGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGAc                                                                                            >  1:2889018/1‑150 (MQ=255)
                                                      gATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACa                                                                                           <  1:668632/150‑1 (MQ=255)
                                                        tACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAg                                                                                         <  1:3519794/150‑1 (MQ=255)
                                                         aCGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAgg                                                                                        <  1:1665536/150‑1 (MQ=255)
                                                         aCGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAgg                                                                                        <  2:2229531/150‑1 (MQ=255)
                                                         aCGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAgg                                                                                        <  2:306205/150‑1 (MQ=255)
                                                              cggcggcgTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTc                                                                                   <  2:384688/150‑1 (MQ=255)
                                                              cggcggcgTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTc                                                                                   <  2:1614828/150‑1 (MQ=255)
                                                              cggcggcgTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTc                                                                                   >  1:3311326/1‑150 (MQ=255)
                                                               ggcggcgTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTcg                                                                                  <  1:3579856/150‑1 (MQ=255)
                                                               ggcggcgTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTcg                                                                                  >  2:4147627/1‑150 (MQ=255)
                                                                gcggcgTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTcgc                                                                                 <  2:3117698/150‑1 (MQ=255)
                                                                  ggcgTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCgcgc                                                                               >  1:4783367/1‑150 (MQ=255)
                                                                  ggcgTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCgcgc                                                                               >  2:2157907/1‑150 (MQ=255)
                                                                     gTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGcgccgc                                                                            >  1:4535173/1‑150 (MQ=255)
                                                                         ggTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATcc                                                                        <  1:3085683/150‑1 (MQ=255)
                                                                         ggTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATcc                                                                        <  1:1889187/150‑1 (MQ=255)
                                                                         ggTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATcc                                                                        <  2:1521043/150‑1 (MQ=255)
                                                                          gTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCg                                                                       >  1:4199859/1‑150 (MQ=255)
                                                                          gTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCg                                                                       <  1:139387/150‑1 (MQ=255)
                                                                           ttGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCgg                                                                      <  2:2553813/150‑1 (MQ=255)
                                                                            tGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGc                                                                     <  1:4406772/150‑1 (MQ=255)
                                                                             gTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAGGGCGCTCGCGCCGCATCCGGCa                                                                    <  1:5124773/150‑1 (MQ=255)
                                                                              tCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCAt                                                                   <  1:1042352/150‑1 (MQ=255)
                                                                               cAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCAtt                                                                  <  1:3835412/150‑1 (MQ=255)
                                                                                agagCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTg                                                                 <  2:431401/150‑1 (MQ=255)
                                                                                 gagCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGt                                                                <  1:4051951/150‑1 (MQ=255)
                                                                                    cAGATACAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCa                                                             >  1:4481259/1‑150 (MQ=255)
                                                                                     aGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCAt                                                            >  1:2336803/1‑150 (MQ=255)
                                                                                     aGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCAt                                                            <  1:1317365/150‑1 (MQ=255)
                                                                                     aGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCAGTGTTCAt                                                            >  2:1935856/1‑150 (MQ=255)
                                                                                          cAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCt                                                       >  2:1261966/1‑150 (MQ=255)
                                                                                            gAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGt                                                     >  2:3206327/1‑150 (MQ=255)
                                                                                              aTGATTTATGGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATCCGACTTATGCGACGCGGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAg                                                   <  2:4597877/150‑1 (MQ=255)
                                                                                               tGATTTATTGCTGCGGGCGACGGGGGGAGGGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGa                                                  >  2:293433/1‑150 (MQ=255)
                                                                                                gATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGAt                                                 >  1:4212581/1‑150 (MQ=255)
                                                                                                  tttATTGCGGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGa                                               <  2:1433827/150‑1 (MQ=255)
                                                                                                    tagtgctgcGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTATATGaca                                             <  2:3221171/147‑1 (MQ=255)
                                                                                                    tATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGaca                                             >  1:587418/1‑150 (MQ=255)
                                                                                                    tATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGaca                                             <  1:3151133/150‑1 (MQ=255)
                                                                                                    tATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGaca                                             <  2:3509701/150‑1 (MQ=255)
                                                                                                        gctgcGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCa                                         >  1:2731821/1‑150 (MQ=255)
                                                                                                         ctgcGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCAt                                        <  2:2520737/150‑1 (MQ=255)
                                                                                                         ctgcGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCAt                                        >  1:72132/1‑150 (MQ=255)
                                                                                                            cGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCa                                     <  2:4913664/150‑1 (MQ=255)
                                                                                                            cGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCa                                     <  1:4900787/150‑1 (MQ=255)
                                                                                                             gggCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCAt                                    <  2:1932290/150‑1 (MQ=255)
                                                                                                                cGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACa                                 <  1:2273073/150‑1 (MQ=255)
                                                                                                                cGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACa                                 >  1:1803739/1‑150 (MQ=255)
                                                                                                                cGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACa                                 <  1:2151290/150‑1 (MQ=255)
                                                                                                                  aCGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGt                               >  2:2173787/1‑150 (MQ=255)
                                                                                                                   cGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTa                              <  1:3412408/150‑1 (MQ=255)
                                                                                                                      ggggAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCg                           <  1:3707153/150‑1 (MQ=255)
                                                                                                                       gggAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCGc                          >  2:660561/1‑150 (MQ=255)
                                                                                                                        ggAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCGcc                         <  2:564511/150‑1 (MQ=255)
                                                                                                                            tgtgtgtgtAAGCGTATACGACTTATGCGACGCTGGTTTCGATTTACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCGCCTGTa                     <  2:2767756/150‑1 (MQ=255)
                                                                                                                            tgtgtgtgtAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCGCCTGTa                     <  2:964318/150‑1 (MQ=255)
                                                                                                                               gtgtgtAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCGCCTGTAcgc                  <  2:2089814/150‑1 (MQ=255)
                                                                                                                                 gtgtAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCGCCTGTACgcgc                >  1:1936568/1‑150 (MQ=255)
                                                                                                                                  tgtAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCGCCTGTACGCGcc               <  1:2456821/150‑1 (MQ=255)
                                                                                                                                   gtAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCGCCTGTACGCGCCa              >  2:4970404/1‑150 (MQ=255)
                                                                                                                                    tAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCGCCTGTACGCGCCAg             >  2:4555282/1‑150 (MQ=255)
                                                                                                                                         gTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCGCCTGTACGCGCCAGTTATc        >  2:3986951/1‑150 (MQ=255)
                                                                                                                                         gTATACGACTGATGCTACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCGCCTGTACGCGCCAGTTATc        >  1:4093988/1‑150 (MQ=255)
                                                                                                                                          tataCGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCGCCTGTACGCGCCAGTTATCa       <  2:409564/150‑1 (MQ=255)
                                                                                                                                              cGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGACAAGTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCGCCTGTACGCGCCAGTTATCAGgct   <  2:4656233/150‑1 (MQ=255)
                                                                                                                                               gACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCGCCTGTACGCGCCAGTTATCAGgctg  >  2:2805792/1‑150 (MQ=255)
                                                                                                                                               gACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGAAGGACAAGGCGCTCGCGACGCATCCGGCATTGTTCATTAGCTGTAGATGACACTCATGCATACAGTATCGCCTGTACGCGACAGTTATCAGgctg  >  1:3329122/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
ATGTTGAAATATGCTTCTCCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCGCCTGTACGCGCCAGTTATCAGGCTG  >  CP009273/4405583‑4405875

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