Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,221,214 C→A 59.1% intergenic (+60/+73) rluF → / ← yjbD 23S rRNA pseudouridine(2604) synthase/DUF3811 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,221,2140CA59.1% 31.5 / inf 255intergenic (+60/+73)rluF/yjbD23S rRNA pseudouridine(2604) synthase/DUF3811 family protein
Reads supporting (aligned to +/- strand):  ref base C (56/48);  new base A (76/74);  total (133/122)
Fisher's exact test for biased strand distribution p-value = 7.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.67e-01

Only 100 of 255 total aligned reads displayed.

ATGGAAAAAACGGCGGAAAAAGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCAGGATCCGGTTTATAAGCCTGC  >  CP009273/4221068‑4221359
                                                                                                                                                  |                                                                                                                                                 
  ggAAAAAACGGCGGAAAAAGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCggg                                                                                                                                                >  2:2076557/1‑148 (MQ=255)
  ggAAAAAACGGCGGAAAAAGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAggg                                                                                                                                                >  2:746017/1‑148 (MQ=255)
    aaaaaaCGGCGGAAAAAGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGaa                                                                                                                                              <  1:1098833/148‑1 (MQ=255)
    aaaaaaCGGCGGAAAAAGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGaa                                                                                                                                              <  2:669142/148‑1 (MQ=255)
      aaaaCGGCGGAAAAAGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGaaaa                                                                                                                                            <  2:48970/148‑1 (MQ=255)
       aaaCGGCGGAAAAAGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAg                                                                                                                                           <  2:3169695/148‑1 (MQ=255)
         aCGGCGGAAAAAGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAggg                                                                                                                                         >  2:1392045/1‑148 (MQ=255)
           ggcggAAAAAGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGtt                                                                                                                                       <  1:1698616/148‑1 (MQ=255)
                aaaaaGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGt                                                                                                                                  >  2:2323636/1‑148 (MQ=255)
                 aaaaGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTg                                                                                                                                 >  2:2221214/1‑148 (MQ=255)
                  aaaGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGa                                                                                                                                >  2:507635/1‑148 (MQ=255)
                  aaaGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGa                                                                                                                                >  1:692251/1‑148 (MQ=255)
                  aaaGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGa                                                                                                                                <  2:1252362/148‑1 (MQ=255)
                   aaGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAg                                                                                                                               >  1:852832/1‑148 (MQ=255)
                    aGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAgg                                                                                                                              >  1:3152934/1‑148 (MQ=255)
                    aGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAgg                                                                                                                              <  2:3013292/148‑1 (MQ=255)
                      gCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAgggg                                                                                                                            >  1:3397666/1‑148 (MQ=255)
                          tCGCCCGGCGTCCACCGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGaaaa                                                                                                                        <  1:2122616/148‑1 (MQ=255)
                           cGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAAt                                                                                                                       <  2:2543194/148‑1 (MQ=255)
                           cGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAAt                                                                                                                       <  2:3518400/148‑1 (MQ=255)
                           cGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAAt                                                                                                                       <  2:1756726/148‑1 (MQ=255)
                             cccGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATtg                                                                                                                     <  1:1392386/148‑1 (MQ=255)
                              ccGGCGTCCAAGGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTTATTAACGCCTCAGACTGATGAAAACGGCCCAACAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATtgt                                                                                                                    <  1:102361/148‑1 (MQ=255)
                              ccGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATtgt                                                                                                                    <  1:150139/148‑1 (MQ=255)
                              ccGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATtgt                                                                                                                    >  2:2825258/1‑148 (MQ=255)
                                ggCGTCCAACGGTAAGCGTTGTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTgtgt                                                                                                                  <  2:2372994/148‑1 (MQ=255)
                                  cGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTc                                                                                                                >  1:2995462/1‑148 (MQ=255)
                                     ccAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTcaca                                                                                                             <  2:1731430/148‑1 (MQ=255)
                                      cAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCacac                                                                                                            >  1:3530780/1‑148 (MQ=255)
                                       aaCGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACt                                                                                                           >  1:2353994/1‑148 (MQ=255)
                                            tAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGa                                                                                                      >  2:862268/1‑148 (MQ=255)
                                             aaGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAg                                                                                                     <  1:510255/148‑1 (MQ=255)
                                             aaGCGATTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAg                                                                                                     <  2:2667552/148‑1 (MQ=255)
                                              aGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGt                                                                                                    <  2:199163/148‑1 (MQ=255)
                                               gCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTg                                                                                                   >  1:72778/1‑148 (MQ=255)
                                                   tttACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGaaaa                                                                                               >  1:3196673/1‑148 (MQ=255)
                                                    ttACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAAt                                                                                              <  2:1988952/148‑1 (MQ=255)
                                                      aCCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAAtct                                                                                            >  1:3274277/1‑148 (MQ=255)
                                                       ccTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAAtctc                                                                                           <  1:547023/148‑1 (MQ=255)
                                                        cTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCa                                                                                          >  2:2473546/1‑148 (MQ=255)
                                                          cGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAgg                                                                                        <  2:2162547/148‑1 (MQ=255)
                                                          cGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAgg                                                                                        >  1:2029545/1‑148 (MQ=255)
                                                               gggCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTg                                                                                   >  1:2576037/1‑148 (MQ=255)
                                                                 gCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGAt                                                                                 <  2:1249790/148‑1 (MQ=255)
                                                                  cGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGAtt                                                                                >  1:1143702/1‑148 (MQ=255)
                                                                   gTAAAAAGATGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATtg                                                                               <  2:3293301/148‑1 (MQ=255)
                                                                   gTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATtg                                                                               <  2:1523470/148‑1 (MQ=255)
                                                                     aaaaaGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATtgtg                                                                             >  1:1042661/1‑148 (MQ=255)
                                                                     aaaaaGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATtgtg                                                                             >  2:3122517/1‑148 (MQ=255)
                                                                      aaaaGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGATGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTgg                                                                            >  1:654563/1‑148 (MQ=255)
                                                                       aaaGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGa                                                                           <  1:1763171/148‑1 (MQ=255)
                                                                         agaagGGGCGCTGATTAGCGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGAtt                                                                         >  1:3354495/1‑148 (MQ=255)
                                                                          gaagGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTa                                                                        >  2:152837/1‑148 (MQ=255)
                                                                              gggCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAgcgc                                                                    <  2:107666/148‑1 (MQ=255)
                                                                               ggCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGcc                                                                   >  2:2773766/1‑148 (MQ=255)
                                                                                  gcTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGc                                                                <  1:1392045/148‑1 (MQ=255)
                                                                                  gcTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGc                                                                <  1:3257130/148‑1 (MQ=255)
                                                                                   cTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGcc                                                               >  1:421555/1‑148 (MQ=255)
                                                                                     gATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCcgc                                                             >  2:1947700/1‑148 (MQ=255)
                                                                                     gATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCcgc                                                             >  1:2316047/1‑148 (MQ=255)
                                                                                     gATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCcgc                                                             >  2:115230/1‑148 (MQ=255)
                                                                                      aTTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCcgcg                                                            >  1:2256144/1‑148 (MQ=255)
                                                                                       ttAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGTGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCcgcgc                                                           >  1:2639318/1‑148 (MQ=255)
                                                                                         aaCGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGACAAGGGTTAGGGTGAGGGGAACATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGt                                                         >  2:3302505/1‑148 (MQ=255)
                                                                                               tCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGt                                                   <  1:1522451/148‑1 (MQ=255)
                                                                                                 aGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTAt                                                 <  1:380164/148‑1 (MQ=255)
                                                                                                  gACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTAtt                                                <  2:410423/148‑1 (MQ=255)
                                                                                                  gACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTAtt                                                >  1:2494093/1‑148 (MQ=255)
                                                                                                    cTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTg                                              <  2:259922/148‑1 (MQ=255)
                                                                                                       atgaAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTg                                           >  1:771895/1‑148 (MQ=255)
                                                                                                       atgaAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTg                                           <  2:617387/148‑1 (MQ=255)
                                                                                                       atgaAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTg                                           <  2:430058/148‑1 (MQ=255)
                                                                                                        tgaAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGa                                          >  2:1658355/1‑148 (MQ=255)
                                                                                                             aCGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGc                                     >  1:458854/1‑148 (MQ=255)
                                                                                                             aCGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGc                                     >  1:3035733/1‑148 (MQ=255)
                                                                                                              cGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCt                                    >  2:424652/1‑148 (MQ=255)
                                                                                                                 cccAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTaaa                                 >  1:3162938/1‑148 (MQ=255)
                                                                                                                  ccAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAAt                                <  1:746017/148‑1 (MQ=255)
                                                                                                                   cAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATg                               >  2:3230450/1‑148 (MQ=255)
                                                                                                                    aaaaGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGa                              <  2:2191168/148‑1 (MQ=255)
                                                                                                                    aaaaGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGa                              >  2:3429938/1‑148 (MQ=255)
                                                                                                                     aaaGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGAc                             >  2:1846060/1‑148 (MQ=255)
                                                                                                                         tGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCtt                         >  2:2237381/1‑148 (MQ=255)
                                                                                                                         tGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCtt                         >  2:469359/1‑148 (MQ=255)
                                                                                                                          gCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTc                        >  1:2186733/1‑148 (MQ=255)
                                                                                                                          gCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTc                        <  1:932674/148‑1 (MQ=255)
                                                                                                                              ggACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCAGGa                    <  1:469403/148‑1 (MQ=255)
                                                                                                                                aCAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCAGGATc                  <  1:3220284/148‑1 (MQ=255)
                                                                                                                                  aGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCAGGATCCg                >  1:3005097/1‑148 (MQ=255)
                                                                                                                                  aGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCAGGATCCg                <  2:3238033/148‑1 (MQ=255)
                                                                                                                                   gTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCAGGATCCgg               >  2:188318/1‑148 (MQ=255)
                                                                                                                                    tCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCAGGATCCGGt              >  2:2476480/1‑148 (MQ=255)
                                                                                                                                     ccccTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCAGGATCCGGtt             <  2:1590612/148‑1 (MQ=255)
                                                                                                                                     ccccTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTACCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCAGGATCCGGtt             <  2:3359430/148‑1 (MQ=255)
                                                                                                                                        cTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCAGGATCCGGTTtat          <  2:690680/148‑1 (MQ=255)
                                                                                                                                             cccTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCAGGATCCGGTTTATAAGcc     >  2:965704/1‑148 (MQ=255)
                                                                                                                                              ccTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCAGGATCCGGTTTATAAGCCt    >  2:2073761/1‑148 (MQ=255)
                                                                                                                                               cTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCAGGATCCGGTTTATAAGCCTg   >  1:2906420/1‑148 (MQ=255)
                                                                                                                                               cTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCAGGATCCGGTTTATAAGCCTg   >  1:1976056/1‑148 (MQ=255)
                                                                                                                                               cTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCAGGATCCGGTTTATAAGCCTg   >  2:858861/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
ATGGAAAAAACGGCGGAAAAAGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCAGGATCCGGTTTATAAGCCTGC  >  CP009273/4221068‑4221359

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: