Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,364,401 T→A 13.1% intergenic (‑350/+45) yjeJ ← / ← epmB uncharacterized protein/EF‑P‑Lys34 lysylation protein; weak lysine 2,3‑aminomutase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,364,4010TA13.1% 269.8 / 23.2 94intergenic (‑350/+45)yjeJ/epmBuncharacterized protein/EF‑P‑Lys34 lysylation protein; weak lysine 2,3‑aminomutase
Reads supporting (aligned to +/- strand):  ref base T (44/37);  new base A (7/6);  total (51/43)
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

GTTAAAATAACATAAATAGCCGTACTCACTCTATATAAACAGAAAAGAAAGGTTAATTATTGGTGTTAGCTATATAAAAAGTAACAACAGCAATGCATATGAATAATTAAAGGTTATGCAACGGGCAAAGATTAAACTTCCCAGTTTAATGACAGGTTGCTGTACTTGCGTTTGTTGCTCTTGCCGGATGCAGCGTAAACGCCTTATCCGGCCTACAGTTCATTGCAAATCCATTCAATTTCAATGCACGATGTAGGCCTGATAAGCGTAGCGAATCAGGCAATTTTAATGTTTAACTTCCCTGTTTAATCAGTTACTGCTGGCGTAGCTGGAGATCCAGCGGCGTTTTGCTGGGTTCGCCGCCAATTTCTCGCGCCAGTTTCGGCACCAGATATCCCGACACCAGTGTCAGCAACTCACGCATAATCTGCCGTGCTTCGTCATCACTCACCATAAAATGCGCCGCGCCCTGTACTTTATCGAGCACATGCAGGTAATACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCC  >  CP009273/4364133‑4364674
                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                 
gttAAAATAAAATAAATAGCCGTACTCACTCTATATAAACAAAAAAGAAAGGTTAATTATTGGTGTTAGATATATAAAAAGTAACAACAGCAATGCATATGAATAATTAAAGGTTATGCAACGGGCAAAGATTAAACTTCCCAGTTTAATGACAGGTTGCTGTACTTGCGTTTGTTGCTCTTGCCGGATGCAGCGTAAACGCCCTATCCGGCCTACAGTTCATTGCAAATCCATTCAATTTCAATGCACGATGTAGGCCTGATAAGCGTAGCGaa                                                                                                                                                                                                                                                                             <  2:673295/275‑1 (MQ=255)
          cATAAATAGCCGTACTCACTCTATATAAACAGAAAAGAAAGGTTAATAATTGGTGTTAGCTATATAAAAAGTAACAACAGCAATGCATATGAATAATTAAAGGTTATGCAACGGGCAAAGATTAAACTTCACAGTTTAATGACAGGTTGCTGTACTTGCGTTTGTTGCTCTTGCCGGATGCAGCGTAAACGCCTTGTCCGGCCTACAGTTCATTGCAAATCCATTCAATTTCAATGCACGATGTAGGCCTGATAAGCGTAg                                                                                                                                                                                                                                                                                 >  1:673304/1‑261 (MQ=255)
                             tctATATAAACAGAAAAGAAAGGTTAATTATTGGTGTTAGCTATATAAAAGTAACAACAGCAATGCATATGAATAATTAAAGGTTATGCAACGGGCAAAGATTAAACTTCCCAGTTTAATGACAGGTTGCTGTACTTGCGTTTGTTGCTCTTGCCGGATGCAGCGTAAACGCCTTATCCGGCCNACAGTTCATTGCAAATCCATTCAATTTCAATGCACGATGTAGGCCTGATAAGCGTAGCGAATCAGGCAATTTTAATGTTTAACTTCCCTg                                                                                                                                                                                                                                                >  1:549573/1‑274 (MQ=255)
                                           aaagaaagGTTAATTATTGGTGTTAGCTATATAAAAAGTAACAACAGCAATGCATATGAATAATTAAAGGTTATGCAACGGGCAAAGATTAAACTTCCCAGTTTAATGACAGGTTGCTGTACTTGCGTTTGTTGCTCTTGCCGGATGCAGCGTAAACGCCTTATCCGGCCTACAGTTCATTGCAAATCCATTCAATTTCAATGCACGATGTAGGCCTGATAAGCGTAGCGAATCAGGCAATTTTAATGTTTAACTTCCCTGTTTAATCAGTTAct                                                                                                                                                                                                                                  >  1:89263/1‑275 (MQ=255)
                                                                               aGTAACAACAGCAATGCATATGAATAATTAAAGGTTATGCAACGGGCAAAGATTAAACTTCCCAGTTTAATGACAGGTTGCTGTACTTGCGTTTGTTGCTCTTGCCGGATGCAGCGTAAACGCCTTATCCGGCCTACAGTTCATTGCAAATCCATTCAATTTCAATGCACGATGTAGGCCTGATAAGCGTAGCGAATCAGGCAATTTTAATGTTTAACTTCCCTGTTTAATCAGTTACTGCTGGCGTAGCTGGAGATCCAGCGGCGTTTTGCTgg                                                                                                                                                                                              >  1:356453/1‑275 (MQ=255)
                                                                                         gCAATGCATATGAATAATTAAAGGTTATGCAACGGGCAAAGATTAAACTTCCCAGTTTAATGACAGGTTGCTGTACTTGCGTTTGTTGCTCTTGCCGGATGCAGCGTAAACGCCTTATCCGGCCTACAGTTCATTGCAAATCCATTCAATTTCAATGCACGATGTAGGCCTGATAAGCGTAGCGAATCAGGCAATTTTAATGTTTAACTTCCCTGTTTAATCAGTTACTGCTGGCg                                                                                                                                                                                                                           >  1:513975/1‑236 (MQ=255)
                                                                                             tGCATATGAATAATTAAAGGTTATTCAACGGGCAAAGATTAAACTTCCCAGTTTAATGACAGGTTGCTGTACTTGCGTTTGTTGCTCTTGCCGGATGCAGCGTAAACGCCTTATCCGGCCTACAGTTCATTGCAAATCCATTCAATTTCAATGCACGATGTAGGCCTGATAAGCGTAGCGAATCAGGCAATTTTAATGTTTAACTTCCCTGTTTAATCAGTTACTGCTGGCGTAGc                                                                                                                                                                                                                       <  2:513968/236‑1 (MQ=255)
                                                                                                atatGAATACTTAAAGGTTATGCAACGGGTAAATATTAAACTTCCCAGTTTAATGACAGGTTGCTGTACTTGCGTTTGTTGCTCTTGCCGGATGCAGCGTAAACGCCTTATCCGGCCTACAGTTCATTGCAAATCCATTCAATTTCAATGCACGATGTAGGCCTGATAAGCGTAGCGAATCAGGCAATTTTAATGTTTAACTTCCCTGTTTAATCAGTTACTGCTGGCGTAGCTGGAGATCCAGCGGCGTTTTGCTGGGTTCGCCGCCAATTtct                                                                                                                                                                             <  1:130623/275‑1 (MQ=255)
                                                                                                    gAATAATTAAAGGTTATGCAACGGGCAAAGATTAATCTTCCCAGTTTAATGACAGGTTGCTGTACTTGCGTTTGTTGCTCTTGCCGGATGCAGCGTAAACGCCTTATCCGGCCTACAGTTCATTGCAAATCCATTCAATTTCAATGCACGATGTAGGCCTGATAAGCGTAGCGAATCAGGCAATTTTAATGTTTAACTTCCCTGTTTAATCAGTTACTGCTGGCGTAGCTGGAGATCCAGCGGCGTTTTGCTGGGTTCGCCGCCAATTTCTcgcg                                                                                                                                                                         >  1:9551/1‑275 (MQ=255)
                                                                                                                    tGCAACGGGCAAAGATTAAACTTCCCAGTTTAATGACAGGTTGCTGTACTTGCGTTTGTTGCTCTTGCCGGATGCAGCGTAAACGCCTTATCCGGCCTACAGTTCATTGCAAATCCATTCAATTTCAATGCACGATGTAGGCCTGATAAGCGAAGCGAATCAGGCAATTTTAATGTTTAACTTCCCTGTTTAATCAGTTACTGCTGGCGTAGCTGGAGATCCAGCGGCGTTTTGCTGGGTTCGCCGCCAATTTCTCGCGCCAGTTTCGGCACCAg                                                                                                                                                         >  1:642427/1‑275 (MQ=255)
                                                                                                                           gtcAAAGTTTAAATTTCCCATTTTTATGACAGGTTTCTGTACTTGCGTTTTTTGCTCTTGCCGGATGCAGCGTAAACGCCTTATCCGGCCTACAGTTCATTGCAAATCCATTCAATTTCAATGCACGATGTAGGCCTGATAAGCGTAGCGAATCAGGCAATTTTAATGTTTAACTTCCCTGTTTAATCAGTTACTGCTGGCGTAGCTGGAGATCCAGCGGCGTTTTGCTGGGTTCGCCGCCAATTTCTCGCGCCAGTTTCGGCACCAGATATccc                                                                                                                                                  <  2:642418/273‑1 (MQ=255)
                                                                                                                           ggCAAAGATTAAACTTCCCAGTTTAATGACAGGTTGCTGTACTTGCGTTTGTTGCTCTTGCCGGATGCAGCGTAAACGCCTTATCCGGCCTACAGTTCATTGCAAATCCATTCAATTTCAATGCACGATGTAGGCCTGATAAGCGTAGCGAATCAGGCAATTTTAATGTTTAACTTCCCTGTTTAATCAGTTACTGCTGGCGTAg                                                                                                                                                                                                                        >  2:440304/1‑205 (MQ=255)
                                                                                                                               aaGATTAAACTTCCCAGTTTAATGACAGGTTGCTGTACTTGCGTTTGTTGCTCTTGCCGGATGCAGCGTAAACGCCTTATCCGGCCTACAGTTCATTGCAAATCCATTCAATTTCAATGCACGATGTAGGCCTGATAAGCGTAGCGAATCAGGCAATTTTAATGTTTAACTTCCCTGTTTAATCAGTTACTGCTGGCGTAGCTgg                                                                                                                                                                                                                    <  1:440311/205‑1 (MQ=255)
                                                                                                                               aaGATTAAACTTCCCAGTTTAATGACAGGTTGCTGTACTTGCGTTTGTTGCTCTTGCCGGATGCAGCGGAAACGCCTTATCCGGCCTACAGTTCATTGCAAATCCATTCAATTTCAATGCACGATGTAGGCCTGATAAGCGTAGCGAATCAGGCAATTTTAATGTTTAACTTCCCTGTTTAATCAGTTACTGCTGGCGTAGCTGGAGATCCAGCGGCGTTTTGCTGGGTTCGCCGCCAATTTCTCGCGCCAGTTTCGGCACCAGATATccc                                                                                                                                                  <  1:307164/271‑1 (MQ=255)
                                                                                                                                gtagtAAACGGCCCCGGTTAATGACAGGGTGCTTTACGTGCGTTTGTTGCTCTTGCCGGATGCAGAGTAAACGCCTTATCCGGCCGACAGTTCATTGCAAATCCATTCAAGTTCAATGCACGATGTAGGCCTGATAAGCGTAGCGAATCAGGCAATTTTAATGTTTAACTTCCCTGTTTAATCAGTTACTGCTGGCGTAGCTGGAGATCCAGCGGCGTTTTGCTGGGTTCGCCGCCAATTTCTCGCGCCAGTTTCGGCACCAGATATCCCGaca                                                                                                                                              <  2:174213/270‑1 (MQ=255)
                                                                                                                                       aCTTCCCAGTTTAATGACAGGTTGCTGTACTTGCGTTTGTTGCTCTTGCCGGATGCAGCGTAAACGCCTTATCCGGCCTACAGTTCATTGCAAATCCATTCAATTTCAATGCACGATGTAGGCCTGATAAGCGTAGCGAATCAGGCAAtt                                                                                                                                                                                                                                                                   >  2:67568/1‑150 (MQ=255)
                                                                                                                                           cccAGTTTAATGACAGGTTGCTGTACTTGCGTTTGTTGCTCTTGCCGGATGCAGCGTAAACGCCTTATCCGGCCTACAGTTCATTGCAAATCCATTCAATTTCAATGCACGATGTAGGCCTGATAAGCGTAGCGAATCAGGCAATtatac                                                                                                                                                                                                                                                               <  1:67569/150‑2 (MQ=255)
                                                                                                                                           cccAGTTTAATGACAGGTTGCTGTACTTGCGTTTGTTGCTCTTGCCGGATGCAGCGTAAACGCCTTATCCGGCCTACAGTTCATTGCAAATCCATTCAATTTCAATGCACGATGTAGGCCTGATAAGCGTAGCGAATCAGGCAATTTTAATGTTTAACTTCCCTGTTTAATCAGTTACTGCTGGCGTAGCTGGAGATCCAGCGGCGTTTTGCTGGGTTCGCCGCCAATTTCTCGCGCCAGTTTCGGCACCAGATATccc                                                                                                                                                  >  1:174216/1‑259 (MQ=255)
                                                                                                                                            ccAGTTTAATGACAGGTTGCTGTACTTGCGTTTGTTGCTCTTGCCGGATGCAGCGTAAACGCCTTATCCGGCCTACAGTTCATTGCAAATCCATTCAATTTCAATGCACGATGTAGGCCTGATAAGCGTAGCGAATCAGGCAATTTTAATGt                                                                                                                                                                                                                                                            >  1:29083/1‑152 (MQ=255)
                                                                                                                                             cggTTTAATGACAGGTTGCTGTACTTGCGTTTGTTGCTCTTGCCGGATGCAGCGTAAACGCCTTATCCGGCCTACAGTTCATTGCAAATCCATTCAATTTCAATGCACGATGTAGGCCTGATAAGCGTAGCGAATCAGGCAATTTTAATGTTTAACTTCCCTGTTTAATCAGTTACTGCTGGCGTAGCTGGAGATCCAGCGGCGTTTTGCTGGGTTCGCCGCCAATTTCTCGCGCCAGTTTCGGCACCAGATATCCCGACACCAGTGTCAGCAAc                                                                                                                                <  2:9551/273‑1 (MQ=255)
                                                                                                                                                tttAATGACAGGTTGCTGTACTTGCGTTTGTTGCTCTTGCCGGATGCAGCGTAAACGCCTTATCCGGCCTACAGTTCATTGCAAATCCATTCAATTTCAATGCACGATGTAGGCCTGATAAGCGTAGCGAATCAGGCAATTTTAATGTTTaa                                                                                                                                                                                                                                                        <  2:29083/152‑1 (MQ=255)
                                                                                                                                                     tGACAGGTTGCTGTACTTGCGTTTGTTGCTCTTGCCGGATGCAGCGTAAACGCCTTATCCGGCCTACAGTTCATTGCAAATCCATTCAATTTCAATGCACGATGTAGGCCTGATAAGCGTAGCGAATCAGGCAATTTTAATGTTTAACTTCCCTGTTTAATCAGTTActgct                                                                                                                                                                                                                               >  2:215872/1‑172 (MQ=255)
                                                                                                                                                      ggcAGGTTGTTGTACTTGCGTTTGTTGCTCTTGCCGGATGCAGCGTAAACGCCTTATCCGGCCTACAGTTCATTGCAAATCCNTTCAATTTNAATGCACGATGTAGGCCTGATAAGCGTAGCGAATCAGGCAATTTTAATGTTTAACTTCCCTGTTTAATCAGTTACTGCTGGCGTAGCTGGAGATCCAGCGGCGTTTTGCTGGGTTCGCCGCCAATTTCTCGCGCCAGTTTCGGCACCAGATATCCCGACACCAGTGTCAGCAACTCACGCATa                                                                                                                       <  1:474101/273‑1 (MQ=255)
                                                                                                                                                         aGGTTGCTGTACTTGCGTTTGTTGCTCTTGCCGGATGCAGCGTAAACGCCTTATCCGGCCTACAGTTCATTGCAAATCCATTCAATTTCAATGCACGATGTAGGCCTGATAAGCGTAGCGAATCAGGCAATTTTAATGTTTAACTTCCCTGTTTAATCAGTTACTGCTGGCg                                                                                                                                                                                                                           <  1:215876/172‑1 (MQ=255)
                                                                                                                                                          ggTTGCGGTACTTGCGTTTGTTGCTCTTGCCGGATGCAGCGTAAACGCCTTATCCGGCCTACAGTTCATTGCAAATCCATTCAATTTCAATGCACGATGTAGGCCTGATAAGCGTAGCGAATCAGGCAATTTTAATGTTTAACTTCCCTGTTTAATCAGTTACTGCTGGCGTAGCTGGAGATCNAGCGGCGTTTTGCTGGGTTCGCCGCCAATTTCTCGCGCCAGTTTCGGCACCAGATATCCCGACACCAGTGTCAGCAACTCAAGCAAAATCt                                                                                                                   >  1:513817/1‑275 (MQ=255)
                                                                                                                                                                            tgttgCTCTTGCCGGATGCAGCGTAAACGCCTTATCCGGCCTACAGTTCATTGCAAATCCATTCAATTTCAATGCACGATGTAGGCCTGATAAGCGTAGCGAATCAGGCAATTTTAATGTTTAACTTCCCTGTTTAATCAGTTACTGCTGGCGTAGCTGGAGATCCAGCGGCGTTTTGCTGGGTTCGCCGCCAATTTCTCGCGCCAGTTTCGGCACCAGat                                                                                                                                                       >  1:269984/1‑221 (MQ=255)
                                                                                                                                                                              ttgCTCTTGCCGGATGCAGCGTAAACGCCTTATCCGGCCTACAGTTCATTGCAAATCCATTCAATTTCAATGCACGATGTANGCCTGATANGCGTAGCGAATCTGGCAATTTTAATGTTTAACTTCCCTGTTTAATCAGTTACTGCGGGCGTAGCTGGAGATCCAGCGGCGTTTTGCTGGGTTCGCCGCCAATTTCTCGCGCCAGTTTCGGCACCAGATATCCCGACACCAGTGTCAGCCACTCACGCATAATCTGCCGTGCTTCGTCAtcact                                                                                                <  1:479433/274‑1 (MQ=255)
                                                                                                                                                                                gCTCTTGCCGGATGCAGCGTAAACGCCTTATCCGGCCTACAGTTCATTGCAAATCCATTCAATTTCAATGCACGATGTAGGCCTGATAAGCGTAGCGAATCAGGCAATTTTAATGTTTAACTTCCCTGTTTAATCAGTTACTGCTGGCGTAGCTGGAGATCCAGCGGCGTTTTGCTGGGTTCGCCGCCAATTTCTCGCGCCAGTTTCGGCACCAGATATcc                                                                                                                                                   <  2:269979/221‑1 (MQ=255)
                                                                                                                                                                                                gaataagagaaggaTCCGGCCTAATGTTCATTGCAAATCCATTCCATTTCAATGCACGATGTAGGCCTGATAAGCGTAGCGAATCAGGCAATTTTAATGTTTAACTTCCCTGTTTAATCAGTTACTGCTGGCGTAGCTGGAGATCCAGCGGCGTTTTGCTGGGTTCGCCGCCAATTTCTCGCGCCAGTTTCGGCACCAGATATCCCGACACCAGTGTCAGCAACTCACGCATAATCTGCCGTGCTTCGTCATCACTCACCATAAAATgcgccg                                                                               <  2:490683/260‑1 (MQ=255)
                                                                                                                                                                                                         ccTTATCCGGCCTACAGTTCATTGCAAATCCATTCAATTTCAATGCACGATGTAGGCCTGATAAGCGTAGCGAAt                                                                                                                                                                                                                                                                            >  2:22907/1‑75 (MQ=255)
                                                                                                                                                                                                           ttATCCGGCCTACAGTTCATTGCAAATCCATTCAATTTCAATGCACGATGTAGGCCTGATAAGCGTAGCGaa                                                                                                                                                                                                                                                                             >  2:34473/1‑72 (MQ=255)
                                                                                                                                                                                                             aTCCGGCCTACAGTTCATTGCAAATCCATTCAATTTCAATGCACGATGTAGGCCTGATAAGCGTAGCGAATCAgg                                                                                                                                                                                                                                                                        <  1:22907/75‑1 (MQ=255)
                                                                                                                                                                                                               ccGGCCTACAGTTCATTGCAAATCCATTCAATTTCAATGCACGATGTAGGCCTGATAAGCGTAGCGAATCAg                                                                                                                                                                                                                                                                         <  1:34473/72‑1 (MQ=255)
                                                                                                                                                                                                                    cTACAGTTCATTGCAAATCCATTCAATTTCAATGCACGATGTAGGCCTGATAAGCGTAGCGAATCAGGCAATTTTAATGTTTAACTTCCCTGTTTAATCAGTTACTGCTGGCGTAGCTGGAGATCCAGCGGCGTTTTGCTGGGTTCGCCGCCNATTTCTCGCGCCAGTTTCGGCACCAGATATCCCGACACCAGTGTCAGCAACTCACGCATAATCTGCCGTGCTTCGTCATCACTCACCATAAAATGCGCCGCGCCCTGTACTTTATCGaa                                                            >  1:607097/1‑271 (MQ=255)
                                                                                                                                                                                                                        aGTTCATTTCTAATTCATTCAATTTCAATGCACGATGTAGGCCTGATAAGCGTAGCGAATCAGGCAATTTTAATNTTTAACTTCCCTGTTTAATCAGTTACTGCTGGCGTAGCTGGAGATCCAGCGGCGTTTTGCTGGGTTCGCCGCCAATTTCTCGCGCCAGTTTCGGCACCAGATATCCCGACACCAGTGTCAGCAACTCACGCATAATCTGCCGTGCTTCGTCATCACTCACCATAAAATGCGCCGCGCCCTGTa                                                                      <  1:558430/258‑1 (MQ=255)
                                                                                                                                                                                                                        aGTTCATTGCAAATCAATTCAATATCAATACACGATGTAGGCCTGATAAGCGTAGCGAATCAGGCAATTTTAATGTTTAACTTCCCTGTTTAATCAGTTACTGCTGGCGTAGCTGGAGATCCAGCGGCGTTTTGCTGGGTTCGCCGCCAATTTCTCGCGCCAGTTTCGGCACCAGATATCCCGNCACCAGTGTCAGCAACTCACGCATAATCTGCCGTGCTTCGTCATCACTCACCATAAAATgc                                                                                   >  1:490690/1‑245 (MQ=255)
                                                                                                                                                                                                                          ttCATTGCAAATCCATTCAATTTCAATGCACGATGTAGGCCTGATAAGCGTAACGAATCAGGCAATTTTAATGttt                                                                                                                                                                                                                                                          >  2:485254/1‑76 (MQ=255)
                                                                                                                                                                                                                            cATTGCAAATCCATTCAATTTCAATGCACGATGTAGGCCTGATAAGCGTAGCGAATCAGGCAATTTTAATGTTTAACTTCACTGTTTAATCAGTTACTGCTGGCGTCTCTGGAGATCCAGCGGCCTTTTGCTGGGTTCGCCGCCAATTTCTCGCGCCAGTTTCGGCACCCGat                                                                                                                                                       >  2:497605/1‑173 (MQ=255)
                                                                                                                                                                                                                              ttGCAAATCCATTCAATTTCAATGCACGATGTAGGCCTGATAAGCGTCGCGAATCAGGCACTTTTAATGTTTAACt                                                                                                                                                                                                                                                      <  1:485261/76‑1 (MQ=255)
                                                                                                                                                                                                                               tGCAAATCCATTCAATTTCAATGCACGATGTAGGCCTGATAAGCGTAGCGAATCAGGCAATTTTAATGTTTAACTTCCCTGTTTAATCAGTTACTGCTGGCGTAGCTGGAGATCCAGCGGCGTTTTGCTGGGTTCGCCGCCAATTTCTCGCGCCAGTTTCGGCACCATATATCCCGACACCAGTGTCAGCAACTCACGCATAATCTTCCGTTCt                                                                                                           >  2:336582/1‑214 (MQ=255)
                                                                                                                                                                                                                               tGCAAATCCATTCAATTTCAATGCACGATGTAGGCCTGATAAGCGTAGCGAATCAGGCAATTTTAATGTTTAACTTCCCTGTTTAATCAGTTACTGCTGGCGTAGCTGGAGATCCAGCGGCGTTTTGCTGGGTTCGCCGCCAATTTCTCGCGCCAGTTTCGGCACCAGATATCCCGACACCAGNGTCAGCAANTCACGCATAATCTGCCGTGCTTCGTCATCACTCACAATAAAATGCGCCGCGCCCTGTCCTTTATCGAGc                                                           >  1:475049/1‑262 (MQ=255)
                                                                                                                                                                                                                                gCAAATCCATTCAATTTCAATGCACGATGTAGGCCTGATAAGCGTAGCGAGTCAGGCAATTTTAATGTTTAACTTCCCTGTTTAATCAGTTACTGCTGGCGTAGCTGGAGATCCAGCGGCGTTTTGCTGGGTTCGCCGCCAATGGCTCCCGCCAGTTTCGGCACCAGATATcc                                                                                                                                                   <  1:497612/173‑1 (MQ=255)
                                                                                                                                                                                                                                   aaTCCATTCACTTTCAATGCACGATGTAGGCCTGATAAGCGTAGCGAATCAGGCAATTTTAATGTTTAACTTCCCTGTTTAATCAGTTACTGCTGGCGTAGCTGGAGATCCAGCGGCGTTTTGCTGGGTTCgccgcctggtt                                                                                                                                                                               >  1:724414/1‑137 (MQ=255)
                                                                                                                                                                                                                                   aaTCCATTCAATTTCAATGCACGATGTAGGCCTGATAAGCGTAGCGAATCAGGCAATTTTAATGTTTAACTTCCCTGTTTAATCAGTTACTGCTGGCGTAGCTGGAGATCCAGCGGCGTTTTGCTGGGTTCGCCGCCAATTTCTCGCGCCAGTTTCGGCACCAGATATCCCGACACCAGTGTCAGCAACTCACGCATAATCTGCCGTGCTTCGt                                                                                                       <  1:336588/214‑1 (MQ=255)
                                                                                                                                                                                                                                    aTCCATTCAATTTCAATGCACGATGTAGGCCTGATaagcgaagcgca                                                                                                                                                                                                                                                                             >  1:572003/1‑45 (MQ=25)
                                                                                                                                                                                                                                    aTCCATTCAATTTCAATGCACGATGTAGGCCTGATAAGCGTAGCGAATc                                                                                                                                                                                                                                                                           >  1:392971/1‑49 (MQ=255)
                                                                                                                                                                                                                                    aTCCATTCAATTTCAATGCACGATGTAGGCCTGATAAGCGTAGCGAATc                                                                                                                                                                                                                                                                           >  2:197513/1‑49 (MQ=255)
                                                                                                                                                                                                                                    aTCCATTCAAATTCAATGCACGATGTAGGCCTGATAAGCGTAGCGAATc                                                                                                                                                                                                                                                                           >  1:723268/1‑49 (MQ=255)
                                                                                                                                                                                                                                        aTTCAATTTCAATGCACGATGTAGGCCTGATAAGCGTAGCGAATCAGGc                                                                                                                                                                                                                                                                       <  1:197516/49‑1 (MQ=255)
                                                                                                                                                                                                                                        aTTCAATTTCAATGCACGATGTAGGCCTGATAAGCGTAGCGAATCAGGc                                                                                                                                                                                                                                                                       <  2:723258/49‑1 (MQ=255)
                                                                                                                                                                                                                                        aTTCAATTTCAATGCACGATGTAGGCCTGATAAGCGTAGCGAATCAGGc                                                                                                                                                                                                                                                                       <  2:392965/49‑1 (MQ=255)
                                                                                                                                                                                                                                        aTTCAATTTCAATGCACGATGTAGGCCTGATAAGCGAAGCGCATCAg                                                                                                                                                                                                                                                                         <  2:571994/47‑1 (MQ=18)
                                                                                                                                                                                                                                         ttcaatttcaatGCACGATGTAGGCCTGATAAGCGTAGCGAATCAGGCAAtttt                                                                                                                                                                                                                                                                 >  1:412422/1‑54 (MQ=255)
                                                                                                                                                                                                                                         ttcaatttcaatGCACGATGTAGGCCTGATAAGCGTAGCGAATCAGGCAATTTTAATGTTTAACTTCCCTGTTTAATCAGTTACTGCTGGCGTAGCTGGAGATCCAGCGGCGTTTTGCTGGGTTCGCCGCCAATTTCTCGCGCCAGTTTCGGCAcc                                                                                                                                                           >  1:101362/1‑156 (MQ=255)
                                                                                                                                                                                                                                             atttcaatGCACGATGTAGGCCTGATAAGCGTAGCGAATCAGGCAATTTTAATg                                                                                                                                                                                                                                                             <  2:412415/54‑1 (MQ=255)
                                                                                                                                                                                                                                             atttcaatGCACGATGTAGGCCTGATAAGCGTAGCGAATCAGGCAATTTTAATGTTTAACTTCCCTGTTTAATCAGTTACTGCTGGCGTAGCTGGAGATCCAGCGGCGTTTTGCTGGGTTCGCCGCCAATTTCTCGCGCCAGTTTCGGCACCAGat                                                                                                                                                       <  2:101361/156‑1 (MQ=255)
                                                                                                                                                                                                                                                    tGCACGATGTAGGCCTGATAAGCGTAGCGAATCAGGCAATTTTAAt                                                                                                                                                                                                                                                              >  1:317666/1‑46 (MQ=255)
                                                                                                                                                                                                                                                        cGATGTAGGCCTGATAAGCGTAGCGAATCAGGCAATTTTAATGttt                                                                                                                                                                                                                                                          <  2:317660/46‑1 (MQ=255)
                                                                                                                                                                                                                                                          ttGTCGGCCTGATCAGAGTAGCGATTCCGGCAGTTTTAATGTTTAATTTCCCTGTTTAATCTGTTACTGCTGTCGTAGCTGGAGATCCAGCGGCGTTTTTCTGGGTTCGCCGCCAATTTCTCGCGCCAGTTTCGGCACCAGATATCCCGACACCAGTGTCAGCAACTCACGCATAATCTGCCGTGCTTCGTCATCACTCACCATAAAATGCGCCGCGCCCTGTACTTTATCGAGCACATGCAGGTAATACGGCATTACGCCGGCATCGAACAACg                   <  2:607088/274‑1 (MQ=255)
                                                                                                                                                                                                                                                           tGTAGGCCTGATAAGCGTAGCGAATCAGGCAATTTTAATg                                                                                                                                                                                                                                                             >  1:633812/1‑40 (MQ=255)
                                                                                                                                                                                                                                                               ggCCTGATAAGCGTAGCGAATCAGGCAATTTTAATGTTTa                                                                                                                                                                                                                                                         <  2:633803/40‑1 (MQ=255)
                                                                                                                                                                                                                                                               ggCCTGATAAGCGAAGCGCATCAGGCAATTTTAATGTTTAACTTCCCTGt                                                                                                                                                                                                                                               >  1:675815/1‑50 (MQ=21)
                                                                                                                                                                                                                                                                gCCTGATAAGCGTAGCGAATCAGGCAATTTTAATGt                                                                                                                                                                                                                                                            >  2:201706/1‑36 (MQ=255)
                                                                                                                                                                                                                                                                gCCTGATAAGCGTAGCGAATCAGGCAATTTTAATGTTTaa                                                                                                                                                                                                                                                        >  2:242557/1‑40 (MQ=255)
                                                                                                                                                                                                                                                                gCCTGATAAGCGTAGCGAATCAGGCAATTTTAATGTTTAACTTCCCTGTTTAATCAGTTACTGCTGGCGTAGCTGGAGATCCAGCGGCGTTTTGCTGGGTTCGCCGCCAATTTCTCGCGCCAGtt                                                                                                                                                                   >  1:483119/1‑125 (MQ=255)
                                                                                                                                                                                                                                                                gCCTGATAAGCGTAGCGAATCAGGCAATTTTAATGTTTAACTTCCCTGTTTAATCAGTTACTGCTGGCGTAGCTGGAGATCCAGCGGCGTTTTGCTGGGTTCGCCGCCAATTTCTCGCGCCAGTTTCGGCACCAGATATCCCGACACCAGTGTCAGCAACTCACGCATAATCTGCCGTGCTTCGTCATCACTCACCATAAAATgc                                                                                   >  1:668909/1‑205 (MQ=255)
                                                                                                                                                                                                                                                                gCCTGATAAGCGTAGCGAATCAGGCAATTTTAATGTTTAACTTCCCTGTTTAATCAGTTACTGCTGGCGTAGCTGGAGATCCAGCGGCGTTTTGCTGGGTTCGCCGCCAATTTCTCGCGCCAGTTTCGGCACCAGATATCCCGACACCAGTGTCAGCAACTCACGCATAATCTGCCGTGCTTCGTCATCACTCACCATAAAATGCGCCGCGCCCTGTACTTTATCGAGCACATGCAGGTAATACGGCATTACGCCGGCACCGAACAACGcccaac             >  2:204668/1‑270 (MQ=255)
                                                                                                                                                                                                                                                                gCCTGATAAGCGTAGCGAATCAGGCAATTTTAATGTTTAACTTCCCTGTTTAATCAGTTACTGCTGGCGTAGCTGGAGATCCAGCGGCGTTTTGCTGGGTTCGCCGCCAATTTCTCGCGCCAGTTTCGGCACCAGATATCCCGACACCAGTGTCAGCAACTCACGCATAATCTGCCGCGCTTCGTCATCACTCACCATAAAATGCGCCGCGCCCTGTACTTAATCGAGCACATGCAGGTAATACGGCATTACGCAGGCATCGAACAACGCATTaa             >  1:381066/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                gCCTGATAAGCGTAGCGAAGCAGGCAATTTTAATGTTTAACTTCCCTGt                                                                                                                                                                                                                                               >  1:566694/1‑49 (MQ=255)
                                                                                                                                                                                                                                                                gCCTGATAAGCGAAGCGCATCAGGCAATTTTAATGTTTAACTTCCCTGtt                                                                                                                                                                                                                                              >  2:538856/1‑50 (MQ=21)
                                                                                                                                                                                                                                                                gCCTGATAAGCGAAGCGCATCAGGCAATTTTAATGTTTAACTTCCCTGTTTAATCAGTTACTGCTGGCGTAGCTGGAGATCCAGCGGCGTTTTGCTGGGTTCGCCGCCAATTTCTCGCGCCAGTTTCGGCACCAGat                                                                                                                                                       >  1:175715/1‑137 (MQ=255)
                                                                                                                                                                                                                                                                gCCTGATAAGCGAAGCGCATCAGGCAATTTTAATGTTTAACTTCCCTGTTTAATCAGTTACTGCTGGCGTAGCTGGAGATCCAGCGGCGTTTTGCTGGGTTCGCCGCCAATTTCTCGCGCCAGTTTCGGCACCAGATATCCCGACACCAGTGTCAGc                                                                                                                                   >  2:215284/1‑157 (MQ=255)
                                                                                                                                                                                                                                                                gCCTGATAAGCGAAGCGAATCAGGCAATTTTAATGTTTAACTTCCCTGt                                                                                                                                                                                                                                               >  1:125850/1‑49 (MQ=255)
                                                                                                                                                                                                                                                                gCCTGATAAGCGAAGCGAATCAGGCAATTTTAATGTTTAACTTCCCTGTTTAATCAGTTACTGCTGGCGTAGCTGGAGATCCAGCGGCGTTTTGCTGGGTTCGCCGCCAATTTCTCGCGCCAGTTTCGGCACCACATATCCCGACACCAGTGTCAGCAACTCACGCCTAATCTGCCGTGCTTCGTCATCACTCCCCATAAAATGCGCCGCGCCCTGCACTTTATCGAGCACATGCAGGTAATACGGCATTACGCCGGCATCAAACAACGCATTa              >  1:259415/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                   tGATAAGCGTAGCGAATCAGGCAATTTGAATGTTTAACTTCCCTGTTTAATCAGTTACTGCTGGCGTAGCTGGAGATCCAGCGGCGTTTTGCTGGGTTCGCCGCCAATTTCTCGCGCCAGTTTCGGCACCAGATATCCCGACACCAGTGTCAGCAACTCACGCATAATCTGCCGTGCTTCGTCATCACTCACCATAAAATGCGCCGCGCCCTGTACTTTATCGAGCACATGCAGGTAATACGGCATTACGCCGGCATCGAACAACGCATTAc             >  1:410600/1‑272 (MQ=255)
                                                                                                                                                                                                                                                                   tGATAAGCGAAGCGCATCAGGCAATTTTAATGTTTAACTTCCCTGTTTaa                                                                                                                                                                                                                                           <  2:675806/50‑1 (MQ=25)
                                                                                                                                                                                                                                                                    gATAAGCGTAGCGAATCAGGCAATTTTAATGTTTaa                                                                                                                                                                                                                                                        <  1:201710/36‑1 (MQ=255)
                                                                                                                                                                                                                                                                    gATAAGCGTAGCGAATCAGGCAATTTTAATGTTTAACTTc                                                                                                                                                                                                                                                    <  1:242562/40‑1 (MQ=255)
                                                                                                                                                                                                                                                                    gATAAGCGTAGCGAATCAGGCAATTTTAATGTTTAACTTCCCTGTTTaa                                                                                                                                                                                                                                           <  2:566685/49‑1 (MQ=255)
                                                                                                                                                                                                                                                                    gATAAGCGTAGCGAATCAGGCAATTTTAATGTTTAACTTCCCTGTTTAATCAGTTACTGCTGGCGTAGCTGGAGATCCAGCGGCGTTTTGCTGGGTTCGCCGCCAATTTCTCGCGCCAGTTTCGGCACCAGATATCCCGACACCAGTGTCAGCAACTCACGCATAATCTGCCGTGCTTCGTCATCACTCACCATAAAATgcgccg                                                                               <  2:668900/205‑1 (MQ=255)
                                                                                                                                                                                                                                                                    gATAAGCGTAGCGAATCAAGCAATGTTAATGTTTAACGTCCCTGTTTAAGCAGTTACTGCTTGCGTAGCTGGAGATCCAACGGCGTTTTGCGGGGTTCGCCGCCAATTTCTCGCGCCAGTTTCgg                                                                                                                                                               <  2:483112/125‑1 (MQ=255)
                                                                                                                                                                                                                                                                    gATAAGCGAAGCGCATCAGGCTATTTTAATGTTTAACTTCCCTGTTTAAt                                                                                                                                                                                                                                          <  1:538864/50‑1 (MQ=255)
                                                                                                                                                                                                                                                                    gATAAGCGAAGCGCATCAGGCAATTTTAATGTTTAACTTCCCTGTTTAATCAGTTACTGCTGGCGTAGCTGGAGATCCAGCGGCGTTTTGCTGGGTTCGCCGCCAATTTCTCGCGCCAGTTTCGGCACCAGATATcc                                                                                                                                                   <  2:175712/137‑1 (MQ=255)
                                                                                                                                                                                                                                                                    gATAAGCGAAGCGCATCAGGCAATTTTAATGTTTAACTTCCCTGTTTAATCAGTTACTGCTGGCGTAGCTGGAGATCCAGCGGCGTTTTGCTGGGTTCGCCGCCAATTTCTCGCGCCAGTTTCGGCACCAGATATCCCGACACCAGTGTCAGCAACt                                                                                                                               <  1:215288/157‑1 (MQ=255)
                                                                                                                                                                                                                                                                    gATAAGCGAAGCGAATCAGGCAATTTTAATGTTTAACTTCCCTGTTTaa                                                                                                                                                                                                                                           <  2:125849/49‑1 (MQ=255)
                                                                                                                                                                                                                                                                     aTAAGCGTAGCGAATCAGGCAATTTTAATGTTTAACTTCCCTGTTTAATCAGTTACTGCTGGCGTAGCTGGAGATCCAGCGGCGTTTTGCTGGGTTCGCCGCCAATTTCTCGCGCCAGTTTCGGCACCAGATATCCCGACACCAGTGTCAGCAAc                                                                                                                                >  2:437082/1‑155 (MQ=255)
                                                                                                                                                                                                                                                                     aTAAGCGTAGCGAATCAGGCAATTTTAATGTTTAACTTCCCTGTTTAATCAGTTACTGCTGGCGTAGCTGGAGATCCAGCGGCGTTTTGCTGGGTTCGCCGCCAATTTCTCGCGCCAGTTTCGGCACCAGATATCCCGACACCAGTGTCAGCAACTc                                                                                                                              >  2:516209/1‑157 (MQ=255)
                                                                                                                                                                                                                                                                     aTAAGCGTAGCGAATCAGGCAATTTTAATGTTTAACTTCCCTGTTTAATCAGTTACTGCTGGCGTAGCTGGAGATCCAGCGGCGTTTTGCTGGGTTCGCCGCCAATTTCTCGCGCCAGTTTCGGCACCAGATATCCCGACACCAGTGTCAGCAACTCACGCATAATCTGCCGTGCTTCGtca                                                                                                     >  1:317172/1‑182 (MQ=255)
                                                                                                                                                                                                                                                                     aTAAGCGTAGCGAATCAGGCAATTTTAATGTTTAACTTCCCTGTTTAATCAGTTACTGCTGGCGTAGCTGGAGATCCAGCGGCGTTTTGCTGGGTTCGCCGCCAATTTCTCGCGCCAGTTTCGGCACCAGATATCCCGACACCAGTGTCAGCAACTCACGCATAATCTGCCGTGCTTCGTCATCACTCACCATAAAATGCGCCGCGCCCTGTACTTTATCGAGCACATGCAGGTAATACGGCATTACGCCGGCATCGAACAACGCATAACTCAgc        >  1:139504/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                         gCGTAGCGAATCAGGCAATTTTAATGTTTAACTTCCCTGTTTAATCAGTTACTGCTGGCGTAGCTGGAGATCCAGCGGCGTTTTGCTGGGTTCGCCGCCAATTTCTCGCGCCAGTTTCGGCACCAGATATCCCGACACCAGTGTCAGCAACTCACGc                                                                                                                          <  1:516216/157‑1 (MQ=255)
                                                                                                                                                                                                                                                                         gCGTAGCGAATCAGGCAATTTTAATGTTTAACTTCCCTGTTTAATCAGTTACTGCTGGCGTAGCTGGAGATCCAGCGGCGTTTTGCTGGGTTCGCCGCCAATTTCTCGCGCCAGTTTCGGCACCAGATATCCCGACACCAGTGTCAGCAACTCACGCATAATCTGCCGTGCTTCGTCAtcac                                                                                                 <  2:317166/182‑1 (MQ=255)
                                                                                                                                                                                                                                                                         gCGTAGCGAATCAGGCAATTTTAATGTTTAACTTCCCTGTTTAATCAGTTACTGCTGGCGTAGCTGGAGATCCAGCGGCGTTTTGCTGGGTTCGCCGCCAATTTCTCGCGCCAGTTTCGGAACCAGATATCCCGACACCAGTGTCAGCAACTCAc                                                                                                                            <  1:437089/155‑1 (MQ=255)
                                                                                                                                                                                                                                                                          cGTAGCGAATCAGGCAATTTTAATGTTTAACTTCCCTGTTTAATCAGTTACTGCTGGCGTAGCTGGAGATCCAGCGGCGTTTTGCTGGTTTCGCCGCCAATTTCTCGCGCCAGTTTCGTCACCCGATATCCCGACACCCGTGTCAGCCACTCACGCATAATCTGCCGTGCTTCTTCATCCCTCACCCTCAAATg                                                                                    >  2:425747/1‑194 (MQ=255)
                                                                                                                                                                                                                                                                          cGTAGCGAATCAGGCAATTTTAATGTTTAACTTCCCTGTTTAATCAGTTACTGCTGGCGTAGCTGGAGATCCAGCGGCGTTTTGCTGGGTTCGCCGCCAATTTCTCGCGCCAGTTTCGGCACCAGATATCCCGACACCAGTGTCAGCAACTCACGCATAATCTGCCGTGCTTCGtcatc                                                                                                   >  2:230363/1‑179 (MQ=255)
                                                                                                                                                                                                                                                                           gTAGCGAATCAGGCAATTTTAATGTTTAACTTCCCTGTTTAATCAGTTACTGCTGGCGTAGCTGGAGATCCAGCGGCGTTTTGCTGGGTTCGCCGCCAATTTCTCGCGCCAGTTTCGGCACCAGATATCCCGACACCAGTGTCAGCAACTCACGCATAATCTGCCGTGCTTCGTCATCACTCACCATAAAATGCGCCGCGCCCTGTACTTTATCGAGCACATGCAGGTAATACGGCCTTACGCCGGCATCGAAAAACGCACTACTCCGGTTTGcc  >  1:255258/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                 
GTTAAAATAACATAAATAGCCGTACTCACTCTATATAAACAGAAAAGAAAGGTTAATTATTGGTGTTAGCTATATAAAAAGTAACAACAGCAATGCATATGAATAATTAAAGGTTATGCAACGGGCAAAGATTAAACTTCCCAGTTTAATGACAGGTTGCTGTACTTGCGTTTGTTGCTCTTGCCGGATGCAGCGTAAACGCCTTATCCGGCCTACAGTTCATTGCAAATCCATTCAATTTCAATGCACGATGTAGGCCTGATAAGCGTAGCGAATCAGGCAATTTTAATGTTTAACTTCCCTGTTTAATCAGTTACTGCTGGCGTAGCTGGAGATCCAGCGGCGTTTTGCTGGGTTCGCCGCCAATTTCTCGCGCCAGTTTCGGCACCAGATATCCCGACACCAGTGTCAGCAACTCACGCATAATCTGCCGTGCTTCGTCATCACTCACCATAAAATGCGCCGCGCCCTGTACTTTATCGAGCACATGCAGGTAATACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCC  >  CP009273/4364133‑4364674

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