Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,957,997 C→A 15.8% D332Y (GAT→TAT)  rhlB ← ATP‑dependent RNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,957,9970CA15.8% 580.7 / 59.9 190D332Y (GAT→TAT) rhlBATP‑dependent RNA helicase
Reads supporting (aligned to +/- strand):  ref base C (88/72);  new base A (21/9);  total (109/81)
Fisher's exact test for biased strand distribution p-value = 1.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 190 total aligned reads displayed.

CGGAATTGAGTGACCAATATAGGTCTCAATAGCAGGCAAATTCAATGCATACTCTTCACACGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGTCGGTGGCAACCAGAATATCCAGATCGCCACGGGTAAATTCATCAAGAATACGCAGACGTTTTTTCTGCGCGACATCGCCTGTCAATAAACCG  >  CP009273/3957851‑3958142
                                                                                                                                                  |                                                                                                                                                 
cGGAATTGAGTGACCAATATAGGTCTCAATAGCAGGCAAATTCAATGCATACTCTTCACACGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATCg                                                                                                                                                  <  2:1460639/148‑1 (MQ=255)
 ggAATTGAGTGACCAATATAGGTCTCAATAGCAGGCAAATTCAATGCATACTCTTCACACGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATCGt                                                                                                                                                 <  1:483509/148‑1 (MQ=255)
 ggAATTGAGTGACCAATATAGGTCTCAATAGCAGGCAAATTCAATGCATACTCTTCACACGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATCGt                                                                                                                                                 <  2:572804/148‑1 (MQ=255)
 ggAATTGAGTGACCAATATAGGTCTCAATAGCAGGCAAATTCAATGCATACTCTTCACACGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATAGt                                                                                                                                                 <  2:1782964/148‑1 (MQ=255)
  gAATTGAGTGACCAATATAGGTCTCAATAGCAGGCAAATTCAATGCATACTCTTCACACGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATCGTa                                                                                                                                                >  1:1837222/1‑148 (MQ=255)
    aTTGAGTGACCAATATAGGTCTCAATAGCAGGCAAATTCAATGCATACTCTTCACACGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATCGTAGt                                                                                                                                              <  2:159034/148‑1 (MQ=255)
      tGAGTGACCAATATAGGTCTCAATAGCAGGCAAATTCAATGCATACTCTTCACACGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATCGTAGTTa                                                                                                                                            >  1:2120857/1‑148 (MQ=255)
      tGAGTGACCAATATAGGTCTCAATAGCAGGCAAATTCAATGCATACTCTTCACACGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATCGTAGTTa                                                                                                                                            <  2:421730/148‑1 (MQ=255)
          tGACCAATATAGGTCTCAATAGCAGGCAAATTCAATGCATACTCTTCACACGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGa                                                                                                                                        <  1:1171431/148‑1 (MQ=255)
          tGACCAATATAGGTCTCAATAGCAGGCAAATTCAATGCATACTCTTCACACGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATCGTAGTAAAAGa                                                                                                                                        <  2:2615641/148‑1 (MQ=255)
            aCCAATATAGGTCTCAATAGCAGGCAAATTCAATGCATACTCTTCACACGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACg                                                                                                                                      >  1:2249706/1‑148 (MQ=255)
                   taGGTCTCAATAGCAGGCAAATTCAATGCATACTCTTCACACGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCa                                                                                                                               >  1:2923892/1‑148 (MQ=255)
                   taGGTCTCAATAGCAGGCAAATTCAATGCATACTCTTCACACGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCa                                                                                                                               <  2:1600914/148‑1 (MQ=255)
                   taGGTCTCAATAGCAGGCAAATTCAATGCATACTCTTCACACGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATAGTAGTTAAAGACGTGCGTCa                                                                                                                               <  2:81808/148‑1 (MQ=255)
                    aGGTCTCAATAGCAGGCAAATTCAATGCATACTCTTCACACGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCAc                                                                                                                              >  2:1577805/1‑148 (MQ=255)
                     ggTCTCAATAGCAGGCAAATTCAATGCATACTCTTCACACGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACt                                                                                                                             >  1:2545155/1‑148 (MQ=255)
                      gTCTCAATAGCAGGCAAATTCAATGCATACTCTTCACACGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTg                                                                                                                            >  2:2862050/1‑148 (MQ=255)
                          cAATAGCAGGCAAATTCAATGCATACTCTTCACACGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATCGTAGTAAAAGACGTGCGTCACTGCCgg                                                                                                                        >  1:1928713/1‑148 (MQ=255)
                             tAGCAGGCAAATTCAATGCATACTCTTCACACGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAat                                                                                                                     <  2:923506/148‑1 (MQ=255)
                                  ggCAAATTCAATGCATACTCTTCACACGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCa                                                                                                                >  1:2568466/1‑148 (MQ=255)
                                   gCAAATTCAATGCATACTCTTCACACGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATAGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCaa                                                                                                               >  2:3091177/1‑148 (MQ=255)
                                        ttCAATGCATACTCTTCACACGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCAc                                                                                                          >  2:2081430/1‑148 (MQ=255)
                                            aTGCATACTCTTCACACGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCAcgcgc                                                                                                      >  1:2098650/1‑148 (MQ=255)
                                            aTGCATACTCTTCACACGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCAcgcgc                                                                                                      >  1:562150/1‑148 (MQ=255)
                                            aTGCATACTCTTCACACGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCAcgcgc                                                                                                      >  2:290874/1‑148 (MQ=255)
                                             tGCATACTCTTCACACGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCAcgcgcg                                                                                                     >  1:1987893/1‑148 (MQ=255)
                                             tGCATACTCTTCACACGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCAcgcgcg                                                                                                     <  1:2110291/148‑1 (MQ=255)
                                              gCATACTCTTCACACGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCgg                                                                                                    >  1:124904/1‑148 (MQ=255)
                                                aTACTCTTCACACGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCa                                                                                                  <  2:1090555/148‑1 (MQ=255)
                                                aTACTCTTCACACGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCa                                                                                                  <  2:2085480/148‑1 (MQ=255)
                                                 tACTCTTCACACGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATAGTCGTTAAAGACGTGCGTCACTGCCGGACTATGAAAACCACGCGCTGCaa                                                                                                 >  2:2553719/1‑148 (MQ=255)
                                                 tACTCTTCACACGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATAGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCaa                                                                                                 >  2:899879/1‑148 (MQ=255)
                                                   ctctTCACACGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATAGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACg                                                                                               >  2:1678516/1‑148 (MQ=255)
                                                   ctctTCACACGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATAGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACg                                                                                               >  2:366768/1‑148 (MQ=255)
                                                    tctTCACACGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGt                                                                                              <  2:2191982/148‑1 (MQ=255)
                                                     ctTCACACGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGTc                                                                                             >  2:679601/1‑148 (MQ=255)
                                                     ctTCACACGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGTc                                                                                             >  2:679405/1‑148 (MQ=255)
                                                     ctTCACACGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATAGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGTc                                                                                             >  2:1528596/1‑148 (MQ=255)
                                                      ttCACACGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGTCg                                                                                            <  1:2302680/148‑1 (MQ=255)
                                                      ttCACACGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATAGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGTCg                                                                                            >  1:1451357/1‑148 (MQ=255)
                                                       tCACACGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGTCgg                                                                                           >  1:1558119/1‑148 (MQ=255)
                                                       tCACACGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGTCgg                                                                                           >  1:1192466/1‑148 (MQ=255)
                                                          cacGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGTCGGTgg                                                                                        >  2:1945109/1‑148 (MQ=255)
                                                          cacGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGTCGGTgg                                                                                        >  2:2984664/1‑148 (MQ=255)
                                                            cGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGACCGTAATCTTCACAGTCATCGGGTAAATCGTAGTAAAAGACGTGCGTCACTGCCGGACTATGCCAACCACGCGCGGCAACGTCGGTGGCa                                                                                      >  1:2416722/1‑148 (MQ=255)
                                                            cGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGTCGGTGGCa                                                                                      <  1:3086109/148‑1 (MQ=255)
                                                            cGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGTCGGTGGCa                                                                                      <  1:1757541/148‑1 (MQ=255)
                                                            cGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATCGTAGTAAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGTCGGTGGCa                                                                                      >  1:764068/1‑148 (MQ=255)
                                                               cAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGTCGGTGGCAAcc                                                                                   >  2:1490345/1‑148 (MQ=255)
                                                                 ggCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGTCGGTGGCAACCAg                                                                                 >  2:1720170/1‑148 (MQ=255)
                                                                  gCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGTCGGTGGCAACCAGa                                                                                >  2:2961033/1‑148 (MQ=255)
                                                                    tGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGTCGGTGGCAACCAGAat                                                                              <  1:259449/148‑1 (MQ=255)
                                                                      aTAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGTCGGTGGCAACCAGAatat                                                                            <  2:2410842/148‑1 (MQ=255)
                                                                      aTAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGTCGGTGGCAACCAGAatat                                                                            >  2:3214029/1‑148 (MQ=255)
                                                                        agagTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGTCGGTGGCAACCAGAATATcc                                                                          <  1:1491714/148‑1 (MQ=255)
                                                                          agTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGTCGGTGGCAACCAGAATATCCAg                                                                        <  2:2660858/148‑1 (MQ=255)
                                                                            tGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGTCGGTGGCAACCAGAATATCCAGAt                                                                      >  1:1445573/1‑148 (MQ=255)
                                                                              aCCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGTCGGTGGCAACCAGAATATCCAGATCg                                                                    >  2:86483/1‑148 (MQ=255)
                                                                               ccGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGTCGGTGGCAACCAGAATATCCAGATCGc                                                                   <  1:1825001/148‑1 (MQ=255)
                                                                                cGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGTCGGTGGCAACCAGAATATCCAGATCGcc                                                                  >  1:292498/1‑148 (MQ=255)
                                                                                  cTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGTCGGTGGCAACCAGAATATCCAGATCGCCAc                                                                >  1:3240245/1‑148 (MQ=255)
                                                                                  cTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGTCGGTGGCAACCAGAATATCCAGATCGCCAc                                                                >  1:3240273/1‑148 (MQ=255)
                                                                                      cgcCTGCGCGACCTGTACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGTCGGTGGCAACCAGAATATCCAGATCGCCACGGGt                                                            >  1:854611/1‑148 (MQ=255)
                                                                                         cTGCGCGACCTGTACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGTCGGTGGCAACCAGAATATCCAGATCGCCACGGGTaaa                                                         <  1:1516494/148‑1 (MQ=255)
                                                                                          tGCGCGACCTGTACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGTCGGTGGCAACCAGAATATCCAGATCGCCACGGGTAAAt                                                        <  2:1973939/148‑1 (MQ=255)
                                                                                           gcgcGACCTGTACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGTCGGTGGCAACCAGAATATCCAGATCGCCACGGGTAAAtt                                                       >  2:3253827/1‑148 (MQ=255)
                                                                                              cgACCTGTACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGTCGGTGGCAACCAGAATATCCAGATCGCCACGGGTAAATtcat                                                    >  1:1063244/1‑148 (MQ=255)
                                                                                              cgACCTGTACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATCGTAGTAAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGTCGGTGGCAACCAGAATATCCAGATCGCCACGGGTAAATtcat                                                    <  1:1361853/148‑1 (MQ=255)
                                                                                                aCCTGTACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGTCGGTGGCAACCAGAATATCCAGATCGCCACGGGTAAATtcatca                                                  >  1:1669299/1‑148 (MQ=255)
                                                                                                aCCTGTACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGTCGGTGGCAACCAGAATATCCAGATCGCCACGGGTAAATtcatca                                                  >  2:642942/1‑148 (MQ=255)
                                                                                                aCCTGTACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGTCGGTGGCAACCAGAATATCCAGATCGCCACGGGTAAATtaatca                                                  <  1:72048/148‑1 (MQ=255)
                                                                                                     tACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGTCGGTGGCAACCAGAATATCCAGATCGCCACGGGTAAATTCATCAAGAAt                                             <  1:3217684/148‑1 (MQ=255)
                                                                                                        gACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGTCGGTGGCAACCAGAATATCCAGATCGCCACGGGTAAATTCATCAAGAATACg                                          >  2:2995799/1‑148 (MQ=255)
                                                                                                          ccAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGTCGGTGGCAACCAGAATATCCAGATCGCCACGGGTAAATTCATCCAGAATACGCa                                        >  2:2377699/1‑148 (MQ=255)
                                                                                                           cAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGTCGGTGGCAACCAGAATATCCAGATCGCCACGGGTAAATTCATCAAGAATACGCAg                                       >  1:1021762/1‑148 (MQ=255)
                                                                                                           cAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATAGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGTCGGTGGCAACCAGAATATCCAGATCGCCACGGGTAAATTCATCAAGAATACGCAg                                       >  2:679799/1‑148 (MQ=255)
                                                                                                             aTACGGTGAACGTAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGTCGGTGGCAACCAGAATATCCAGATCGCCACGGGTAAATTCATCAAGAATACGCAGAc                                     >  1:2687316/1‑148 (MQ=255)
                                                                                                              tACGGTGAACGTAATCTTCACAGTCATCTGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGTCGGTGGCAACCAGAATATCCAGATCGCCACGGGTAAATTCATCAAGAATACGCAGACg                                    <  1:2827174/148‑1 (MQ=255)
                                                                                                              tACGGTGAACGTAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGTCGGTGGCAACCAGAATATCCAGATCGCCACGGGTAAATTCATCAAGAATACGCAGACg                                    >  1:1030001/1‑148 (MQ=255)
                                                                                                                 ggTGAACGTAATCTTCGCAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGTCGGTGGCAACCAGAATATCCAGATCGCCACGGGTAAATTCATCAAGAATACGCAGACGttt                                 <  1:1802518/148‑1 (MQ=255)
                                                                                                                 ggTGAACGTAATCTTCACAGTCATCGGGTAAATAGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGTCGGTGGCAACCAGAATATCCAGATCGCCACGGGTAAATTCATCAAGAATACGCAGACGttt                                 <  1:1611181/148‑1 (MQ=255)
                                                                                                                   tGAACGTAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGTCGGTGGCAACCAGAATATCCAGATCGCCACGGGTAAATTCATCAAGAATACGCAGACGttttt                               <  1:2786292/148‑1 (MQ=255)
                                                                                                                   tGAACGTAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGTCGGTGGCAACCAGAATATCCAGATCGCCACGGGTAAATTCATCAAGAATACGCAGACGttttt                               <  1:662706/148‑1 (MQ=255)
                                                                                                                    gAACGTAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGTCGGTGGCAACCAGAATATCCAGATCGCCACGGGTAAATTCATCAAGAATACGCAGACGtttttt                              <  2:2402213/148‑1 (MQ=255)
                                                                                                                     aaCGTAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGTCGGTGGCAACCAGAATATCCAGATCGCCACGGGTAAATTCATCAAGAATACGCAGACGTTTTTTc                             <  2:1837222/148‑1 (MQ=255)
                                                                                                                     aaCGTAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGTCGGTGGCAACCAGAATATCCAGATCGCCACGGGTAAATTCATCAAGAATACGCAGACGTTTTTTc                             >  2:2472136/1‑148 (MQ=255)
                                                                                                                          aaTCTTCACAGTCATCGGGTAAATAGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGTCGGTGGCAACCAGAATATCCAGATCGCCACGGGTCAATTCATCAAGAATACGCAGACGTTTTTTCTgcgc                        >  2:1361302/1‑148 (MQ=255)
                                                                                                                            tCTTCACAGTCATCGGGTAAATAGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGTCGGTGGCAACCAGAATATCCAGATCGCCACGGGTAAATTCATCAAGAATACGCAGACGTTTTTTCTGCGCGa                      >  1:2492776/1‑148 (MQ=255)
                                                                                                                                cacaGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGTCGGTGGCAACCAGAATATCCAGATCGCCACGGGTAAATTCATCAAGAATACGCAGACGTTTTTTCTGCGCGACATc                  >  2:115187/1‑148 (MQ=255)
                                                                                                                                   aGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGTCGGTGGCAACCAGAATATCCAGATCGCCACGGGTAAATTCATCAAGAATACGCAGACGTTTTTTCTGCGCGACATCGcc               <  2:1653349/148‑1 (MQ=255)
                                                                                                                                      cATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGTCGGTGGCAACCAGAATATCCAGATCGCCACGGGTAAATTCATCAAGAATACGCAGACGTTTTTTCTGCGCGACATCGCCTGt            >  2:108692/1‑148 (MQ=255)
                                                                                                                                      cATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGTCGGTGGCAACCAGAATATCCAGATCGCCACGGGTAAATTCATCAAGAATACGCAGACGTTTTTTCTGCGCGACATCGCCTGt            >  2:62735/1‑148 (MQ=255)
                                                                                                                                      cATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGTCGGTGGCAACCAGAATATCCAGATCGCCACGGGTAAATTCATCAAGAATACGCAGACGTTTTTTCTGCGCGACATCGCCTGt            <  2:771764/148‑1 (MQ=255)
                                                                                                                                        tCGGGTAAATCGTAGTTTAAGACGTGCGTAACTGCCGGAATATGCAAACCACGCGCGGCAACGTCGGTGGCAACAAGAATATCCAGATCGCCACGGGTAAATTCATCAAGAATACGCAGACGTTTTTTCTGCGCGACATCGCCTGTCa          <  1:3150031/148‑1 (MQ=255)
                                                                                                                                             tAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGTCGGTGGCAACCAGAATATCCAGATCGCCACGGGTAAATTCATCAAGAATACGCAGACGTTTTTTCTGCGCGACATCGCCTGTCAATaaa     <  1:1397401/148‑1 (MQ=255)
                                                                                                                                             tAAATAGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGTCGGTGGCAACCAGAATATCCAGATCGCCACGGGTAAATTCATCAAGAATACGCAGACGTTTTTTCTGCGCGACATCGCCTGTCAATaaa     <  2:2236018/148‑1 (MQ=255)
                                                                                                                                              aaaTCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGTCGGTGGCAACCAGAATATCCAGATCGCCACGGGTAAATTCATCAAGAATACGCAGACGTTTTTTCTGCGCGACATCGCCTGTCAATAAAc    >  2:750588/1‑148 (MQ=255)
                                                                                                                                              aaaTCGTAGTAAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGTCGGTGGCAACCAGAATATCCAGATCGCCACGGGTAAATTCATCAAGAATACGCAGACGTTTTTTCTGCGCGACATCGCCTGTCAATAAAc    >  2:844647/1‑148 (MQ=255)
                                                                                                                                              aaaTAGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGTCGGTGGCAACCAGAATATCCAGATCGCCACGGGTAAATTCATCAAGAATACGCAGACGTTTTTTCTGCGCGACATCGCCTGTCAATAAAc    >  2:2352486/1‑148 (MQ=255)
                                                                                                                                               aaTCGTAGTTAAAGACGAGCGTCACTGCCGGAATATGCAAAACACGCGCGGCAACGTCGGTGGAAACCAGAATATCCAGATCGCCACGGGTAAATACATCAAGAATACGCAGACGTTTTTTCTCCTCGACATCGCCTGACAATaaaac   >  1:193655/1‑146 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
CGGAATTGAGTGACCAATATAGGTCTCAATAGCAGGCAAATTCAATGCATACTCTTCACACGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGTCGGTGGCAACCAGAATATCCAGATCGCCACGGGTAAATTCATCAAGAATACGCAGACGTTTTTTCTGCGCGACATCGCCTGTCAATAAACCG  >  CP009273/3957851‑3958142

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: