Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 503,475 C→A 8.3% W21L (TGG→TTG)  ybaP ← TraB family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273503,4750CA 8.3% 604.2 / 21.2 169W21L (TGG→TTG) ybaPTraB family protein
Reads supporting (aligned to +/- strand):  ref base C (84/71);  new base A (6/8);  total (90/79)
Fisher's exact test for biased strand distribution p-value = 5.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

Only 100 of 171 total aligned reads displayed.

ACATCCGCCTCGACGATCAGCGCATCGGCGTTTTTGAGCTTTTTGAGCAAACGGGTGGGCAGGGGAGCCATATCGTGGCTACCCATATGAATACTGCCAATCAGATGAAAATGGCGATTGCCGGGGAGGGTGATATCGATGGCAGGCCAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGTTTTTACCCGGTACAACAGATCCATACGACCTCCCTTTGTGAAATATCATGCTAGCGCGCGGTGAGGGATGGCGCAAATGGAGTTGGAGACTTGTTTAATGTGTTTGTATG  >  CP009273/503328‑503622
                                                                                                                                                   |                                                                                                                                                   
aCATCCGCCTCGACGATCAGCGCATCGGCGTTTTTGAGCTTTTTGAGCAAACGGGTGGGCAGGGGAGCCATATCGTGGCTACCCATATGAATACTGCCAATCAGATGAAAATGGCGATTGCCGGGGAGGGTGATATCGATGGCAGGcc                                                                                                                                                     <  2:2080103/148‑1 (MQ=255)
 cATCCGCCTCGACGATCAGCGCATCGGCGTTTTTGAGCTTTTTGTGCAAACGGGTGGGCAGGGGAGCCATATCGTGGCTACCCATATGAATACTGCCAATCAGATGAAAATGGCGATTGCCGGGGAGGGTGATATAGATGGCAGGCCa                                                                                                                                                    >  1:290428/1‑148 (MQ=255)
 cATCCGCCTCGACGATCAGCGCATCGGCGTTTTTGAGCTTTTTGAGCAAACGGGTGGGCAGGGGAGCCATATCGTGGCTACCCATATGAATACTGCCAATCAGATGAAAATGGCGATTGCCGGGGAGGGTGATATCGATGGCAGGCCa                                                                                                                                                    <  1:769269/148‑1 (MQ=255)
 cATCCGCCTCGACGATCAGCGCATCGGCGTTTTTGAGCTTTTTGAGCAAACGGGTGGGCAGGGGAGCCATATCGTGGCTACCCATATGAATACTGCCAATCAGATGAAAATGGCGATTGCCGGGGAGGGTGATATCGATGGCAGGCCa                                                                                                                                                    <  1:769239/148‑1 (MQ=255)
 cATCCGCCTCGACGATCAGCGCATCGGCGTTTTTGAGCTTTTTGAGCAAACGGGTGGGCAGGGGAGCCATATCGTGGCTACCCATATGAATACTGCCAATCAGATGAAAATGGCGATTGCCGGGGAGGGTGATATCGATGGCAGGCCa                                                                                                                                                    <  1:768083/148‑1 (MQ=255)
 cATCCGCCTCGACGATCAGCGCATCGGCGTGTTTGAGCTTTTTGAGCAAACGGGTGGGCAGGGGAGCCATATCGTGGCTACCCATATGAATACTGCCAATCAGATGAAAATGGCGATTGCCGGGGAGGGTGATATCGATGGCAGGCCa                                                                                                                                                    >  1:1170665/1‑148 (MQ=255)
        cTCGACGATCAGCGCATCGGCGTTTTTGAGCTTTTTGAGCAAACGGGTGGGCAGGGGAGCCATATCGTGGCTACCCATATGAATACTGCCAATCAGATGAAAATGGCGATTGCCGGGGAGGGTGATATCGATGGCAGGCCAGGTGTaa                                                                                                                                             <  2:2235936/148‑1 (MQ=255)
          cgacgaTCAGCGCATCGGCGTTTTTGAGCTTTTTGAGCAAACGGGTGGGCAGGGGAGCCATATCGTGGCTACCCATATGAATACTGCCAATCAGATGAAAATGGCGATTGCCGGGGAGGGTGATATCGATGGCAGGCCAGGTGTAATg                                                                                                                                           >  2:790643/1‑148 (MQ=255)
           gacgaTCAGCGCATCGGCGTTTTTGAGCTTTTTGAGCAAACGGGTGGGCAGGGGAGCCATATCGTGGCTACCCATATGAATACTGCCAATCAGATGAAAATGGCGATTGCCGGGGAGGGTGATATCGATGGCAGGCCAGGTGTAATGa                                                                                                                                          <  1:1119525/148‑1 (MQ=255)
           gacgaTCAGCGCATCGGCGTTTTTGAGCTTTTTGAGCAAACGGGTGGGCAGGGGAGCCATATCGTGGCTACCCATATGAATACTGCCAATCAGATGAAAATGGCGATTGCCGGGGAGGGTGATATCGATGGCAGGCCAGGTGTAATGa                                                                                                                                          <  1:494400/148‑1 (MQ=255)
            acgaTCAGCGCATCGGCGTTTTTGAGCTTTTTGAGCAAACGGGTGGGCAGGGGAGCCATATCGTGGCTACCCATATGAATACTGCCAATCAGATGAAAATGGCGATTGCCGGGGAGGGTGATATCGATGGCAGGCAAGGTGTAATGAt                                                                                                                                         <  2:3268078/148‑1 (MQ=255)
            acgaTCAGCGCATCGGCGTTTTTGAGCTTTTTGAGCAAACGGGTGGGCAGGGGAGCCATATCGTGGCTACCCATATGAATACTGCCAATCAGATGAAAATGGCGATTGCCGGGGAGGGTGATATCGATGGCAGGCAAGGTGTAATGAt                                                                                                                                         <  2:3268066/148‑1 (MQ=255)
             cgaTCAGCGCATCGGCGTTTTTGAGCTTTTTGAGCAAACGGGTGGGCAGGGGAGCCATATCGTGGCTACCCATATGAATACTGCCAATCAGATGAAAATGGCGATTGCCGGGGAGGGTGATATCGATGGCAGGCCAGGTGTAATGAtt                                                                                                                                        <  1:746545/148‑1 (MQ=255)
              gaTCAGCGCATCGGCGTTTTTGAGCTTTTTGAGCAAACGGGTGGGCAGGGGAGCCATATCGTGGCTACCCATATGAATACTGCCAATCAGATGAAAATGGCGATTGCCGGGGAGGGTGATATCGATGGCAGGCCAGGTGTAATGATTa                                                                                                                                       <  1:616235/148‑1 (MQ=255)
                tCAGCGCATCGGCGTTTTTGAGCTTTTTGAGCAAACGGGTGGGCAGGGGAGCCATATCGTGGCTACCCATATGAATACTGCCAATCAGATGAAAATGGCGATTGCCGGGGAGGGTGATATCGATGGCAGGCCAGGTGTAATGATTAcc                                                                                                                                     >  2:852086/1‑148 (MQ=255)
                          ggCGTTTTTGAGCTTTTTGAGCAAACGGGTGGGCAGGGGAGCCATATCGTGGCTACCCATATGAATACTGCCAATCAGATGAAAATGGCGATTGCCGGGGAGGGTGATATCGATGGCAGGCCACGTGTAATGATTACCGCGCAGcgcg                                                                                                                           >  1:674230/1‑148 (MQ=255)
                             gTTTTTGAGCTTTTTGAGCAAACGGGTGGGCAGGGGAGCCATATCGTGGCTACCCATATGAATACTGCCAATCAGATGAAAATGGCGATTGCCGGGGAGGGTGATATCGATGGCAGGCCAGGTGTAATGATTACCGCGCAGCGCGGcc                                                                                                                        <  2:1169073/148‑1 (MQ=255)
                                tttgagctttttgagcAAACGGGTGGGCAGGGGAGCCATATCGTGGCTACCCATATGAATACTGCCAATCAGATGAAAATGGCGATTGCCGGGGAGGGTGATATCGATGGCAGGCCAGGTGTAATGATTACCGCGCAGCGCGGCCCaa                                                                                                                     >  2:2814573/1‑148 (MQ=255)
                                tttgagctttttgagcAAACGGGTGGGCAGGGGAGCCATATCGTGGCTACCCATATGAATACTGCCAATCAGATGAAAATGGCGATTGCCGGGGAGGGTGATATCGATGGCAGGCCAGGTGTAATGATTACCGCGCAGCGCGGCCCaa                                                                                                                     >  2:2815024/1‑148 (MQ=255)
                                 ttgagctttttgagcAAACGGGTGGGCAGGGGAGCCATATCGTGGCTACCCATATGAATACTGCCAATCAGATGAAAATGGCGATTGCCGGGGAGGGTGATATCGATGGCAGGCAAGGTGTAATGATTACCGCGCAGCGCGGCCCaaa                                                                                                                    <  1:2577883/148‑1 (MQ=255)
                                  tgagctttttgagcAAACGGGTGGGCAGGGGAGCCATATCGTGGCTACCCATATGAATACTGCCAATCAGATGAAAATGGCGATTGCCGGGGAGGGTGATATCGATGGCAGGCCAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAg                                                                                                                   <  2:2143060/148‑1 (MQ=255)
                                   gagctttttgagcAAACGGGTGGGCAGGGGAGCCATATCGTGGCTACCCATATGAATACTGCCAATCAGATGAAAATGGCGATTGCCGGGGAGGGTGATATCGATGGCAGGCCAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGc                                                                                                                  >  1:3024893/1‑148 (MQ=255)
                                    agctttttgagcAAACGGGTGGGCAGGGGAGCCATATCGTGGCTACCCATATGAATACTGCCAATCAGATGAAAATGGCGATTGCCGGGGAGGGTGATATCGATGGCAGGCCAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCg                                                                                                                 >  2:2863078/1‑148 (MQ=255)
                                    agctttttgagcAAACGGGTGGGCAGGGGAGCCATATCGTGGCTACCCATATGAATACTGCCAATCAGATGAAAATGGCGATTGCCGGGGAGGGTGATATCGATGGCAGGCCAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCg                                                                                                                 >  2:2060969/1‑148 (MQ=255)
                                    agctttttgagcAAACGGGTGGGCAGGGGAGCCATATCGTGGCTACCCATATGAATACTGCCAATCAGATGAAAATGGCGATTGCCGGGGAGGGTGATATCGATGGCAGGCAAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCg                                                                                                                 <  2:456028/148‑1 (MQ=255)
                                        ttttgagcAAACGGGTGGGCAGGGGAGCCATATCGTGGCTACCCATATGAATACTGCCAATCAGATGAAAATGGCGATTGCCGGGGAGGGTGATATCGATGGCAGGCAAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGtttt                                                                                                             >  2:1810933/1‑148 (MQ=255)
                                            gagcAAACGGGTGGGCAGGGGAGCCATATCGTGGCTACCCATATGAATACTGCCAATCAGATGAAAATGGCGATTGCCGGGGAGGGTGATATCGATGGCAGGCCAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGTTTTTAcc                                                                                                         <  1:2001250/148‑1 (MQ=255)
                                               cAAACGGGTGGGCAGGGGAGCCATATCGTGGCTACCCATATGAATACTGCCAATCAGATGAAAATGGCGATTGCCGGGGAGGGTGATATCGATGGCAGGCCAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGTTTTTACCCgg                                                                                                      <  1:320262/148‑1 (MQ=255)
                                               cAAACGGGTGGGCAGGGGAGCCATATCGTGGCTACCCATATGAATACTGCCAATCAGATGAAAATGGCGATTGCCGGGGAGGGTGATATCGATGGCAGGCCAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGTTTTTACCCgg                                                                                                      <  1:2550374/148‑1 (MQ=255)
                                               cAAACGGGTGGGCAGGGGAGCCATATCGTGGCTACCCATATGAATACTGCCAATCAGATGAAAATGGCGATTGCCGGGGAGGGTGATATCGATGGCAGGCAAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGTTTTTACCCgg                                                                                                      >  1:3228304/1‑148 (MQ=255)
                                                 aaCGGGTGGGCAGGGGAGCCATATCGTGGCTACCCATATGAATACTGCCAATCAGATGAAAATGGCGATTGCCGGGGAGGGTGATATCGATGGCAGGCCAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGTTTTTACCCGGTa                                                                                                    >  2:1393006/1‑148 (MQ=255)
                                                         ggCAGGGGAGCCATATCGTGGCTACCCATATGAATACTGCCAATCAGATGAAAATGGCGATTGCCGGGGAGGGTGATATCGATGGCAGGCCAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGTTTTTACCCGGTACAACAGAt                                                                                            <  2:2994681/148‑1 (MQ=255)
                                                         ggCAGGGGAGCCATATCGTGGCTACCCATATGAATACTGCCAATCAGATGAAAATGGCGATTGCCGGGGAGGGTGATATCGATGGCAGGCCAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGTTTTTACCCGGTACAACAGAt                                                                                            <  2:2552777/148‑1 (MQ=255)
                                                         ggCAGGGGAGCCATATCGTGGCTACCCATATGAATACTGCCAATCAGATGAAAATGGCGATTGCCGGGGAGGGTGATATCGATGGCAGGCCAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGTTTTTACCCGGTACAACAGAt                                                                                            >  1:1070623/1‑148 (MQ=255)
                                                         ggCAGGGGAGCCATATCGTGGCTACCCATATGAATACTGCCAATCAGATGAAAATGGCGATTGCCGGGGAGGGTGATATCGATGGCAGGCCAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGTTTTTACCCGGTACAACAGAt                                                                                            <  1:2977961/148‑1 (MQ=255)
                                                         ggCAGGGGAGCCATATCGTGGCTACCAATATGAATACTGCCAATCAGATGAAAATGGCGATTGCCGGGGAGGGTGATATCGATGGCAGGCCAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGTTTTTACCCGGTACAACAGAt                                                                                            <  2:2994639/148‑1 (MQ=255)
                                                          gCAGGGGAGCCATATCGTGGCTACCCATATGAATACTGCCAATCAGAGGAAAATGGCGATTGCCGGGGAGGGTGATCTCGAAGGCAGGCCAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGTTTTTACCCGGTACAACAGATc                                                                                           >  2:2388489/1‑148 (MQ=255)
                                                           cAGGGGAGCCATATCGTGGCTACCCATATGAATACTGCCAATCAGATGAAAATGGCGATTGCCGGGGAGGGTGATATCGATGGCAGGCCAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGTTTTTACCCGGTACAACAGATcc                                                                                          <  1:3250831/148‑1 (MQ=255)
                                                           cAGGGGAGCCATATCGTGGCTACCCATATGAATACTGCCAATCAGATGAAAATGGCGATTGCCGGGGAGGGTGATATCGATGGCAGGCCAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGTTTTTACCCGGTACAACAGATcc                                                                                          <  1:1479263/148‑1 (MQ=255)
                                                              gggAGCCATATCGTGGCTACCCATATGAATACTGCCAATCAGATGAAAATGGCGATTGCCGGGGAGGGTGATATCGATGGCAGGCCAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGTTTTTACCCGGTACAACAGATCCATa                                                                                       >  2:469679/1‑148 (MQ=255)
                                                                     atatCGTGGCTACCCATATGAATACTGCCAATCAGATGAAAATGGCGATTGCCGGGGAGGGTGATATCGATGGCAGGCCAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGTTTTTACCCGGTACAACAGATCCATACGACCTc                                                                                >  2:1024248/1‑148 (MQ=255)
                                                                     atatCGTGGCTACCCATATGAATACTGCCAATCAGATGAAAATGGCGATTGCCGGGGAGGGTGATAACGATGGCAGGCCAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGTTTTTACCCGGTACAACAGATCCATACGACCTc                                                                                >  1:376276/1‑148 (MQ=255)
                                                                       atCGTGGCTACCCATATGAATACTGCCAATCAGATGAAAATGGCGATTGCCGGGGAGGGTGATATCGATGGCAGGCCAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGTTTTTACCCGGTACAACAGATCCATACGACCTccc                                                                              >  2:2469188/1‑148 (MQ=255)
                                                                        tCGTGGCTACCCATATGAATACTGCCAATCAGATGAAAATGGCGATTGCCGGGGAGGGTGATATCGATGGCAGGCCAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGTTTTTACCCGGTACAACAGATCCATACGACCTCCCt                                                                             <  2:740308/148‑1 (MQ=255)
                                                                        tCGTGGCTACCCATATGAATACTGCCAATCAGATGAAAATGGCGATTGCCGGGGAGGGTGATATCGATGGCAGGCCAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGTTTTTACCCGGTACAACAGATCCATACGACCTCCCt                                                                             <  2:2209148/148‑1 (MQ=255)
                                                                               tACCCATATGAATACTGCCAATCAGATGAAAATGGCGATTGCCGGGGAGGGTGATATCGATGGCAGGCCAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGTTTTTACCCGGTACAACAGATCCATACGACCTCCCTTTGTGaa                                                                      <  1:1699285/148‑1 (MQ=255)
                                                                                aCCCATATGAATACTGCCAATCAGATGAAAATGGCGATTGCCGGGGAGGGTGATATCGATGGCAGGCCAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGTTTTTACCCGGTACAACAGATCCATACGACCTCCCTTTGTGaaa                                                                     >  2:530704/1‑148 (MQ=255)
                                                                                  ccATATGAATACTGCCAATCAGATGAAAATGGCGATTGCCGGGGAGGGTGATATCGATGGCAGGCCAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGTTTTTACCCGGTACAACAGATCCATACGACCTCCCTTTGTGAAata                                                                   >  2:30369/1‑148 (MQ=255)
                                                                                   cATATGAATACTGCCAATCAGATGAAAATGGCGATTGCCGGGGAGGGTGATATCGATGGCAGGCAAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGTTTTTACCCGGTACAACAGATCCATACGACCTCCCTTTGTGAAatat                                                                  >  1:996946/1‑148 (MQ=255)
                                                                                    atatGAATACTGCCAATCAGATGAAAATGGCGATTGCCGGGGAGGGTGATATCGATGGCAGGCCAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGTTTTTACCCGGTACAACAGATCCATACGACCTCCCTTTGTGAAATATc                                                                 <  2:169580/148‑1 (MQ=255)
                                                                                     tatGAATACTGCCAATCAGATGAAAATGGCGATTGCCGGGGAGGGTGATATCGATGGCAGGCCAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGTTTTTACCCGGTACAACAGATCCATACGACCTCCCTTTGTGAAATATCa                                                                <  2:1716336/148‑1 (MQ=255)
                                                                                      atGAATACTGCCAATCAGATGAAAATGGCGATTGCCGGGGAGGGTGATATCGATGGCAGGCCAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGTTTTTACCCGGTACAACAGATCCATACGACCTCCCTTTGTGAAATATCAt                                                               >  2:1876714/1‑148 (MQ=255)
                                                                                       tGAATACTGCCAATCAGATGAAAATGGCGATTGCCGGGGAGGGTGATATCGATGGCAGGCCAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGTTTTTACCCGGTACAACAGATCCATACGACCTCCCTTTGTGAAATATCATg                                                              <  2:3521532/148‑1 (MQ=255)
                                                                                        gAATACTGCCAATCAGATGAAAATGGCGATTGCCGGGGAGGGTGATATCGATGGCAGGCCAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGTTTTTACCCGGTACAACAGATCCATACGACCTCCCTTTGTGAAATATCATGc                                                             >  2:850459/1‑148 (MQ=255)
                                                                                         aaTACTGCCAATCAGATGAAAATGGCGATTGCCGGGGAGGGTGATATCGATGGCAGGCCAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGTTTTTACCCGGTACAACAGATCCATACGACCTCCCTTTGTGAAATATCATGCt                                                            <  1:79570/148‑1 (MQ=255)
                                                                                          aTACTGCCAATCAGATGAAAATGGCGATTGCCGGGGAGGGTGATATCGATGGCAGGCCAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGTTTTTACCCGGTACAACAGATCCATACGACCTCCCTTTGTGAAATATCATGCTa                                                           <  1:2411875/148‑1 (MQ=255)
                                                                                          aTACTGCCAATCAGATGAAAATGGCGATTGCCGGGGAGGGTGATATCGATGGCAGGCCAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGTTTTTACCCGGTACAACAGATCCATACGACCTCCCTTTGTGAAATATCATGCTa                                                           >  1:2976820/1‑148 (MQ=255)
                                                                                          aTACTGCCAATCAGATGAAAATGGCGATTGCCGGGGAGGGTGATATCGATGGCAGGCCAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGTTTTTACCCGGTACAACAGATCCATACGACCTCCCTTTGTGAAATATCATGCTa                                                           >  2:1115406/1‑148 (MQ=255)
                                                                                            aCTGCCAATCAGATGAAAATGGCGATTGCCGGGGAGGGTGATATCGATGGCAGGCCAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGTTTTTACCCGGTACAACAGATCCATACGACCTCCCTTTGTGAAATATCATGCTAgc                                                         >  1:3155644/1‑148 (MQ=255)
                                                                                            aCTGCCAATCAGATGAAAATGGCGATTGCCGGGGAGGGTGATATCGATGGCAGGCCAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGTTTTTACCCGGTACAACAGATCCATACGACCTCCCTTTGTGAAATATCATGCTAgc                                                         >  1:289543/1‑148 (MQ=255)
                                                                                               gCCAATCAGATGAAAATGGCGATTGCCGGGGAGGGTGATATCGATGGCAGGCCAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGTTTTTACCCGGTACAACAGATCCATACGACCTCCCTTTGTGAAATATCATGCTAgcgcg                                                      >  2:2659244/1‑148 (MQ=255)
                                                                                               gCCAATCAGATGAAAATGGCGATTGCCGGGGAGGGTGATATCGATGGCAGGCCAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGTTTTTACCCGGTACAACAGATCCATACGACCTCCCTTTGTGAAATATCATGCTAgcgcg                                                      >  2:2659060/1‑148 (MQ=255)
                                                                                               gCCAATCAGATGAAAATGGCGATTGCCGGGGAGGGTGATATCGATGGCAGGCCAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGTTTTTACCCGGTACAACAGATCCATACGACCTCCCTTTGTGAAATATCATGCTAgcgcg                                                      >  2:2658967/1‑148 (MQ=255)
                                                                                                  aaTCAGATGAAAATGGCGATTGCCGGGGAGGGTGATATCGATGGCAGGCCAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGTTTTTACCCGGTACAACAGATCCATACGACCTCCCTTTGTGAAATATCATGCTAGCGCGCgg                                                   >  1:3427922/1‑148 (MQ=255)
                                                                                                   aTCAGATGAAAATGGCGATTGCCGGGGAGGGTGATATCGATGGCAGGCCAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGTTTTTACCCGGTACAACAGATCCATACGACCTCCCTTTGTGAAATATCATGCTAGCGCGCGGt                                                  >  1:525293/1‑148 (MQ=255)
                                                                                                   aTCAGATGAAAATGGCGATTGCCGGGGAGGGTGATATCGATGGCAGGCCAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGTTTTTACCCGGTACAACAGATCCATACGACCTCCCTTTGTGAAATATCATGCTAGCGCGCGGt                                                  >  1:2318253/1‑148 (MQ=255)
                                                                                                         tGAAAATGGCGATTGCCGGGGAGGGTGATATCGATGGCAGGCCAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGTTTTTACCCGGTACAACAGATCCATACGACCTCCCTTTGTGAAATATCATGCTAGCGCGCGGTGAGGGa                                            <  2:1249387/148‑1 (MQ=255)
                                                                                                          gAAAATGGCGATTGCCGGGGAGGGTGATATCGATGGCAGGCCAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGTTTTTACCCGGTACAACAGATCCATACGACCTCCCTTTGTGAAATATCATGCTAGCGCGCGGTGAGGGAt                                           >  1:2729546/1‑148 (MQ=255)
                                                                                                              aTGGCGATTGCCGGGGAGGGTGATATCGATGGCAGGCCAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGTTTTTACCCGGTACAACAGATCCATACGACCTCCCTTTGTGAAATATCATGCTAGCGCGCGGTGAGGGATGgcg                                       >  2:1491884/1‑148 (MQ=255)
                                                                                                                ggCGATTGCCGGGGAGGGTGATATCGATGGCAGGCCAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGTTTTTACCCGGTACAACAGATCCATACGACCTCCCTTTGTGAAATATCATGCTAGCGCGCGGTGAGGGATGGCGCa                                     >  1:1963114/1‑148 (MQ=255)
                                                                                                                ggCGATTGCCGGGGAGGGTGATATCGATGGCAGGCCAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGTTNTTACCCGGTACAACAGATCCATACGACCTCCCTTTGTGAAATATCATGCTAGCGCGCGGTGAGGGATGGCGCa                                     <  1:1808085/148‑1 (MQ=255)
                                                                                                                 gCGATTGCCGGGGAGGGTGATATCGATGGCAGGCCAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGTTTTTACCCGGTACAACAGATCCATACGACCTCCCTTTGTGAAATATCATGCTAGCGCGCGGTGAGGGATGGCGCaa                                    >  1:1738244/1‑148 (MQ=255)
                                                                                                                 gCGATTGCCGGGGAGGGTGATATCGATGGCAGGCCAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGTTTTTACCCGGTACAACAGATCCATACGACCTCCCTTTGTGAAATATCATGCTAGCGCGCGGTGAGGGATGGCGCaa                                    >  1:81607/1‑148 (MQ=255)
                                                                                                                  cGATTGCCGGGGAGGGTGATATCGATGGCAGGCCAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGTTTTTACCCGGTACAACAGATCCATACGACCTCCCTTTGTGAAATATCATGCTAGCGCGCGGTGAGGGATGGCGCaaa                                   <  2:3541336/148‑1 (MQ=255)
                                                                                                                   gATTGCCGGGGAGGGTGATATCGATGGCAGGCCAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGTTTTTACCCGGTACAACAGATCCATACGACCTCCCTTTGTGAAATATCATGCTAGCGCGCGGTGAGGGATGGCGCAAAt                                  >  2:341743/1‑148 (MQ=255)
                                                                                                                   gATTGCCGGGGAGGGTGATATCGATGGCAGGCCAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGTTTTTACCCGGTACAACAGATCCATACGACCTCCCTTTGTGAAATATCATGCTAGCGCGCGGTGAGGGATGGCGCAAAt                                  <  2:1070623/148‑1 (MQ=255)
                                                                                                                    aTTGCCGGGGAGGGTGATATCGATGGCAGGCCAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGTTTTTACCCGGTACAACAGATCCATACGACCTCCCTTTGTGAAATATCATGCTAGCGCGCGGTGAGGGATGGCGCAAATg                                 >  1:2269310/1‑148 (MQ=255)
                                                                                                                    aTTGCCGGGGAGGGTGATATCGATGGCAGGCCAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGTTTTTACCCGGAACAACAGATCCATACGACCTCCCTTTGTGAAATATCATGCTAGCGCGCGGTGAGGGATGGCGCAAATg                                 >  1:55479/1‑148 (MQ=255)
                                                                                                                      tGCCGGGGAGGGTGATATCGATGGCAGGCCAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGTTTTTACCCGGTACAACAGATCCATACGACCTCCCTTTGTGAAATATCATGCTAGCGCGCGGTGAGGGATGGCGCAAATGGa                               >  1:260284/1‑148 (MQ=255)
                                                                                                                         cGGGGAGGGTGATATCGATGGCAGGCCAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGTTTTTACCCGGTACAACAGATCCATACGACCTCCCTTTGTGAAATATCATGCTAGCGCGCGGTGAGGGATGGCGCAAATGGAGtt                            <  2:1238736/148‑1 (MQ=255)
                                                                                                                            ggAGGGTGATATCGATGGCAGGCCAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGTTTTTACCCGGTACAACAGATCCATACGACCTCCCTTTGTGAAATATCATGCTAGCGCGCGGTGAGGGATGGCGCAAATGGAGTTGga                         >  1:3328832/1‑148 (MQ=255)
                                                                                                                            ggAGGGTGATATCGATGGCAGGCCAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGTTTTTACCCGGTACAACAGATCCATACGACCTCCCTTTGTGAAATATCATGCTAGCGCGCGGTGAGGGATGGCGCAAATGGAGTTGga                         >  1:477270/1‑148 (MQ=255)
                                                                                                                            ggAGGGTGATATCGATGGCAGGCAAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGTTTTTACCCGGTACAACAGATCCATACGACCTCCCTTTGTGAAATATCATGCTAGCGCGCGGTGAGGGATGGCGCAAATGGAGTTGga                         >  2:587981/1‑148 (MQ=255)
                                                                                                                               gggTGATATCGATGGCAGGCCAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGTTTTTACCCGGTACAACAGATCCATACGACCTCCCTTTGTGAAATATCATGCTAGCGCGCGGTGAGGGATGGCGCAAATGGAGTTGGAGAc                      <  2:706745/148‑1 (MQ=255)
                                                                                                                               gggTGATATCGATGGCAGGCCAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGTTTTTACCCGGTACAACAGATCCATACGACCTCCCTTTGTGAAATATCATGCTAGCGCGCGGTGAGGGATGGCGCAAATGGAGTTGGAGAc                      <  2:2612348/148‑1 (MQ=255)
                                                                                                                               gggTGATATCGATGGCAGGCCAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGTTTTTACCCGGTACAACAGATCCATACGACCTCCCTTTGTGAAATATCATGCTAGCGCGCGGTGAGGGATGGCGCAAATGGAGTTGGAGAc                      <  1:1823776/148‑1 (MQ=255)
                                                                                                                               gggTGATATCGATGGCAGGACAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGTTTTTACCCGGTACAACAGATCCATACGACCTCCCTTTGTGAAATATCATGCTAGCGCGCGGTGAGGGATGGCGCAAATGGAGTTGGAGAc                      <  2:1959023/148‑1 (MQ=255)
                                                                                                                                   gATATCGATGGCAGGCCAGGTGTAATGATTACCGCGCAGCGCGTCCCAAAGCGTTTTTACCCGGTACAACAGATCCATACGACCTCCCTTTGTGAAATATCATGCTAGCGCGCGGTGAGGGATGGCGCAAATGGAGTTGGAGACTTGt                  <  1:38786/148‑1 (MQ=255)
                                                                                                                                   gATATCGATGGCAGGCCAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGTTTTTACCCGGTACAACAGATCCATACGACCTCCCTTTGTGAAATATCATGCTAGCGCGCGGTGAGGGATGGCGCAAATGGAGTTGGAGACTTGt                  <  2:2542309/148‑1 (MQ=255)
                                                                                                                                   gATATCGATGGCAGGCCAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGTTTTTACCCGGTACAACAGATCCATACGACCTCCCTTTGTGAAATATCATGCTAGCGCGCGGTGAGGGATGGCGCAAATGGAGTTGGAGACTTGt                  <  1:1631170/148‑1 (MQ=255)
                                                                                                                                     tatCGATGGCAGGCCAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGTTTTTACCCGGTACAACAGATCCATACGACCTCCCTTTGTGAAATATCATGCTAGCGCGCGGTGAGGGATGGCGCAAATGGAGTTGGAGACTTGttt                <  1:2814573/148‑1 (MQ=255)
                                                                                                                                     tatCGATGGCAGGCCAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGTTTTTACCCGGTACAACAGATCCATACGACCTCCCTTTGTGAAATATCATGCTAGCGCGCGGTGAGGGATGGCGCAAATGGAGTTGGAGACTTGttt                <  1:2815024/148‑1 (MQ=255)
                                                                                                                                        cGATGGCAGGCCAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGTTTTTACCCGGTACAACAGATCCATACGACCTCCCTTTGTGAAATATCATGCTAGCGCGCGGTGAGGGATGGCGCAAATGGAGTTGGAGACTTGTTTAAt             >  1:1351198/1‑148 (MQ=255)
                                                                                                                                        cGATGGCAGGCCAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGTTTTTACCCGGTACAACAGATCCATACGACCTCCCTTTGTGAAATATCATGCTAGCGCGCGGTGAGGGATGGCGCAAATGGAGTTGGAGACTTGTTTAAt             >  1:3337073/1‑148 (MQ=255)
                                                                                                                                           tGGCAGGCAAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGTTTTTACCCGGTACAACAGATCCATACGACCTCCCTTTGTGAAATATCATGCTAGCGCGCGGTGAGGGATGGCGCAAATGGAGTTGGAGACTTGTTTAAtgtg          >  2:1787890/1‑148 (MQ=255)
                                                                                                                                            ggCAGGCCAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGTTTTTACCCGGTACAACAGATCCATACGACCTCCCTTTGTGAAATATCATGCTAGCGCGCGGTGAGGGATGGCGCAAATGGAGTTGGAGACTTGTTTAATgtgt         >  2:1669970/1‑148 (MQ=255)
                                                                                                                                                ggCCAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGTTTTTACCCGGTACAACAGATCCATACGACCTCCCTTTGTGAAATATCATGCTAGCGCGCGGTGAGGGATGGCGCAAATGGAGTTGGAGACTTGTTTAATGTGTTTGt     <  1:2060969/148‑1 (MQ=255)
                                                                                                                                                ggCCAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGTTTTTACCCGGTACAACAGATCCATACGACCTCCCTTTGTGAAATATCATGCTAGCGCGCGGTGAGGGATGGCGCAAATGGAGTTGGAGACTTGTTTAATGTGTTTGt     <  2:2729546/148‑1 (MQ=255)
                                                                                                                                                 gCCAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGTTTTTACCCGGTACAACAGATCCATACGACCTCCCTTTGTGAAATATCATGCTAGCGCGCGGTGAGGGATGGCGCAAATGGAGTTGGAGACTTGTTTAATGTGTTTGTa    >  2:3143334/1‑148 (MQ=255)
                                                                                                                                                 gCCAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGTTTTTACCCGGTACAACAGATCCATACGACCTCCCTTTGTGAAATATCATGCTAGCGCGCGGTGAGGGATGGCGCAAATGGAGTTGGAGACTTGTTTAATGTGTTTGTa    >  2:1369733/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ACATCCGCCTCGACGATCAGCGCATCGGCGTTTTTGAGCTTTTTGAGCAAACGGGTGGGCAGGGGAGCCATATCGTGGCTACCCATATGAATACTGCCAATCAGATGAAAATGGCGATTGCCGGGGAGGGTGATATCGATGGCAGGCCAGGTGTAATGATTACCGCGCAGCGCGGCCCAAAGCGTTTTTACCCGGTACAACAGATCCATACGACCTCCCTTTGTGAAATATCATGCTAGCGCGCGGTGAGGGATGGCGCAAATGGAGTTGGAGACTTGTTTAATGTGTTTGTATG  >  CP009273/503328‑503622

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: