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,7300GT97.6% 643.0 / 6.0 206intergenic (+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 (2/3);  new base T (86/115);  total (88/118)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 206 total aligned reads displayed.

GTTGAAATATGCTTCTCCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCGCCTGTACGCGCCAGTTATCAGGCT  >  CP009273/4405585‑4405874
                                                                                                                                                 |                                                                                                                                                
    aaaTATGCTTCTCCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGa                                                                                                                                            >  1:1621297/1‑148 (MQ=255)
      atatGCTTCTCCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACg                                                                                                                                          >  2:3004915/1‑148 (MQ=255)
      atatGCTTCTCCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACg                                                                                                                                          >  2:1540932/1‑148 (MQ=255)
       tatGCTTCTCCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGc                                                                                                                                         >  1:32575/1‑148 (MQ=255)
       tatGCTTCTCCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGc                                                                                                                                         <  1:1301843/148‑1 (MQ=255)
             tctcCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGttt                                                                                                                                   >  1:2170405/1‑148 (MQ=255)
             tctcCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGttt                                                                                                                                   >  2:2861524/1‑148 (MQ=255)
             tctcCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGttt                                                                                                                                   <  1:1786065/148‑1 (MQ=255)
              ctcCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTc                                                                                                                                  >  1:1859484/1‑148 (MQ=255)
                ccGCCTATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGa                                                                                                                                <  2:614608/148‑1 (MQ=255)
                ccGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGa                                                                                                                                >  2:1635158/1‑148 (MQ=255)
                ccGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGa                                                                                                                                <  2:1479420/148‑1 (MQ=255)
                 cGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGAt                                                                                                                               <  1:174696/148‑1 (MQ=255)
                   ccAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTa                                                                                                                             <  1:2351463/148‑1 (MQ=255)
                    cAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTaa                                                                                                                            <  2:1692448/148‑1 (MQ=255)
                           gcAGGCGTGGTGTCAGGTGATGTTAGATTCGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGGGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATg                                                                                                                     <  2:2447012/148‑1 (MQ=255)
                             aGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGaa                                                                                                                   <  2:837167/148‑1 (MQ=255)
                              ggCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGaaa                                                                                                                  <  1:2561139/148‑1 (MQ=255)
                              ggCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGaaa                                                                                                                  <  1:3043074/148‑1 (MQ=255)
                               gCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAat                                                                                                                 >  1:720674/1‑148 (MQ=255)
                                   ggtgTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATAtgt                                                                                                             >  1:3523508/1‑148 (MQ=255)
                                    gtgtCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATAtgtg                                                                                                            <  1:606667/148‑1 (MQ=255)
                                    gtgtCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATAtgtg                                                                                                            <  1:3395224/148‑1 (MQ=255)
                                    gtgtCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATAtgtg                                                                                                            <  2:551072/148‑1 (MQ=255)
                                       tCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGaaa                                                                                                         >  1:446251/1‑148 (MQ=255)
                                       tCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGaaa                                                                                                         <  1:1217372/148‑1 (MQ=255)
                                        cAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGaaaa                                                                                                        <  2:2112946/148‑1 (MQ=255)
                                            tGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGt                                                                                                    <  1:621554/148‑1 (MQ=255)
                                             gATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTa                                                                                                   <  2:1554420/148‑1 (MQ=255)
                                             gATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTa                                                                                                   <  1:2564476/148‑1 (MQ=255)
                                             gATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTa                                                                                                   <  1:2115604/148‑1 (MQ=255)
                                               tGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAgg                                                                                                 >  1:199180/1‑148 (MQ=255)
                                                  tAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTCACTAAATGAAATATGTGACAATTGTAGGCCg                                                                                              >  2:2913796/1‑148 (MQ=255)
                                                  tAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCg                                                                                              >  2:2689483/1‑148 (MQ=255)
                                                     aTACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTGATGCTACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGAc                                                                                           >  1:3562540/1‑148 (MQ=255)
                                                      tACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACa                                                                                          >  2:2160129/1‑148 (MQ=255)
                                                      tACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACa                                                                                          <  1:2041619/148‑1 (MQ=255)
                                                        cgcgcgGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGTCGGACAAg                                                                                        <  1:2350469/148‑1 (MQ=255)
                                                           gcggcggcgTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGgcg                                                                                     <  1:84668/148‑1 (MQ=255)
                                                             ggcggcgTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCt                                                                                   <  1:2380970/148‑1 (MQ=255)
                                                                    tACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGcgccg                                                                            >  2:3385249/1‑148 (MQ=255)
                                                                         ttGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATcc                                                                       >  2:3416438/1‑148 (MQ=255)
                                                                         ttGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATAGGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATcc                                                                       >  2:721064/1‑148 (MQ=255)
                                                                            tCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGc                                                                    <  1:2780232/148‑1 (MQ=255)
                                                                             cAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCa                                                                   <  1:2028227/148‑1 (MQ=255)
                                                                               gagCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCAtt                                                                 <  1:1601888/148‑1 (MQ=255)
                                                                                 gCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAAAGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCAGCATCCGGCAgtgt                                                               >  2:1278998/1‑148 (MQ=255)
                                                                                  cAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCAATGtt                                                              >  1:3081495/1‑148 (MQ=255)
                                                                                     aTCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCAt                                                           >  1:2222569/1‑148 (MQ=255)
                                                                                     aTCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCAt                                                           >  1:2889427/1‑148 (MQ=255)
                                                                                       ccAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCAttt                                                         <  2:1379090/148‑1 (MQ=255)
                                                                                        cAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTg                                                        <  1:1776823/148‑1 (MQ=255)
                                                                                         aGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTACGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGc                                                       <  1:784095/148‑1 (MQ=255)
                                                                                         aGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGc                                                       <  1:2540493/148‑1 (MQ=255)
                                                                                           aaTGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTg                                                     >  1:2904712/1‑148 (MQ=255)
                                                                                            aTGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGt                                                    <  1:604552/148‑1 (MQ=255)
                                                                                             tGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTa                                                   >  2:2664671/1‑148 (MQ=255)
                                                                                              gATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAg                                                  <  2:2827014/148‑1 (MQ=255)
                                                                                              gATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAg                                                  >  1:1992117/1‑148 (MQ=255)
                                                                                               aTTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGa                                                 >  1:947149/1‑148 (MQ=255)
                                                                                                 ttATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATg                                               >  1:3160646/1‑148 (MQ=255)
                                                                                                 ttATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATg                                               <  1:581208/148‑1 (MQ=255)
                                                                                                   aTTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGac                                             <  1:2576184/148‑1 (MQ=255)
                                                                                                       ctgcGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTc                                         <  2:1495702/148‑1 (MQ=255)
                                                                                                       ctgcGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTc                                         <  2:1495656/148‑1 (MQ=255)
                                                                                                       ctgcGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGATGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTc                                         <  2:2429627/148‑1 (MQ=255)
                                                                                                        tgcGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATTTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCa                                        <  2:210465/148‑1 (MQ=255)
                                                                                                        tgcGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCa                                        <  1:2278296/148‑1 (MQ=255)
                                                                                                           gggCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGc                                     <  1:3306236/148‑1 (MQ=255)
                                                                                                             gCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTATATGACACTCATGCAt                                   <  2:3213913/148‑1 (MQ=255)
                                                                                                             gCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCAt                                   <  1:1120923/148‑1 (MQ=255)
                                                                                                              cGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATa                                  <  2:463273/148‑1 (MQ=255)
                                                                                                               gACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATAc                                 <  2:2483289/148‑1 (MQ=255)
                                                                                                               gACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATAc                                 <  2:2565248/148‑1 (MQ=255)
                                                                                                               gACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATAc                                 <  2:2565407/148‑1 (MQ=255)
                                                                                                                aCGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACa                                >  1:980185/1‑148 (MQ=255)
                                                                                                                  ggggggAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGt                              <  1:2441763/148‑1 (MQ=255)
                                                                                                                   gggggAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTa                             <  2:3108621/148‑1 (MQ=255)
                                                                                                                       gAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCGc                         >  1:1525229/1‑148 (MQ=255)
                                                                                                                       gAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCGc                         <  1:3474640/148‑1 (MQ=255)
                                                                                                                                tgtAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCGCCTGTAcgcg                >  1:3048329/1‑148 (MQ=255)
                                                                                                                                 gtAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCGCCTGTACgcgc               <  2:2071392/148‑1 (MQ=255)
                                                                                                                                 gtAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCGCCTGTACgcgc               <  2:1192425/148‑1 (MQ=255)
                                                                                                                                 gtAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCGCCTGTACgcgc               <  2:340531/148‑1 (MQ=255)
                                                                                                                                  tAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCGCCTGTACGCGcc              <  2:3085051/148‑1 (MQ=255)
                                                                                                                                      cGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCGCCTGTACGCGCCAGtt          <  2:2170413/148‑1 (MQ=255)
                                                                                                                                       gTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCGCCTGTACGCGCCAGTTa         <  1:798385/148‑1 (MQ=255)
                                                                                                                                         ataCGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCGCCTGTACGCGCCAGTTATc       >  1:2598447/1‑148 (MQ=255)
                                                                                                                                         ataCGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCGCCTGTACGCGCCAGTTATc       >  2:2438940/1‑148 (MQ=255)
                                                                                                                                          taCGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGTGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCGCCTGTACGCGCCAGTTATCa      >  1:3107262/1‑148 (MQ=255)
                                                                                                                                          taCGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCGCCTGTACGCGCCAGTTATCa      <  1:3091400/148‑1 (MQ=255)
                                                                                                                                          taCGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCGCCTGTACGCGCCAGTTATCa      >  1:908178/1‑148 (MQ=255)
                                                                                                                                          taCGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCGCCTGTACGCGCCAGTTATCa      >  2:966260/1‑148 (MQ=255)
                                                                                                                                          taCGACTGATGCTACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCGCCTGTACGCGCCAGTTATCa      <  2:3562540/148‑1 (MQ=255)
                                                                                                                                           aCGACTGATGCTACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCGCCTGTACGCGCCAGTTATCAg     >  1:1379740/1‑148 (MQ=255)
                                                                                                                                            cGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCGCCTGTACGCGCCAGTTATCAgg    <  1:825194/148‑1 (MQ=255)
                                                                                                                                            cGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCGCCTGTACGCGCCAGTTATCAgg    <  2:3523508/148‑1 (MQ=255)
                                                                                                                                             gACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCGCCTGTACGCGCCAGTTATCAGgc   <  1:1017394/148‑1 (MQ=255)
                                                                                                                                              aCTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCGCCTGTACGCGCCAGTTATCAGgct  >  2:364999/1‑148 (MQ=255)
                                                                                                                                              aCTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCGCCTGTACGCGCCAGTTATCAGgct  >  2:1830180/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GTTGAAATATGCTTCTCCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCGCCTGTACGCGCCAGTTATCAGGCT  >  CP009273/4405585‑4405874

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: