Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,220,948 G→T 7.4% D223Y (GAT→TAT)  rluF → 23S rRNA pseudouridine(2604) synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,220,9480GT 7.4% 656.4 / 9.0 217D223Y (GAT→TAT) rluF23S rRNA pseudouridine(2604) synthase
Reads supporting (aligned to +/- strand):  ref base G (89/104);  new base T (4/20);  total (93/124)
Fisher's exact test for biased strand distribution p-value = 7.64e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 221 total aligned reads displayed.

CCGTTTGTCTTCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTCGTTAAGATGGAAAAAACGGCGGAAAAAGGCGG  >  CP009273/4220801‑4221093
                                                                                                                                                   |                                                                                                                                                 
  gTTTGTCTTCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGAt                                                                                                                                                 <  1:2511187/148‑1 (MQ=255)
  gTTTGTCTTCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGAt                                                                                                                                                 >  1:2442874/1‑148 (MQ=255)
  gTTTGTCTTCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGAt                                                                                                                                                 <  2:3509590/148‑1 (MQ=255)
    ttGTCTTCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGACTGGCGCGAttt                                                                                                                                               >  1:3490468/1‑148 (MQ=255)
    ttGTCTTCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGAttt                                                                                                                                               >  2:2797620/1‑148 (MQ=255)
     tGTCTTCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTa                                                                                                                                              <  2:1131648/148‑1 (MQ=255)
      gTCTTCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTaa                                                                                                                                             <  2:1706107/148‑1 (MQ=255)
        cTTCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAAcc                                                                                                                                           >  1:2081845/1‑148 (MQ=255)
              cATTACCCTGGTGCTGGGGCTGAACCGTCAGTTCCGGCGCATGTGCGAGAATTTGGGCTTTTAAGTGAAAAAGCTGTAACGCACGCGCATCATGAACGTTAGATTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGAt                                                                                                                                     <  2:2776591/148‑1 (MQ=255)
              cATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGAt                                                                                                                                     >  2:2868028/1‑148 (MQ=255)
                ttACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGa                                                                                                                                   >  1:32440/1‑148 (MQ=255)
                ttACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGa                                                                                                                                   >  2:3107385/1‑148 (MQ=255)
                 gaCCCTGGTGCTGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAg                                                                                                                                  <  2:1785531/147‑1 (MQ=255)
                  aCCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGt                                                                                                                                 <  2:1108896/148‑1 (MQ=255)
                    ccTGGTGCAGGGGGTGAACCGTAAGATCGGGCGCATGTGCGAGCATTTCGGCTTTGAAGTGAAAAAGCTGGAACGAACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTa                                                                                                                               <  2:328017/148‑1 (MQ=255)
                    ccTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTa                                                                                                                               >  2:764461/1‑148 (MQ=255)
                     cTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTaa                                                                                                                              <  2:2480388/148‑1 (MQ=255)
                     cTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTaa                                                                                                                              <  1:1751661/148‑1 (MQ=255)
                      tGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAAt                                                                                                                             >  1:3280035/1‑148 (MQ=255)
                      tGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAAt                                                                                                                             <  1:1123707/148‑1 (MQ=255)
                      tGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAAt                                                                                                                             >  2:815804/1‑148 (MQ=255)
                          gCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGAc                                                                                                                         <  2:1532307/148‑1 (MQ=255)
                           cAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGAcc                                                                                                                        <  2:1629402/148‑1 (MQ=255)
                             ggggCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACctc                                                                                                                      >  1:26207/1‑148 (MQ=255)
                              gggCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACctct                                                                                                                     <  2:1271318/148‑1 (MQ=255)
                              gggCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACctct                                                                                                                     >  1:918530/1‑148 (MQ=255)
                              gggCTGAACCGTCAGATCCGGCGAATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAAGGTTAGCTTAAGCGGCATTCCGCTGGGGGATTGGCGCGATTTAACCGACGATGAGTTAATCGACctct                                                                                                                     <  1:745296/148‑1 (MQ=255)
                                 cTGACCCGTCCGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAACGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTa                                                                                                                  <  2:2963673/148‑1 (MQ=255)
                                   gAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGATATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAg                                                                                                                <  1:3041822/148‑1 (MQ=255)
                                      ccGTCAGATCCGGCGCATGTGCGAGCATTTCGGATATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTc                                                                                                             <  1:520119/148‑1 (MQ=255)
                                       cGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCa                                                                                                            <  2:656991/148‑1 (MQ=255)
                                          cAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACTCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACATCTTTAAGCTCATTg                                                                                                         >  1:1815090/1‑148 (MQ=255)
                                          cAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGACTGGCGCTATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTg                                                                                                         >  2:3252440/1‑148 (MQ=255)
                                          cAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTg                                                                                                         >  2:2886076/1‑148 (MQ=255)
                                          cAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTg                                                                                                         >  1:82030/1‑148 (MQ=255)
                                          cAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTg                                                                                                         >  1:934176/1‑148 (MQ=255)
                                                cGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAAtt                                                                                                   <  2:2256020/148‑1 (MQ=255)
                                                  gcgcATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTcc                                                                                                 >  1:356867/1‑148 (MQ=255)
                                                   cgcATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCct                                                                                                <  1:2219892/148‑1 (MQ=255)
                                                     cATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCctct                                                                                              <  1:185123/148‑1 (MQ=255)
                                                      aTGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCtt                                                                                             <  1:1360216/148‑1 (MQ=255)
                                                       tgtgCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTc                                                                                            >  1:2577878/1‑148 (MQ=255)
                                                         tgCGAGCATTTCGGCTATGAAGTTAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCg                                                                                          <  1:352930/148‑1 (MQ=255)
                                                              gCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTa                                                                                     <  2:662471/148‑1 (MQ=255)
                                                                  ttCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAAc                                                                                 >  1:2030967/1‑148 (MQ=255)
                                                                      gCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATTAACGTTCGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTaa                                                                             <  1:1453404/148‑1 (MQ=255)
                                                                      gCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTaa                                                                             >  1:2460337/1‑148 (MQ=255)
                                                                      gCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTaa                                                                             <  2:1322119/148‑1 (MQ=255)
                                                                       cTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTaaa                                                                            <  2:1215638/148‑1 (MQ=255)
                                                                        tATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATGCCGCTGGGGGATTGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAg                                                                           <  1:1570864/148‑1 (MQ=255)
                                                                         aTGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGc                                                                          <  2:3294213/148‑1 (MQ=255)
                                                                          tGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCg                                                                         >  1:2216451/1‑148 (MQ=255)
                                                                           gAAGTGAAAAAGCTGGAACGCACGCGCATAATGAAGGTTAGCTTAAGGGGCTTTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGa                                                                        <  1:27698/148‑1 (MQ=255)
                                                                              ttGAAAAAGCTGGAACGCACGCGCATCATTAACGTTAGCTTAAGCGTCATTCCGCTGGGGGAATGGCTCTATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAgg                                                                     <  2:3382041/147‑1 (MQ=255)
                                                                              gTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAgg                                                                     <  2:3382774/148‑1 (MQ=255)
                                                                              gTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAgg                                                                     <  2:471288/148‑1 (MQ=255)
                                                                               ttaaaaaGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGc                                                                    <  2:1917017/146‑1 (MQ=255)
                                                                                gAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGcc                                                                   >  1:963616/1‑148 (MQ=255)
                                                                                  aaaaGCTGCAACGCACGCGAATCATGAACGTTAGCTTACGCGGCTTTCCGGTGGGGGATTGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCaa                                                                 <  1:3256484/148‑1 (MQ=255)
                                                                                   aaaGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCaaa                                                                >  1:643573/1‑148 (MQ=255)
                                                                                   aaaGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCaaa                                                                >  2:2586049/1‑148 (MQ=255)
                                                                                   aaaGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCaaa                                                                >  1:2632822/1‑148 (MQ=255)
                                                                                    aaGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAAc                                                               >  2:933135/1‑148 (MQ=255)
                                                                                     aGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAAcc                                                              >  1:1978316/1‑148 (MQ=255)
                                                                                     aGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAAcc                                                              >  1:71698/1‑148 (MQ=255)
                                                                                       ctgcaccgcacgcGAATCTTGAAGGTTAGTTTAAGGGGCTTTCCGTTGGGGGATTGGCGCTATTAAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGa                                                            <  1:1236304/142‑1 (MQ=255)
                                                                                       cTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGa                                                            <  2:3167895/148‑1 (MQ=255)
                                                                                          gAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGaaaa                                                         <  2:567123/148‑1 (MQ=255)
                                                                                          gAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGaaaa                                                         <  2:1365372/148‑1 (MQ=255)
                                                                                            acgcacgcGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCTATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGACAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACa                                                       >  2:1305400/1‑148 (MQ=255)
                                                                                             cgcacgcGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGATTGGGGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAg                                                      <  1:3353360/148‑1 (MQ=255)
                                                                                               cacgcGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCg                                                    <  1:815594/148‑1 (MQ=255)
                                                                                               cacgcGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCg                                                    >  1:657661/1‑148 (MQ=255)
                                                                                               cacgcGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCg                                                    >  1:2420507/1‑148 (MQ=255)
                                                                                                 cgcgCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCggg                                                  >  2:2914248/1‑148 (MQ=255)
                                                                                                   cgcATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCa                                                >  2:2146809/1‑148 (MQ=255)
                                                                                                   cgcATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCa                                                <  2:2458176/148‑1 (MQ=255)
                                                                                                           gAACGTTAGTTTAAGGGGATTTCGGGTGGGGTATTGGCGCTATTTAACCGAGGATGATTTAATCGACATATTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGt                                        <  1:2797620/148‑1 (MQ=255)
                                                                                                           gAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGt                                        >  2:267256/1‑148 (MQ=255)
                                                                                                              cGTTAGCTTAAGCGGCATTCCGCTGGGGTAATGGCGCTATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCa                                     <  2:2391002/148‑1 (MQ=255)
                                                                                                                ttAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGt                                   >  1:1003410/1‑148 (MQ=255)
                                                                                                                ttAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGt                                   >  1:1015039/1‑148 (MQ=255)
                                                                                                                ttAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGt                                   >  1:3246211/1‑148 (MQ=255)
                                                                                                                 tAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTc                                  >  1:2656687/1‑148 (MQ=255)
                                                                                                                   gCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCACACGTCCAGTCGt                                >  1:1883044/1‑148 (MQ=255)
                                                                                                                     ttAAGCGGCATTCCGCTGGGGTAATGGCGCTATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTCGTTa                              <  2:1844320/148‑1 (MQ=255)
                                                                                                                     ttAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTCGTTa                              <  2:1444376/148‑1 (MQ=255)
                                                                                                                         gCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTCGTTAAGAt                          >  2:3175342/1‑148 (MQ=255)
                                                                                                                         gCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTCGTTAAGAt                          <  2:401471/148‑1 (MQ=255)
                                                                                                                             cATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTCGTTAAGATGGaa                      >  1:211579/1‑148 (MQ=255)
                                                                                                                              tttCCGCTGGGGGAATGGCGAGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTAATTTAAAATTACTCTTCCTAGGTAAAACATAAAGAGAAGTCCAAACCGAAAACAGCTGGCATAAAACGTCCAGTCGTTAAGATGGaaa                     <  2:2211994/147‑1 (MQ=255)
                                                                                                                                tCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTCGTTAAGATGGaaaaa                   >  1:2172833/1‑148 (MQ=255)
                                                                                                                                tCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTCGTTAAGATGGaaaaa                   >  1:2629778/1‑148 (MQ=255)
                                                                                                                                    cTGGGGGAAGGGCGCTATTAACCCGACCATGAGTAAATCGACCTCTTTAAGCTCATTTAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGCCCAAACCGAAAACCGCGGGCATCAAACGTCCAGTCGTTAAGATGGAAAAAAcgg               <  1:2868028/148‑1 (MQ=255)
                                                                                                                                       ggggAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTCGTTAAGATGGAAAAAAcggcgg            >  1:1905263/1‑148 (MQ=255)
                                                                                                                                           atgggCGCTATTAAACCGACTATGATTTAATCTACCTCTTTACGCTATTTTAAAATTCCTCTTCCGAGGTAAAACATAACGCGACGGCCAAACCGAAAACCGCGGGCATCAAACGTCCAGTCGTTAAGATGGAAAAAACGGCGGaaaa        <  1:1627927/145‑1 (MQ=255)
                                                                                                                                             tGGCGCGATTTACCCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTCGTTAAGATGGAAAAAACGGCGGAAAAAg      <  2:396534/148‑1 (MQ=255)
                                                                                                                                             tGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTCGTTAAGATGGAAAAAACGGCGGAAAAAg      >  1:2764588/1‑148 (MQ=255)
                                                                                                                                              ggCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTCGTTAAGATGGAAAAAACGGCGGAAAAAgg     <  2:3221779/148‑1 (MQ=255)
                                                                                                                                                cgctATTAAACCGACGATGAGTTAATCGACATCTTTAAGCACATTGAAAATTCCTCTTCCGAGTTAAAACCTAAAGCGAAGGCAAAACCGAAAACAGCGGGCATCAAACGTCCAGTCGTTAAGATGGAAAAAACGGCGGAAAAGGCgg  <  1:76998/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CCGTTTGTCTTCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTCGTTAAGATGGAAAAAACGGCGGAAAAAGGCGG  >  CP009273/4220801‑4221093

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