Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,221,214 C→A 26.4% 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,2140CA26.4% 507.3 / 133.1 222intergenic (+60/+73)rluF/yjbD23S rRNA pseudouridine(2604) synthase/DUF3811 family protein
Reads supporting (aligned to +/- strand):  ref base C (66/97);  new base A (27/32);  total (93/129)
Fisher's exact test for biased strand distribution p-value = 5.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 223 total aligned reads displayed.

GATGGAAAAAACGGCGGAAAAAGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCAGGATCCGGTTTATAAGCCTGC  >  CP009273/4221067‑4221359
                                                                                                                                                   |                                                                                                                                                 
  tGGAAAAAACGGCGGAAAAAGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAgg                                                                                                                                                 <  2:1687401/148‑1 (MQ=255)
       aaaaCGGCGGAAAAAGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGaaaa                                                                                                                                            <  2:871786/148‑1 (MQ=255)
       aaaaCGGCGGAAAAAGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGaaaa                                                                                                                                            >  2:2399281/1‑148 (MQ=255)
        aaaCGGCGGAAAAAGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAg                                                                                                                                           <  2:2985441/148‑1 (MQ=255)
         accggcggAAAAAGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAgg                                                                                                                                          <  2:459094/146‑1 (MQ=255)
         aaCGGCGGAAAAAGGCGGTCGCCCGGCGTCCAACGTTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGTAAAAgg                                                                                                                                          <  2:568534/148‑1 (MQ=255)
          aCGGCGGAAAAAGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAggg                                                                                                                                         <  1:2868437/148‑1 (MQ=255)
           cggcggAAAAAGGCGGTCGCCCGGCCTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGt                                                                                                                                        >  1:3288486/1‑148 (MQ=255)
                 aaaaaGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGt                                                                                                                                  <  1:1755116/148‑1 (MQ=255)
                    aaGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAg                                                                                                                               <  1:1311818/148‑1 (MQ=255)
                    aaGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAg                                                                                                                               <  1:2959873/148‑1 (MQ=255)
                    aaGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGGTAGGGTGAg                                                                                                                               >  2:1356934/1‑148 (MQ=255)
                      ggcgggcgCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAggg                                                                                                                             <  1:2260518/148‑1 (MQ=255)
                       gCGGTCGCCCGGCGTCCACCGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAgggg                                                                                                                            <  1:721565/148‑1 (MQ=255)
                       gCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAgggg                                                                                                                            >  2:893730/1‑148 (MQ=255)
                          gTCGCCCGGCGTCCACCGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGaaa                                                                                                                         <  2:1665122/148‑1 (MQ=255)
                          gTCGCCCGGCGTCCAACGGTTAGCGTTTGACCTCGCCGGGGCGTTAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGaaa                                                                                                                         <  2:713646/148‑1 (MQ=255)
                           tCGCCCGGCTTCCAACGGTAAGCTTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGaaaa                                                                                                                        <  1:2515063/148‑1 (MQ=255)
                           tCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAACACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGaaaa                                                                                                                        <  2:3053178/148‑1 (MQ=255)
                           tCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGaaaa                                                                                                                        <  1:1340872/148‑1 (MQ=255)
                                cGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATtgtg                                                                                                                   <  2:1324956/148‑1 (MQ=255)
                                 ggCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTgtgt                                                                                                                  >  1:766502/1‑148 (MQ=255)
                                   cGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAACTTGTGTTc                                                                                                                >  2:3182623/1‑148 (MQ=255)
                                         ccGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTg                                                                                                          <  1:1459562/147‑1 (MQ=255)
                                         aCGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCGTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTg                                                                                                          >  1:2125917/1‑148 (MQ=255)
                                          cGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGa                                                                                                         <  2:483921/148‑1 (MQ=255)
                                              aaGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAg                                                                                                     <  1:2809004/148‑1 (MQ=255)
                                              aaGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAg                                                                                                     >  2:33648/1‑148 (MQ=255)
                                              aaGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAg                                                                                                     >  1:566665/1‑148 (MQ=255)
                                               aGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGt                                                                                                    >  1:1080716/1‑148 (MQ=255)
                                                 cGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGa                                                                                                  >  2:1866687/1‑148 (MQ=255)
                                                    tttACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGaaaa                                                                                               >  1:1612453/1‑148 (MQ=255)
                                                     ttACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAAt                                                                                              >  1:3037657/1‑148 (MQ=255)
                                                     ttACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAAt                                                                                              <  2:1591733/148‑1 (MQ=255)
                                                      tACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAAtc                                                                                             >  1:63282/1‑148 (MQ=255)
                                                       aCCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAAtct                                                                                            >  1:856389/1‑148 (MQ=255)
                                                       aCCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAAtct                                                                                            >  2:1178101/1‑148 (MQ=255)
                                                       aCCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAAtct                                                                                            <  2:2605164/148‑1 (MQ=255)
                                                        ccTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAAtctc                                                                                           >  2:3120968/1‑148 (MQ=255)
                                                        ccTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAAtctc                                                                                           >  1:2996502/1‑148 (MQ=255)
                                                         cTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCa                                                                                          <  2:1143896/148‑1 (MQ=255)
                                                             ccGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCg                                                                                      <  2:2492/148‑1 (MQ=255)
                                                             ccGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCg                                                                                      <  2:3257073/148‑1 (MQ=255)
                                                              cGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGt                                                                                     <  1:1916152/148‑1 (MQ=255)
                                                                gggCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTg                                                                                   >  2:727563/1‑148 (MQ=255)
                                                                 ggCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGa                                                                                  >  1:2112165/1‑148 (MQ=255)
                                                                  gCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGAt                                                                                 <  2:1903933/148‑1 (MQ=255)
                                                                  gCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGAt                                                                                 <  2:1618210/148‑1 (MQ=255)
                                                                  gCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGAt                                                                                 <  2:2717993/148‑1 (MQ=255)
                                                                      aaaaaGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATtgtg                                                                             >  1:1136387/1‑148 (MQ=255)
                                                                      aaaaaGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATtgtg                                                                             >  1:1136378/1‑148 (MQ=255)
                                                                      aaaaaGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAACGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATATCAGGCGTTGATtgtg                                                                             >  2:2009699/1‑148 (MQ=255)
                                                                        aaaGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGa                                                                           >  2:217652/1‑148 (MQ=255)
                                                                        aaaGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGa                                                                           <  2:937648/148‑1 (MQ=255)
                                                                              ggggCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAgcg                                                                     >  2:2707391/1‑148 (MQ=255)
                                                                               gggCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAgcgc                                                                    <  2:1989431/148‑1 (MQ=255)
                                                                                       aTTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAATGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCcgcg                                                            >  1:2539218/1‑148 (MQ=255)
                                                                                       aTTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCcgcg                                                            >  1:800885/1‑148 (MQ=255)
                                                                                       aTTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCcgcg                                                            <  1:3203655/148‑1 (MQ=255)
                                                                                       aTTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGACAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGAGTGTGGATTAGCGCCTGCcgcg                                                            >  2:2412612/1‑148 (MQ=255)
                                                                                       aTTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCcgcg                                                            >  1:2412306/1‑148 (MQ=255)
                                                                                          aaCGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGt                                                         <  2:93252/148‑1 (MQ=255)
                                                                                           aCGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTc                                                        >  1:1022460/1‑148 (MQ=255)
                                                                                           aCGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTc                                                        >  1:1231693/1‑148 (MQ=255)
                                                                                            cGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCg                                                       <  1:182386/148‑1 (MQ=255)
                                                                                                tCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGt                                                   <  2:297688/148‑1 (MQ=255)
                                                                                                  aGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTAt                                                 >  2:2785721/1‑148 (MQ=255)
                                                                                                  aGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTAt                                                 >  2:1250828/1‑148 (MQ=255)
                                                                                                   gACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTAtt                                                >  1:1655333/1‑148 (MQ=255)
                                                                                                     cTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTg                                              <  1:2869432/148‑1 (MQ=255)
                                                                                                     cTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTg                                              <  1:1172975/148‑1 (MQ=255)
                                                                                                      tgatgaAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGc                                             >  2:2872685/1‑148 (MQ=255)
                                                                                                        atgaAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTg                                           <  1:892623/148‑1 (MQ=255)
                                                                                                         tgaAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTGAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGa                                          >  1:1804704/1‑148 (MQ=255)
                                                                                                             aaCGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAg                                      <  2:620857/148‑1 (MQ=255)
                                                                                                             aaCGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAg                                      >  1:2639064/1‑148 (MQ=255)
                                                                                                                  cccAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTaaa                                 <  2:841503/148‑1 (MQ=255)
                                                                                                                  cccAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTaaa                                 >  1:1527967/1‑148 (MQ=255)
                                                                                                                   ccAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAAt                                <  1:2399281/148‑1 (MQ=255)
                                                                                                                    cAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATg                               <  2:2218138/148‑1 (MQ=255)
                                                                                                                    cAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATg                               <  2:1882542/148‑1 (MQ=255)
                                                                                                                     aaaaGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGa                              <  1:1722282/148‑1 (MQ=255)
                                                                                                                       aaGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACg                            <  2:2960201/148‑1 (MQ=255)
                                                                                                                       aaGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACg                            <  2:827275/148‑1 (MQ=255)
                                                                                                                       aaGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGACTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACg                            <  2:3220079/148‑1 (MQ=255)
                                                                                                                         gTGCCCGGACAGTCCCCTCGCCCCTTCGGGAACAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCt                          >  2:2365753/1‑148 (MQ=255)
                                                                                                                          tGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCtt                         <  1:2585133/148‑1 (MQ=255)
                                                                                                                           gCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTc                        <  1:2179258/148‑1 (MQ=255)
                                                                                                                            cccGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCa                       <  2:1805321/148‑1 (MQ=255)
                                                                                                                             ccGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCAg                      >  1:2236361/1‑148 (MQ=255)
                                                                                                                               ggACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCAggc                    >  2:1337245/1‑147 (MQ=255)
                                                                                                                                        ccTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCAGGATCCGGTTta           >  1:840870/1‑148 (MQ=255)
                                                                                                                                         cTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCAGGATCCGGTTtat          <  1:2707391/148‑1 (MQ=255)
                                                                                                                                          tCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCAGGATCCGGTTtata         >  1:1249597/1‑148 (MQ=255)
                                                                                                                                          tCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCAGGATCCGGTTtata         >  1:550308/1‑148 (MQ=255)
                                                                                                                                            gCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCAGGATCCGGTTTATAAg       <  1:994920/148‑1 (MQ=255)
                                                                                                                                             ccccTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCAGGATCCGGTTTATAAGc      >  1:2030335/1‑148 (MQ=255)
                                                                                                                                             ccccTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCAGGATCCGGTTTATAAGc      <  2:2284749/148‑1 (MQ=255)
                                                                                                                                              cccTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCAGGATCCGGTTTATAAGcc     <  1:689545/148‑1 (MQ=255)
                                                                                                                                                 ttCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCAGGATCCGGTTTATAAGCCTGc  >  1:1504644/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GATGGAAAAAACGGCGGAAAAAGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCAGGATCCGGTTTATAAGCCTGC  >  CP009273/4221067‑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: