Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,405,730 G→T 90.5% 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,7300GT90.5% 448.4 / 38.9 173intergenic (+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 (7/10);  new base T (62/94);  total (69/104)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 173 total aligned reads displayed.

GTTGAAATATGCTTCTCCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCGCCTGTACGCGCCAGTTATCAGGCT  >  CP009273/4405585‑4405874
                                                                                                                                                 |                                                                                                                                                
  tGAAATATGCTTCTCCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGc                                                                                                                                              <  2:1144812/148‑1 (MQ=255)
      atatGCTTCTCCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACg                                                                                                                                          <  2:2881981/148‑1 (MQ=255)
      atatGCTTCTCCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACg                                                                                                                                          <  1:1378786/148‑1 (MQ=255)
      atatGCTTCTCCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACg                                                                                                                                          <  2:1232891/148‑1 (MQ=255)
         tGCTTCTCCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTg                                                                                                                                       >  1:1381352/1‑148 (MQ=255)
           cTTCTCCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGt                                                                                                                                     <  1:2570665/148‑1 (MQ=255)
           cTTCTCCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGt                                                                                                                                     <  2:2283033/148‑1 (MQ=255)
           cTTCTCCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGt                                                                                                                                     <  2:1762754/148‑1 (MQ=255)
            ttCTCCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGtt                                                                                                                                    <  1:3172454/148‑1 (MQ=255)
              ctcCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTc                                                                                                                                  >  1:241846/1‑148 (MQ=255)
                ccGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGa                                                                                                                                <  2:555379/148‑1 (MQ=255)
                 cGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGAt                                                                                                                               <  2:531512/148‑1 (MQ=255)
                   ccAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTa                                                                                                                             >  2:1774663/1‑148 (MQ=255)
                     aaTGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAAc                                                                                                                           >  2:2805495/1‑148 (MQ=255)
                      aTGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACt                                                                                                                          >  2:2550786/1‑148 (MQ=255)
                      aTGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACt                                                                                                                          <  2:90196/148‑1 (MQ=255)
                       tGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTa                                                                                                                         >  1:3131195/1‑148 (MQ=255)
                        ggCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTaa                                                                                                                        <  1:539714/148‑1 (MQ=255)
                         gcgcAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTaaa                                                                                                                       >  2:2898514/1‑148 (MQ=255)
                           gggggCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATg                                                                                                                     >  1:1247663/4‑148 (MQ=255)
                              ggCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGaaa                                                                                                                  <  1:999417/148‑1 (MQ=255)
                               gCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAat                                                                                                                 >  2:626405/1‑148 (MQ=255)
                                 gtggtgTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAatat                                                                                                               <  1:2482761/148‑1 (MQ=255)
                                  tggtgTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATAtg                                                                                                              <  2:163739/148‑1 (MQ=255)
                                       tCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGaaa                                                                                                         <  1:864723/148‑1 (MQ=255)
                                       tCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTGATGCTACGCTGGTTTCGATTAACTAAATGAAATATGTGaaa                                                                                                         <  2:1098639/148‑1 (MQ=255)
                                       tCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTGATGCTACGCTGGTTTCGATTAACTAAATGAAATATGTGaaa                                                                                                         <  2:1098505/148‑1 (MQ=255)
                                       tCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTGATGCTACGCTGGTTTCGATTAACTAAATGAAATATGTGaaa                                                                                                         <  2:2751383/148‑1 (MQ=255)
                                        cAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGaaaa                                                                                                        <  1:10827/148‑1 (MQ=255)
                                         aGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAAt                                                                                                       <  1:2716551/148‑1 (MQ=255)
                                         aGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAAt                                                                                                       <  1:42140/148‑1 (MQ=255)
                                          ggTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAAtt                                                                                                      <  1:2314168/148‑1 (MQ=255)
                                          ggTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAAtt                                                                                                      <  1:2314147/148‑1 (MQ=255)
                                          ggTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAAtt                                                                                                      <  1:3077648/148‑1 (MQ=255)
                                           gTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTg                                                                                                     <  1:1836115/148‑1 (MQ=255)
                                           gTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTGATGCTACGCTGGCTTCGATTAACTAAATGAAATATGTGAAAATTg                                                                                                     >  2:910988/1‑148 (MQ=255)
                                            tGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGt                                                                                                    <  1:1235273/148‑1 (MQ=255)
                                             gATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTa                                                                                                   >  1:606418/1‑148 (MQ=255)
                                             gATGTTAGAGACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTa                                                                                                   <  2:789288/148‑1 (MQ=255)
                                              aTGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAg                                                                                                  <  2:1796565/148‑1 (MQ=255)
                                              aTGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAg                                                                                                  <  2:3192372/148‑1 (MQ=255)
                                              aTGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTGATGCTACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAg                                                                                                  >  1:1449422/1‑148 (MQ=255)
                                                gTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGc                                                                                                <  2:748200/148‑1 (MQ=255)
                                                 ttAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTGATGCTACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGcc                                                                                               <  1:1991392/148‑1 (MQ=255)
                                                  tAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCg                                                                                              <  1:998191/148‑1 (MQ=255)
                                                  tAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCg                                                                                              <  1:2926983/148‑1 (MQ=255)
                                                    gATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGa                                                                                            >  1:1595945/1‑148 (MQ=255)
                                                     aTACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAAt                                                                                                       >  1:1672616/1‑136 (MQ=255)
                                                     aTACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAAt                                                                                                       <  2:1672616/136‑1 (MQ=255)
                                                        cgcgcgGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAg                                                                                        <  1:976305/148‑1 (MQ=255)
                                                        cgcgcgGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAg                                                                                        >  1:2662264/1‑148 (MQ=255)
                                                        cgcgcgGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAg                                                                                        >  1:2663946/1‑148 (MQ=255)
                                                         gcgcgGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAgg                                                                                       <  1:1451235/148‑1 (MQ=255)
                                                          cgcgGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGgc                                                                                      >  1:3043157/1‑148 (MQ=255)
                                                            cggcggcgTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGgcgc                                                                                    <  1:1409270/148‑1 (MQ=255)
                                                            cggcggcgTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGgcgc                                                                                    <  2:26953/148‑1 (MQ=255)
                                                            cggcggcgTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGgcgc                                                                                    <  2:1872911/148‑1 (MQ=255)
                                                               cggcgTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTGAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTcg                                                                                 <  1:3226456/148‑1 (MQ=255)
                                                                  cgTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCgcgc                                                                              >  1:2860254/1‑148 (MQ=255)
                                                                    tACGGTTGTCATAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGcgccg                                                                            <  2:1892042/148‑1 (MQ=255)
                                                                       ggTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCAt                                                                         <  2:2239483/148‑1 (MQ=255)
                                                                         ttGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATcc                                                                       >  2:727723/1‑148 (MQ=255)
                                                                         ttGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATcc                                                                       >  2:1483998/1‑148 (MQ=255)
                                                                           gacAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCgg                                                                     <  2:2074368/146‑1 (MQ=255)
                                                                            tCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGc                                                                    >  1:97185/1‑148 (MQ=255)
                                                                            tCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGc                                                                    <  1:2887671/148‑1 (MQ=255)
                                                                            tCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGc                                                                    >  1:97160/1‑148 (MQ=255)
                                                                            tCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGc                                                                    >  1:799146/1‑148 (MQ=255)
                                                                             cAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCa                                                                   >  2:41560/1‑148 (MQ=255)
                                                                              agagCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCAt                                                                  >  1:111294/1‑148 (MQ=255)
                                                                                  cAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGtt                                                              <  1:127632/148‑1 (MQ=255)
                                                                                    gATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCa                                                            <  1:225500/148‑1 (MQ=255)
                                                                                    gATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATGGTTCa                                                            >  2:2287840/1‑148 (MQ=255)
                                                                                     aTCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCAt                                                           <  1:2949629/148‑1 (MQ=255)
                                                                                      tCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTGATGCTACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCAtt                                                          >  2:2524179/1‑148 (MQ=255)
                                                                                             tGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATACGGCATTGTTCATTTGCTGTa                                                   >  1:1733487/1‑148 (MQ=255)
                                                                                              gATTTAGTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAg                                                  <  2:354421/148‑1 (MQ=255)
                                                                                               aTTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGa                                                 >  1:2111214/1‑148 (MQ=255)
                                                                                                    ttGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGaca                                            <  2:3268630/148‑1 (MQ=255)
                                                                                                    ttGCTGCGGGCAACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGaca                                            >  2:2315081/1‑148 (MQ=255)
                                                                                                     tgctgcGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGacac                                           <  2:2050431/148‑1 (MQ=255)
                                                                                                       ctgcGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTc                                         <  2:703752/148‑1 (MQ=255)
                                                                                                       ctgcGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTGATGCGACGCTGGTTTCGAGTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTc                                         <  2:284574/148‑1 (MQ=255)
                                                                                                              cGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATa                                  >  2:2013106/1‑148 (MQ=255)
                                                                                                              cGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATa                                  <  1:271957/148‑1 (MQ=255)
                                                                                                              cGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATa                                  <  1:3138426/148‑1 (MQ=255)
                                                                                                                aCGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACa                                >  1:2177301/1‑148 (MQ=255)
                                                                                                                aCGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACa                                <  1:474039/148‑1 (MQ=255)
                                                                                                                aCGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACa                                >  1:355879/1‑148 (MQ=255)
                                                                                                                    ggggAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTAt                            <  2:1461669/148‑1 (MQ=255)
                                                                                                                       gAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCGc                         >  1:1724128/1‑148 (MQ=255)
                                                                                                                        aGTGTGTGTGTAAGCGTATACGACTGATGCTACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCGcc                        <  1:1270755/148‑1 (MQ=255)
                                                                                                                          tgtgtgtgtAAGCGTATACGACTGATGCTACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCGCCTg                      >  2:902501/1‑148 (MQ=255)
                                                                                                                           gtgtgtgtAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCGCCTGt                     <  2:2453367/148‑1 (MQ=255)
                                                                                                                            tgtgtgtAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCGCCTGTa                    <  2:3057813/148‑1 (MQ=255)
                                                                                                                                   aaGCGTATATGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCGCCTGTACGCGCCa             >  1:1053052/1‑148 (MQ=255)
                                                                                                                                        tataCGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCGCCTGTACGCGCCAGTTAt        <  2:926473/148‑1 (MQ=255)
                                                                                                                                         ataCGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCGCCTGTACGCGCCAGTTATc       >  1:2538361/1‑148 (MQ=255)
                                                                                                                                          taCGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCGCCTGTACGCGCCAGTTATCa      <  2:2663946/148‑1 (MQ=255)
                                                                                                                                          taCGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCGCCTGTACGCGCCAGTTATCa      <  2:2662264/148‑1 (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: