Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,122,778 A→C 100% pseudogene (725/1101 nt) yghO ← pseudogene, DNA‑binding transcriptional regulator homology

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,122,7780AC96.6% 803.3 / 4.5 203pseudogene (725/1101 nt)yghOpseudogene, DNA‑binding transcriptional regulator homology
Reads supporting (aligned to +/- strand):  ref base A (4/3);  new base C (91/105);  total (95/108)
Fisher's exact test for biased strand distribution p-value = 7.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 204 total aligned reads displayed.

GTTCGCACACCGCTGACTTTCAAACGCCACAGCAATTTTGCCCAGCCGAAGGGAAAGAGCGATCCGTTCAGATCGGCAATCGCCTCGTTGATGTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGATCTGCATCTCTTCGGCAAACCGCT  >  CP009273/3122631‑3122919
                                                                                                                                                   |                                                                                                                                             
   cGCACACCGCTGACTTTCAAACGCCACAGCAATTTTGCCCAGCCGAAGGGAAAGAGCGATCCGTTCAGATCGGCAATCGCCTCGTTGATGTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCa                                                                                                                                            <  2:811923/148‑1 (MQ=255)
         ccGCTGACTTTCAAACGCCACAGCAATTTTGCCCAGCCGAAGGGAAAGAGCGATCCGTTCAGATCGGCAATCGCCTCGTTGATGTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGt                                                                                                                                      <  1:3505199/148‑1 (MQ=255)
         ccGCTGACTTTCAAACGCCACAGCAATTTTGCCCAGCCGAAGGGAAAGAGCGATCCGTTCAGATCGGCAATCGCCTCGTTGATGTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGt                                                                                                                                      <  1:1909586/148‑1 (MQ=255)
         ccGCTGACTTTCAAACGCCACAGCAATTTTGCCCAGCCGAAGGGAAAGAGCGATCCGTTCAGATCGGCAATCGCCTCGTTGATGTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGt                                                                                                                                      <  2:2796783/148‑1 (MQ=255)
          cGCTGACTTTCAAACGCCACAGCAATTTTGCCCAGCCGAAGGGAAAGAGCGATCCGTTCAGATCGGCAATCGCCTCGTTGATGTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTc                                                                                                                                     <  1:1294836/148‑1 (MQ=255)
          cGCTGACTTTCAAACGCCACAGCAATTTTGCCCAGCCGAAGGGAAAGAGCGATCCGTTCAGATCGGCAATCGCCTCGTTGATGTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTc                                                                                                                                     <  1:514571/148‑1 (MQ=255)
           gCTGACTTTCAAACGCCACAGCAATTTTGCCCAGCCGAAGGGAAAGAGCGATCCGTTCAGATCGGCAATCGCCTCGTTGATGTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTcc                                                                                                                                    >  2:2655253/1‑148 (MQ=255)
            cTGACTTTCAAACGCCACAGCAATTTTGCCCAGCCGAAGGGAAAGAGCGATCCGTTCAGATCGGCAATCGCCTCGTTGATGTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCg                                                                                                                                   >  1:3414452/1‑148 (MQ=255)
              gACTTTCAAACGCCACAGCAATTTTGCCCAGCCGAAGGGAAAGAGCGATCCGTTCAGATCGGCAATCGCCTCGTTGATGTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGc                                                                                                                                 >  2:2620224/1‑148 (MQ=255)
                 tttCAAACGCCACAGCAATTTTGCCCAGCCGAAGGGAAAGAGCGATCCGTTCAGATCGGCAATCGCCTCGTTGATGTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCAcc                                                                                                                              <  2:2858287/148‑1 (MQ=255)
                 tttCAAACGCCACAGCAATTTTGCCCAGCCGAAGGGAAAGAGCGATCCGTTCAGATCGGCAATCGCCTCGTTGATGTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCAcc                                                                                                                              >  2:2504670/1‑148 (MQ=255)
                  ttCAAACGCCACAGCAATTTTGCCCAGCCGAAGGGAAAGAGCGATCCGTTCAGATCGGCAATCGCCTCGTTGATGTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCa                                                                                                                             >  1:794907/1‑148 (MQ=255)
                   tCAAACGCCACAGCAATTTTGCCCAGCCGAAGGGAAAGAGCGATCCGTTCAGATCGGCAATCGCCTCGTTGATGTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCAg                                                                                                                            <  2:348901/148‑1 (MQ=255)
                   tCAAACGCCACAGCAATTTTGCCCAGCCGAAGGGAAAGAGCGATCCGTTCAGATCGGCAATCGCCTCGTTGATGTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCAg                                                                                                                            <  1:471413/148‑1 (MQ=255)
                    cAAACGCCACAGCAATTTTGCCCAGCCGAAGGGAAAGAGCGATCCGTTCAGATCGGCAATCGCCTCGTTGATGTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCAGa                                                                                                                           <  2:851788/148‑1 (MQ=255)
                     aaaCGCCACAGCAATTTTGCCCAGCCGAAGGGAAAGAGCGATCCGTTCAGATCGGCAATCGCCTCGTTGATGTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCAGat                                                                                                                          <  2:3414452/148‑1 (MQ=255)
                      aaCGCCACAGCAATTTTGCCCAGCCGAAGGGAAAGAGCGATCCGTTCAGATCGGCAATCGCCTCGTTGATGTTCGGCAAGCCTACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCAGata                                                                                                                         <  2:1422258/148‑1 (MQ=255)
                      aaCGCCACAGCAATTTTGCCCAGCCGAAGGGAAAGAGCGATCCGTTCAGATCGGCAATCGCCTCGTTGATGTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCAGata                                                                                                                         >  2:2994395/1‑148 (MQ=255)
                       aCGCCACAGCAATTTTGCCCAGCCGAAGGGAAAGAGCGATCCGTTCAGATCGGCAATCGCCTCGTTGATGTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCAGatat                                                                                                                        >  2:2805661/1‑148 (MQ=255)
                       aCGCCACAGCAATTTTGCCCAGCCGAAGGGAAAGAGCGATCCGTTCAGATCGGCAATCGCCTCGTTGATGTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCAGatat                                                                                                                        >  2:2331166/1‑148 (MQ=255)
                       aCGCCACAGCAATTTTGCCCAGCCGAAGGGAAAGAGCGATCCGTTCAGATCGGCAATCGCCTCGTTGATGTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCAGatat                                                                                                                        >  2:152887/1‑148 (MQ=255)
                         gCCACAGCAATTTTGCCCAGCCGAAGGGAAAGAGCGATCCGTTCAGATCGGCAATCGCCTCGTTGATGTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCAGATAttt                                                                                                                      <  1:3032300/148‑1 (MQ=255)
                           cacaGCAATTTTGCCCAGCCGAAGGGAAAGAGCGATCCGTTCAGATCGGCAATCGCCTCGTTGATGTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCAGATATTTaa                                                                                                                    <  2:1626540/148‑1 (MQ=255)
                           cacaGCAATTTTGCCCAGCCGAAGGGAAAGAGCGATCCGTTCAGATCGGCAATCGCCTCGTTGATGTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCAGATATTTaa                                                                                                                    <  1:3292275/148‑1 (MQ=255)
                             caGCAATTTTGCCCAGCCGAAGGGAAAGAGCGATCCGTTCAGATCGGCAATCGCCTCGTTGATGTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCAGATATTTAAGt                                                                                                                  >  2:605424/1‑148 (MQ=255)
                             caGCAATTTTGCCCAGCCGAAGGGAAAGAGCGATCCGTTCAGATCGGCAATCGCCTCGTTGATGTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGt                                                                                                                  >  1:2784919/1‑148 (MQ=255)
                              aGCAATTTTGCCCAGCCGAAGGGAAAGAGCGATCCGTTCAGATCGGCAATCGCCTCGTTGATGTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCAGATATTTAAGtt                                                                                                                 >  2:1967825/1‑148 (MQ=255)
                               gCAATTTTGCCCAGCCGAAGGGAAAGAGCGATCCGTTCAGATCGGCAATCGCCTCGTTGATGTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCAGATATTTAAGTTg                                                                                                                <  2:2456185/148‑1 (MQ=255)
                                cAATTTTGCCCAGCCGAAGGGAAAGAGCGATCCGTTCAGATCGGCAATCGCCTCGTTGATGTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCAGATATTTAAGTTGa                                                                                                               <  2:798084/148‑1 (MQ=255)
                                 aaTTTTGCCCAGCCGAAGGGAAAGAGCGATCCGTTCAGATCGGCAATCGCCTCGTTGATGTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCAGATATTTAAGTTGAt                                                                                                              <  1:3282487/148‑1 (MQ=255)
                                 aaTTTTGCCCAGCCGAAGGGAAAGAGCGATCCGTTCAGATCGGCAATCGCCTCGTTGATGTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCAGATATTTAAGTTGAt                                                                                                              <  1:3325591/148‑1 (MQ=255)
                                           aGCCGAAGGGAAAGAGCGATCCGTTCAGATCGGCAATCGCCTCGTTGATGTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGt                                                                                                    >  1:2876653/1‑148 (MQ=255)
                                            gCCGAAGGGAAAGAGCGATCCGTTCAGATCGGCAATCGCCTCGTTGATGTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTc                                                                                                   >  1:264401/1‑148 (MQ=255)
                                             ccGAAGGGAAAGAGCGATCCGTTCAGATCGGCAATCGCCTCGTTGATGTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTcg                                                                                                  >  1:2069247/1‑148 (MQ=255)
                                                aaGGGAAAGAGCGATCCGTTCAGATCGGCAATCGCCTCGTTGATGTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGa                                                                                               <  2:1809879/148‑1 (MQ=255)
                                                  gggAAAGAGCGATCCGTTCAGATCGGCAATCGCCTCGTTGATGTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGAAt                                                                                             >  2:2916818/1‑148 (MQ=255)
                                                  gggAAAGAGCGATCCGTTCAGATCGGCAATCGCCTCGTTGATGTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGAAt                                                                                             >  2:2682336/1‑148 (MQ=255)
                                                  gggAAAGAGCGATCCGTTCAGATCGGCAATCGCCTCGTTGATGTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGAAt                                                                                             <  1:2619597/148‑1 (MQ=255)
                                                   ggAAAGAGCGATCCGTTCAGATCGGCAATCGCCTCGTTGATGTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGAAtt                                                                                            <  1:1157302/148‑1 (MQ=255)
                                                    gAAAGAGCGATCCGTTCAGATCGGCAATCGCCTCGTTGATGTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGAATTc                                                                                           <  1:2421593/148‑1 (MQ=255)
                                                    gAAAGAGCGATCCGTTCAGATCGGCAATCGCCTCGTTGATGTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGAATTc                                                                                           <  2:901888/148‑1 (MQ=255)
                                                    gAAAGAGCGATCCGTTCAGATCGGCAATCGCCTCGTTGATGTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCAGATATGTAAGTTGATCGCCCATGGTCGCGACTTc                                                                                           >  2:2228559/1‑148 (MQ=255)
                                                     aaaGAGCGATCCGTTCAGATCGGCAATCGCCTCGTTGATGTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGAATTCa                                                                                          <  2:2694454/148‑1 (MQ=255)
                                                       agagTGATCCGTTCAGATCGGCAATCGCCTCGTTGATGTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGAATTCATg                                                                                        <  1:2699706/148‑1 (MQ=255)
                                                       agagCGATCCGTTCAGATCGGCAATCGCCTCGTTGATGTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGAATTCATg                                                                                        >  1:2983504/1‑148 (MQ=255)
                                                       agagCGATCCGTTCAGATCGGCAATCGCCTCGTTGATGTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGAATTCATg                                                                                        <  1:1300745/148‑1 (MQ=255)
                                                        gagCGATCCGTTCAGATCGGCAATCGCCTCGTTGATGTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCCTATCGTCCGGCACCCGATATTTAAGTTGATCGCCCATGGTCGCGAATTCATGt                                                                                       >  2:363851/1‑148 (MQ=255)
                                                         agCGATCCGTTCAGATCGGCAATCGCCTCGTTGATGTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGAATTCATGtt                                                                                      <  1:918440/148‑1 (MQ=255)
                                                           cGATCCGTTCAGATCGGCAATCGCCTCGTTGATGTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGAATTCATGTTCg                                                                                    <  1:1675549/148‑1 (MQ=255)
                                                              tCCGTTCAGATCGGCAATCGCCTCGTTGATGTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGACTTCATGTACGGTa                                                                                 >  2:2124620/1‑148 (MQ=255)
                                                              tCCGTTCAGATCGGCAATCGCCTCGTTGATGTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGAATTCATGTTCGGTa                                                                                 >  2:3016151/1‑148 (MQ=255)
                                                              tCCGTTCAGATCGGCAATCGCCTCGTTGATGTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGAATTCATGTTCGGTa                                                                                 >  1:2105864/1‑148 (MQ=255)
                                                              tCCGTTCAGATCGGCAATCGCCTCGTTGATGTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGAATTCATGTTCGGTa                                                                                 >  2:1150883/1‑148 (MQ=255)
                                                              tCCGTTCAGATCGGCAATCGCCTCGTTGATGTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGAATTCATGTTCGGTa                                                                                 >  1:335675/1‑148 (MQ=255)
                                                               ccGTTCAGATCGGCAATCGCCTCGTTGATGTTCGACAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGAATTCATGTTCGGTaa                                                                                >  2:558185/1‑148 (MQ=255)
                                                                  ttCAGATCGGCAATCGCCTCGTTGATGTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGAATTCATGTTCGGTAAACg                                                                             <  2:1586519/148‑1 (MQ=255)
                                                                  ttCAGATCGGCAATCGCCTCGTTGATGTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGAATTCATGTTCGGTAAACg                                                                             <  2:507717/148‑1 (MQ=255)
                                                                  ttCAGATCGGCAATCGCCTCGTTGATGTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGAATTCATGTTCGGTAAACg                                                                             <  2:2242318/148‑1 (MQ=255)
                                                                   tCAGATCGGCAATCGCCTCGTTGATGTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCCGCGATATAGCTCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGAATTCATGTTCGGTAAACgg                                                                            <  1:676986/148‑1 (MQ=255)
                                                                   tCAGATCGGCAATCGCCTCGTTGATGTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGAATTCATGTTCGGTAAACgg                                                                            <  1:676962/148‑1 (MQ=255)
                                                                        tCGGCAATCGCCTCGTTGATGTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGAATTCATGTTCGGTAAACGGCACaa                                                                       >  1:2545262/1‑148 (MQ=255)
                                                                         cGGCAATCGCCTCGTTGATGTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGAATTCATGTTCGGTAAACGGCACaaa                                                                      <  2:2547402/148‑1 (MQ=255)
                                                                          ggCAATCGCCTCGTTGATGTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAAt                                                                     <  1:1867061/148‑1 (MQ=255)
                                                                            cAATCGCCTCGTTGATGTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATcc                                                                   <  1:545919/148‑1 (MQ=255)
                                                                               tCGCCTCGTTGATGTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCa                                                                <  1:1724079/148‑1 (MQ=255)
                                                                               tCGCCTCGTTGATGTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCa                                                                >  1:2918354/1‑148 (MQ=255)
                                                                                   cTCGTTGATGTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTtg                                                            <  2:1244677/148‑1 (MQ=255)
                                                                                    tCGTTGATGTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTtgt                                                           <  2:3150078/148‑1 (MQ=255)
                                                                                     cGTTGATGTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTtgtg                                                          >  1:2765882/1‑148 (MQ=255)
                                                                                        tGATGTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGAATTCATGTTCGGTACACGGCACAAATCCCCAGTTGtgctg                                                       <  2:1561766/148‑1 (MQ=255)
                                                                                          aTGTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTgcc                                                     >  1:594838/1‑148 (MQ=255)
                                                                                            gTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTgccag                                                   >  2:1995071/1‑148 (MQ=255)
                                                                                                  cAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCagag                                             <  1:2505440/148‑1 (MQ=255)
                                                                                                  cAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCagag                                             <  1:202050/148‑1 (MQ=255)
                                                                                                     gCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTg                                          >  2:3305952/1‑148 (MQ=255)
                                                                                                     gCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTg                                          <  2:1072321/148‑1 (MQ=255)
                                                                                                     gCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTg                                          >  1:1799819/1‑148 (MQ=255)
                                                                                                     gCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTg                                          >  1:1800380/1‑148 (MQ=255)
                                                                                                     gCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTg                                          >  1:2006834/1‑148 (MQ=255)
                                                                                                     gCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTg                                          >  1:234497/1‑148 (MQ=255)
                                                                                                     gCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTg                                          >  2:1720050/1‑148 (MQ=255)
                                                                                                       cGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGaa                                        <  1:1618727/148‑1 (MQ=255)
                                                                                                         aCAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGaaaa                                      >  2:419999/1‑148 (MQ=255)
                                                                                                          cAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGGCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAAt                                     <  1:2688630/148‑1 (MQ=255)
                                                                                                             tAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAAtctc                                  <  2:537025/148‑1 (MQ=255)
                                                                                                             tAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAAtctc                                  >  1:1525629/1‑148 (MQ=255)
                                                                                                             tAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAAtctc                                  >  1:317553/1‑148 (MQ=255)
                                                                                                               aaCGCGCAGGGTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCAc                                >  1:761706/1‑148 (MQ=255)
                                                                                                                       gggTGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGACTTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGagc                        >  2:1951067/1‑146 (MQ=255)
                                                                                                                          tGCAGAATCAATCTCAGCGATATAGCTCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCACCAGTTGTGCTGCCAGACAGAGTTGAAAATCTCACGCAGGCTCTGc                     >  2:1888290/1‑148 (MQ=255)
                                                                                                                               aatcaatcTCAGCGATATAGCTCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGATCTGCAtctc                >  2:3181471/1‑148 (MQ=255)
                                                                                                                               aatcaatcTCAGCGATATAGCTCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGATCTGCAtctc                >  2:1394011/1‑148 (MQ=255)
                                                                                                                                 tcaatcTCAGCGATATAGCTCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGATCTGCATCTCtt              <  2:2619376/148‑1 (MQ=255)
                                                                                                                                 tcaatcTCAGCGATATAGCTCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGAATTCATGTTCGGTAAAAGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGATCTGCATCTCtt              >  2:2165567/1‑148 (MQ=255)
                                                                                                                                 tcaatcTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGATCTGCATCTCtt              <  2:2213532/148‑1 (MQ=255)
                                                                                                                                     tctcAGCGATATAGCTCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATCGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGATCTGCATCTCTTCGGc          >  2:2176399/1‑148 (MQ=255)
                                                                                                                                      ctcAGCGATATAGCTCATATCGTCCGGCCCCAGATATTTAAGTTGATCGCCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGATCTGCATCTCTTCGGCa         >  2:3293759/1‑148 (MQ=255)
                                                                                                                                        cAGCGATATAGCTCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGATCTGCATCTCTTCGGCaaa       >  2:784846/1‑148 (MQ=255)
                                                                                                                                         aGCGATATAGCTCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGATCTGCATCTCTTCGGCAAAc      <  2:794907/148‑1 (MQ=255)
                                                                                                                                             atatAGCTCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGATCTGCATCTCTTCGGCAAACCGCt  >  2:691900/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                             
GTTCGCACACCGCTGACTTTCAAACGCCACAGCAATTTTGCCCAGCCGAAGGGAAAGAGCGATCCGTTCAGATCGGCAATCGCCTCGTTGATGTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGATCTGCATCTCTTCGGCAAACCGCT  >  CP009273/3122631‑3122919

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: