Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,221,214 C→A 100% 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,2140CA96.4% 593.8 / 11.8 203intergenic (+60/+73)rluF/yjbD23S rRNA pseudouridine(2604) synthase/DUF3811 family protein
Reads supporting (aligned to +/- strand):  ref base C (2/6);  new base A (88/107);  total (90/113)
Fisher's exact test for biased strand distribution p-value = 3.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 203 total aligned reads displayed.

GAAAAAACGGCGGAAAAAGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCAGGATCCGGTTTATAAGCCT  >  CP009273/4221071‑4221357
                                                                                                                                               |                                                                                                                                               
gAAAAAACGGCGGAAAAAGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTGAGGGa                                                                                                                                             >  2:1520240/1‑148 (MQ=255)
 aaaaaaCGGCGGAAAAAGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGaa                                                                                                                                            >  2:1754230/1‑148 (MQ=255)
  aaaaaCGGCGGAAAAAGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGaaa                                                                                                                                           >  1:957330/1‑148 (MQ=255)
   aaaaCGGCGGAAAAAGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGaaaa                                                                                                                                          <  2:657661/148‑1 (MQ=255)
   aaaaCGGCGGAAAAAGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGaaaa                                                                                                                                          <  2:2420507/148‑1 (MQ=255)
      aCGGCGGAAAAAGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAggg                                                                                                                                       <  2:3436908/148‑1 (MQ=255)
       cggcggAAAAAGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGt                                                                                                                                      <  1:3213744/148‑1 (MQ=255)
       cggcggAAAAAGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGt                                                                                                                                      >  1:729365/1‑148 (MQ=255)
        ggcggAAAAAGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGtt                                                                                                                                     <  1:2336779/148‑1 (MQ=255)
        ggcggAAAAAGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGtt                                                                                                                                     <  2:2702054/148‑1 (MQ=255)
        ggcggAAAAAGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGtt                                                                                                                                     <  1:1063542/148‑1 (MQ=255)
        ggcggAAAAAGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATGAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGtt                                                                                                                                     <  2:1883121/148‑1 (MQ=255)
         gcggAAAAAGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTa                                                                                                                                    <  2:1614673/148‑1 (MQ=255)
          ggaaaaaaGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGAAAAAACAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAg                                                                                                                                   <  1:1561826/145‑1 (MQ=255)
           ggAAAAAGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAgg                                                                                                                                  <  2:3493651/148‑1 (MQ=255)
            gAAAAAGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAggg                                                                                                                                 >  2:799804/1‑148 (MQ=255)
            gAAAAAGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAggg                                                                                                                                 <  1:103770/148‑1 (MQ=255)
             aaaaaGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGt                                                                                                                                >  1:939122/1‑148 (MQ=255)
             aaaaaGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGt                                                                                                                                >  1:2569952/1‑148 (MQ=255)
              aaaaGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTg                                                                                                                               >  1:1680473/1‑148 (MQ=255)
                      gTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTCGGGTGAGGGGaaa                                                                                                                       >  2:1212914/1‑148 (MQ=255)
                      gTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGaaa                                                                                                                       <  2:1905970/148‑1 (MQ=255)
                      gTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGaaa                                                                                                                       <  2:2632822/148‑1 (MQ=255)
                              gCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGtt                                                                                                               <  1:2694735/148‑1 (MQ=255)
                                gTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTca                                                                                                             <  2:71698/148‑1 (MQ=255)
                                    aaCGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACt                                                                                                         <  2:3137851/148‑1 (MQ=255)
                                    aaCGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACt                                                                                                         >  2:587629/1‑148 (MQ=255)
                                     aCGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTg                                                                                                        >  2:807711/1‑148 (MQ=255)
                                     aCGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTg                                                                                                        <  2:2629778/148‑1 (MQ=255)
                                        gTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACg                                                                                                     >  1:1659585/1‑148 (MQ=255)
                                         tAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGa                                                                                                    <  2:1074158/148‑1 (MQ=255)
                                          aaGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGACAATTGTGTTCACACTGACGAg                                                                                                   >  2:2340309/1‑148 (MQ=255)
                                           aGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGt                                                                                                  >  2:2318346/1‑148 (MQ=255)
                                            gCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCGTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTg                                                                                                 >  2:415276/1‑148 (MQ=255)
                                                tttACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACGGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGaaaa                                                                                             >  2:2611463/1‑148 (MQ=255)
                                                   aCCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAAtct                                                                                          <  2:2764661/148‑1 (MQ=255)
                                                    ccTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAAtctc                                                                                         <  2:2640688/148‑1 (MQ=255)
                                                      tCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAg                                                                                       >  1:2496784/1‑148 (MQ=255)
                                                         ccGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCg                                                                                    <  2:3246244/148‑1 (MQ=255)
                                                         ccGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCg                                                                                    <  2:3246208/148‑1 (MQ=255)
                                                           ggggCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGtt                                                                                  >  2:1560252/1‑148 (MQ=255)
                                                               cGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGAtt                                                                              <  2:788513/148‑1 (MQ=255)
                                                                  aaaaaGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATtgtg                                                                           <  2:832187/148‑1 (MQ=255)
                                                                    aaaGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGa                                                                         >  2:352173/1‑148 (MQ=255)
                                                                      agaagGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGAtt                                                                       >  2:2937099/1‑148 (MQ=255)
                                                                       gaagGGGCGCTGATTAACGCCTCAGACTGATGAACACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTa                                                                      <  1:1305400/148‑1 (MQ=255)
                                                                       gaagGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTa                                                                      <  1:1305864/148‑1 (MQ=255)
                                                                       gaagGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTa                                                                      >  1:3180086/1‑148 (MQ=255)
                                                                       gaagGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTa                                                                      <  1:2891406/148‑1 (MQ=255)
                                                                        aagGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAg                                                                     >  1:1554875/1‑148 (MQ=255)
                                                                            ggCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGcc                                                                 <  2:2137094/148‑1 (MQ=255)
                                                                            ggCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGcc                                                                 >  1:2633960/1‑148 (MQ=255)
                                                                                 tGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCcg                                                            <  1:3354177/148‑1 (MQ=255)
                                                                                    ttAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCcgcgc                                                         >  1:658830/1‑148 (MQ=255)
                                                                                     tAACGCCTCAGACTGATGAACACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCcgcgcg                                                        <  1:3398945/148‑1 (MQ=255)
                                                                                      aaCGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGt                                                       <  1:2178845/148‑1 (MQ=255)
                                                                                        cGCCTCAGCCTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCg                                                     <  2:2119958/148‑1 (MQ=255)
                                                                                        cGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCg                                                     >  1:1163501/1‑148 (MQ=255)
                                                                                            tCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGt                                                 <  2:211579/148‑1 (MQ=255)
                                                                                              aGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTAt                                               <  1:785533/148‑1 (MQ=255)
                                                                                               gACTGATGAAAACGGCGCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTAtt                                              >  2:989503/1‑148 (MQ=255)
                                                                                               gACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTAtt                                              <  1:1830623/148‑1 (MQ=255)
                                                                                               gACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTAtt                                              >  1:3491947/1‑148 (MQ=255)
                                                                                                aCTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTAttt                                             >  1:2238194/1‑148 (MQ=255)
                                                                                                 cTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTg                                            <  1:166465/148‑1 (MQ=255)
                                                                                                 cTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTg                                            <  1:2971752/148‑1 (MQ=255)
                                                                                                 cTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTg                                            <  2:2569952/148‑1 (MQ=255)
                                                                                                  tgatgaAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGc                                           <  2:20112/148‑1 (MQ=255)
                                                                                                  tgatgaAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGc                                           <  2:1773581/148‑1 (MQ=255)
                                                                                                    atgaAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGACAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTg                                         >  2:2414144/1‑148 (MQ=255)
                                                                                                    atgaAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTg                                         >  2:2413901/1‑148 (MQ=255)
                                                                                                    atgaAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGGTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTg                                         >  2:368032/1‑148 (MQ=255)
                                                                                                     tgaAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGa                                        >  1:597913/1‑148 (MQ=255)
                                                                                                         aaCGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAg                                    <  1:1009139/148‑1 (MQ=255)
                                                                                                          aCGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGc                                   >  1:2134529/1‑148 (MQ=255)
                                                                                                          aCGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGc                                   >  2:1290346/1‑148 (MQ=255)
                                                                                                           cGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCt                                  >  1:917399/1‑148 (MQ=255)
                                                                                                             gCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTaa                                <  1:1536523/148‑1 (MQ=255)
                                                                                                             gCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTaa                                <  1:306956/148‑1 (MQ=255)
                                                                                                             gCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATAGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTaa                                >  1:1111389/1‑148 (MQ=255)
                                                                                                              cccAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTaaa                               <  1:2777951/148‑1 (MQ=255)
                                                                                                               ccAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAAt                              >  1:402721/1‑148 (MQ=255)
                                                                                                                cAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATg                             >  2:1344661/1‑148 (MQ=255)
                                                                                                                cAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTAGTTGCTGACCAGCTAAATg                             >  1:1345288/1‑148 (MQ=255)
                                                                                                                   aaGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACg                          <  1:1115950/148‑1 (MQ=255)
                                                                                                                    aGTGCCCGGACAGTTCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGc                         <  1:876046/148‑1 (MQ=255)
                                                                                                                      tGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCtt                       <  2:729365/148‑1 (MQ=255)
                                                                                                                       gCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTc                      <  1:569880/148‑1 (MQ=255)
                                                                                                                        cccGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCa                     <  1:1754230/148‑1 (MQ=255)
                                                                                                                        cccGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCa                     <  2:418932/148‑1 (MQ=255)
                                                                                                                          cGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCAgg                   <  1:2870643/148‑1 (MQ=255)
                                                                                                                            gACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCAGGAt                 >  1:1348144/1‑148 (MQ=255)
                                                                                                                             aCAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCAGGATc                <  1:1474890/148‑1 (MQ=255)
                                                                                                                                  ccccTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCAGGATCCGGtt           <  2:1304251/148‑1 (MQ=255)
                                                                                                                                   cccTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCAGGATCCGGttt          <  1:2355801/148‑1 (MQ=255)
                                                                                                                                   cccTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCAGGATCCGGttt          <  2:1338923/148‑1 (MQ=255)
                                                                                                                                   cccTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCAGGATCCGGttt          >  2:3255396/1‑148 (MQ=255)
                                                                                                                                      tCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCAGGATCCGGTTtata       >  2:3019328/1‑148 (MQ=255)
                                                                                                                                      tCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCAGGATCCGGTTtata       >  1:3002417/1‑148 (MQ=255)
                                                                                                                                      tCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCAGGATCCGGTTtata       >  2:1800559/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                               
GAAAAAACGGCGGAAAAAGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCAGGATCCGGTTTATAAGCCT  >  CP009273/4221071‑4221357

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: