Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 958,594 G→T 17.6% G382C (GGC→TGC)  rpsA → 30S ribosomal subunit protein S1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273958,5940GT17.6% 469.0 / 60.0 165G382C (GGC→TGC) rpsA30S ribosomal subunit protein S1
Reads supporting (aligned to +/- strand):  ref base G (71/65);  new base T (13/16);  total (84/81)
Fisher's exact test for biased strand distribution p-value = 5.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

Only 100 of 166 total aligned reads displayed.

CTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGGCGACGAAATCGCTGCAGTTGTTCTGCAGGTTGACGCAGAACGTGAACGTATCT  >  CP009273/958447‑958738
                                                                                                                                                   |                                                                                                                                                
cTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCg                                                                                                                                                  <  2:2634880/148‑1 (MQ=255)
   gATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCt                                                                                                                                               <  1:755630/148‑1 (MQ=255)
   gATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCt                                                                                                                                               <  2:2773529/148‑1 (MQ=255)
       tCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGa                                                                                                                                           >  1:240963/1‑148 (MQ=255)
        cgacgaAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGAc                                                                                                                                          >  2:2649298/1‑148 (MQ=255)
        cgacgaAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGAc                                                                                                                                          >  2:1254310/1‑148 (MQ=255)
              agaaCGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACggcggc                                                                                                                                    >  2:570530/1‑148 (MQ=255)
                aaCGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCTGCTTGGACGGCGGCAt                                                                                                                                  >  2:661078/1‑148 (MQ=255)
                aaCGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCAt                                                                                                                                  >  2:2198992/1‑148 (MQ=255)
                       tcgtATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGc                                                                                                                           >  2:1378780/1‑148 (MQ=255)
                       tcgtATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGc                                                                                                                           >  1:1395598/1‑148 (MQ=255)
                        cgtATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGcc                                                                                                                          <  1:253123/148‑1 (MQ=255)
                          tATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACNGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCTGCTTGGACGGCGGCATCGACGGCCTg                                                                                                                        >  1:1679786/1‑148 (MQ=255)
                           aTCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTgg                                                                                                                       >  2:371933/1‑148 (MQ=255)
                            tctcCCTGGGTCTGAAACAGTGCAAAGCTAACGAGTGGCAGCAGTTCGCGGAAACCCAAAACAAGGGCGAGCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGt                                                                                                                      <  2:1685658/148‑1 (MQ=255)
                            tctcCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGt                                                                                                                      <  2:1685603/148‑1 (MQ=255)
                                ccTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCAc                                                                                                                  <  2:998339/148‑1 (MQ=255)
                                ccTGGGTCTGAAACAGTGCAAAGCTAACCCGAGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCAc                                                                                                                  >  1:1070092/1‑148 (MQ=255)
                                 cTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCAcc                                                                                                                 <  1:1714470/148‑1 (MQ=255)
                                  tGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCt                                                                                                                >  1:323338/1‑148 (MQ=255)
                                     gTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTc                                                                                                             <  2:1592579/148‑1 (MQ=255)
                                         gAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCTGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGAc                                                                                                         <  2:2540774/148‑1 (MQ=255)
                                         gAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGAc                                                                                                         >  1:1318415/1‑148 (MQ=255)
                                            aCAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACAtc                                                                                                      >  2:1683311/1‑148 (MQ=255)
                                              aGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACAtctc                                                                                                    <  2:2269239/148‑1 (MQ=255)
                                              aGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACAtctc                                                                                                    <  1:142079/148‑1 (MQ=255)
                                               gTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTcc                                                                                                   <  2:922610/148‑1 (MQ=255)
                                               gTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTcc                                                                                                   >  1:1984222/1‑148 (MQ=255)
                                               gTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTcc                                                                                                   >  1:1897607/1‑148 (MQ=255)
                                                    aaGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGaa                                                                                              >  1:3026379/1‑148 (MQ=255)
                                                     aGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAAc                                                                                             >  2:603386/1‑148 (MQ=255)
                                                     aGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAAc                                                                                             >  2:3297096/1‑148 (MQ=255)
                                                      gCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACg                                                                                            >  2:2851176/1‑148 (MQ=255)
                                                        tAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGtt                                                                                          <  2:948646/148‑1 (MQ=255)
                                                         aaCCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCTGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTg                                                                                         <  1:661078/148‑1 (MQ=255)
                                                         aaCCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTg                                                                                         >  1:62440/1‑148 (MQ=255)
                                                         aaCCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTg                                                                                         <  2:1251190/148‑1 (MQ=255)
                                                          aCCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGc                                                                                        >  1:2775259/1‑148 (MQ=255)
                                                          aCCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGc                                                                                        >  1:1698058/1‑148 (MQ=255)
                                                           cccGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCa                                                                                       >  1:2970368/1‑148 (MQ=255)
                                                           cccGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCa                                                                                       >  2:193811/1‑148 (MQ=255)
                                                                ggCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCga                                                                                  <  1:2161555/148‑1 (MQ=255)
                                                                  cagcagTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCgaag                                                                                >  2:1002264/1‑148 (MQ=255)
                                                                    gcagTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCgaagaa                                                                              >  1:2686013/1‑148 (MQ=255)
                                                                     cagTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGaagaag                                                                             >  2:1628088/1‑148 (MQ=255)
                                                                     cagTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGaagaag                                                                             >  1:939444/1‑148 (MQ=255)
                                                                      agTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGc                                                                            <  1:2563570/148‑1 (MQ=255)
                                                                            cgGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCg                                                                      <  1:1020612/148‑1 (MQ=255)
                                                                                aaCCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCTGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGaa                                                                  >  1:1994102/1‑148 (MQ=255)
                                                                                aaCCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGaa                                                                  >  1:1687720/1‑148 (MQ=255)
                                                                                  cccACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATa                                                                <  1:2198992/148‑1 (MQ=255)
                                                                                       aacaaGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACaaaa                                                           <  1:1532957/148‑1 (MQ=255)
                                                                                        acaaGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACaaaaa                                                          <  1:3092826/148‑1 (MQ=255)
                                                                                         caaGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGACTACaaaaaa                                                         >  2:1968804/1‑148 (MQ=255)
                                                                                         caaGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACaaaaaa                                                         <  2:724654/148‑1 (MQ=255)
                                                                                           aGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAgg                                                       <  1:988941/148‑1 (MQ=255)
                                                                                            gggCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGGc                                                      <  1:345631/148‑1 (MQ=255)
                                                                                             ggCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGGcg                                                     <  1:2759593/148‑1 (MQ=255)
                                                                                                  ccGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGGCGACGaa                                                >  1:119746/1‑148 (MQ=255)
                                                                                                    gtgtTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGGCGACGAAAt                                              <  2:1111281/148‑1 (MQ=255)
                                                                                                      gtTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCTGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGGCGACGAAATCg                                            >  1:2955252/1‑148 (MQ=255)
                                                                                                      gtTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGGCGACGAAATCg                                            <  2:2508920/148‑1 (MQ=255)
                                                                                                       ttGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGGCGACGAAATCGc                                           >  2:2977377/1‑148 (MQ=255)
                                                                                                       ttGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGGCGACGAAATCGc                                           >  2:1328759/1‑148 (MQ=255)
                                                                                                            ggTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGGCGACGAAATCGCTGCAg                                      <  2:2514645/148‑1 (MQ=255)
                                                                                                            ggTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGGCGACGAAATCGCTGCAg                                      <  2:2513569/148‑1 (MQ=255)
                                                                                                             gTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGGCGACGAAATCGCTGCAgt                                     <  2:2070615/148‑1 (MQ=255)
                                                                                                             gTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGGCGACGAAATCGCTGCAgt                                     <  1:570530/148‑1 (MQ=255)
                                                                                                             gTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGGCGACGAAATCGCTGCAgt                                     >  2:2951704/1‑148 (MQ=255)
                                                                                                               aaaaTCAAGTCTATCACTGACTTCGGTATCTTCATCTGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCCGTTCGTGAATACAGAAAAGGCGACGAAATCGCTGCAgttg                                   <  1:1549506/148‑1 (MQ=255)
                                                                                                                 aaTCAAGTCTATCACTGACTTCGGTATCTTCATCTGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGGCGACGAAATCGCTGCAgttgtt                                 <  2:3085372/148‑1 (MQ=255)
                                                                                                                 aaTCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGGCGACGAAATCGCTGCAgttgtt                                 >  1:2054420/1‑148 (MQ=255)
                                                                                                                 aaTCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGGCGACGAAATCGCTGCAgttgtt                                 <  2:141323/148‑1 (MQ=255)
                                                                                                                 aaTCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGGCGACGAAATCGCTGCAgttgtt                                 >  1:3170227/1‑148 (MQ=255)
                                                                                                                  aTCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGGCGACGAAATCGCTGCAGTTGTTc                                >  1:652150/1‑148 (MQ=255)
                                                                                                                    cAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACACAAAAGGCGACGAAATCGCTGCAGTTGTTCTg                              >  2:155600/1‑148 (MQ=255)
                                                                                                                    cAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGGCGACGAAATCGCTGCAGTTGTTCTg                              >  2:2489470/1‑148 (MQ=255)
                                                                                                                     aaGTCTATCACTGACTTCGGTATCTTCATCTGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGGCGACGAAATCGCTGCAGTTGTTCTGc                             >  2:223669/1‑148 (MQ=255)
                                                                                                                      aGTCTATCACTGACTTCGGTATCTTCATCTGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGGCGACGAAATCGCTGCAGTTGTTCTGCa                            <  1:2492468/148‑1 (MQ=255)
                                                                                                                      aGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGGCGACGAAATCGCTGCAGTTGTTCTGCa                            >  1:1331024/1‑148 (MQ=255)
                                                                                                                        tCTATCACTGACTTCGGTATCTTCATCTGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGGCGACGAAATCGCTGCAGTTGTTCTGCAgg                          <  2:502615/148‑1 (MQ=255)
                                                                                                                        tCTATCACTGACTTCGGTATCTTCATCTGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGGCGACGAAATCGCTGCAGTTGTTCTGCAgg                          <  2:2806739/148‑1 (MQ=255)
                                                                                                                        tCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGGCGACGAAATCGCTGCAGTTGTTCTGCAgg                          >  1:313534/1‑148 (MQ=255)
                                                                                                                         cTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGGCGACGAAATCGCTGCAGTTGTTCTGCAGGt                         <  2:506820/148‑1 (MQ=255)
                                                                                                                           aTCACTGACTTCGGTATCTTCATCTGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGGCGACGAAATCGCTGCAGTTGTTCTGCAGGTTg                       >  1:1211993/1‑148 (MQ=255)
                                                                                                                            tCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGGCGACGAAATCGCTGCAGTTGTTCTGCAGGTTGa                      <  1:3309776/148‑1 (MQ=255)
                                                                                                                             cACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGGCGACGAAATCGCTGCAGTTGTTCTGCAGGTTGAc                     >  1:1352039/1‑148 (MQ=255)
                                                                                                                             cACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGGCGACGAAATCGCTGCAGTTGTTCTGCAGGTTGAc                     >  2:1861698/1‑148 (MQ=255)
                                                                                                                                tGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGGCGACGAGATCGCTGCAGTTGTTCTGCAGGTTGACGCa                  >  2:939187/1‑148 (MQ=255)
                                                                                                                                  aCTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGGCGACGAACTCGCTGCAGTTGTTCTGCAGGTTGACGCAga                >  2:2949564/1‑148 (MQ=255)
                                                                                                                                     tCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGGCGACGAAATCGCTGCAGTTGTTCTGCAGGTTGACGCAgaacg             >  2:3136043/1‑148 (MQ=255)
                                                                                                                                         tATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGGCGACGAAATCGCTGCAGTTGTTCTGCAGGTTGACGCAgaacgtgaa         >  1:1987942/1‑148 (MQ=255)
                                                                                                                                          aTCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGGCGACGAAATCGCTGCAGTTGTTCTGCAGGTTGACGCAgaacgtgaac        <  2:240963/148‑1 (MQ=255)
                                                                                                                                           tCTTCATCTGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGGCGACGAAATCGCTGCGGTTGTTCTGCAGGTTGACGCAgaacgtgaacg       <  2:900749/148‑1 (MQ=255)
                                                                                                                                           tCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGGCGACGAAATCGCTGCAGTTGTTCTGCAGGTTGACGCAgaacgtgaacg       >  2:2285677/1‑148 (MQ=255)
                                                                                                                                           tCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGGCGACGAAATCGCTGCAGTTGTTCTGCAGGTTGACGCAgaacgtgaacg       <  2:2702173/148‑1 (MQ=255)
                                                                                                                                            cTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGGCGACGAAATCGCTGCAGTTGTTCTGCAGGTTGACGCAgaacgtgaacgt      <  1:183755/148‑1 (MQ=255)
                                                                                                                                            cTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGGCGACGAAATCGCTGCAGTTGTTCTGCAGGTTGACGCAgaacgtgaacgt      <  1:183557/148‑1 (MQ=255)
                                                                                                                                            cTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGGCGACGAAATCGCTGCAGTTGTTCTGCAGGTTGACGCAgaacgtgaacgt      <  2:2377549/148‑1 (MQ=255)
                                                                                                                                              tCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGGCGACGAAATCGCTGCAGTTGTTCTGCAGGTTGACGCAGAACGTGAACGTAt    <  1:193811/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
CTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGGCGACGAAATCGCTGCAGTTGTTCTGCAGGTTGACGCAGAACGTGAACGTATCT  >  CP009273/958447‑958738

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: