Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,218,928 G→A 8.2% E424K (GAA→AAA)  yjbB → putative Na+/Pi‑cotransporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,218,9280GA 8.2% 719.1 / 34.9 208E424K (GAA→AAA) yjbBputative Na+/Pi‑cotransporter
Reads supporting (aligned to +/- strand):  ref base G (98/93);  new base A (7/10);  total (105/103)
Fisher's exact test for biased strand distribution p-value = 4.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 210 total aligned reads displayed.

ACGTTCTCTACACCGCCATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTT  >  CP009273/4218782‑4219075
                                                                                                                                                  |                                                                                                                                                   
aCGTTCTCTACACCGCCATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGa                                                                                                                                                    >  2:883294/1‑148 (MQ=255)
aCGTTCTCTACACCGCCATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGa                                                                                                                                                    >  2:1714369/1‑148 (MQ=255)
aCGTTCTCTACACCGCCATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCAGGGGCAGCGa                                                                                                                                                    >  2:3093842/1‑148 (MQ=255)
 cGTTCTCTACACCGCCATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGaa                                                                                                                                                   <  2:288501/148‑1 (MQ=255)
 cGTTCTCTACACCGCCATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGaa                                                                                                                                                   <  1:2209721/148‑1 (MQ=255)
 cGTTCTCTACACCGCCATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGaa                                                                                                                                                   >  2:2046598/1‑148 (MQ=255)
  gTTCTATACACCGCCATTAAGCTTTATCTGGCGCGGATGCCAAAAGAAGAGCTGTCGGAGGAAGAGTCGCGCCTCGGGGCAGAGTTCATCGAAATGTCGCTCAAACTTGAACAGGCCTCCGATATCGTAGAGCGCATGGGCAGCGaaa                                                                                                                                                  <  2:449102/148‑1 (MQ=255)
    tctcTACACCGCCATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAAtt                                                                                                                                                >  1:126058/1‑148 (MQ=255)
     ctctACACCGCCATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTg                                                                                                                                               <  1:2842147/148‑1 (MQ=255)
      tctACACCGCCATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGc                                                                                                                                              <  1:1718939/148‑1 (MQ=255)
       ctACACCGCCATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCt                                                                                                                                             <  1:3225885/148‑1 (MQ=255)
       ctACACCGCCATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCt                                                                                                                                             >  1:690199/1‑148 (MQ=255)
          cacCGCCATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGAc                                                                                                                                          >  1:2532001/1‑148 (MQ=255)
          cacCGCCATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGAc                                                                                                                                          >  1:2182820/1‑148 (MQ=255)
             cGCCATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACaaa                                                                                                                                       >  2:2298238/1‑148 (MQ=255)
               ccATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGCGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATc                                                                                                                                     >  2:2344016/1‑148 (MQ=255)
                cATTAAGCTGTATCTGGCGCGGATTCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTTGTCAGAGATCATCGAAATGTCGCTCAAACTTGAACAGGCCTACGATATCGTCGAGCGGATGGGCATCGAAATTGCTGACAAATCa                                                                                                                                    <  2:1054588/148‑1 (MQ=255)
                 aTTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGAGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCAc                                                                                                                                   <  1:300275/148‑1 (MQ=255)
                 aTTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGAGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCAc                                                                                                                                   <  1:2942443/148‑1 (MQ=255)
                  ttAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACt                                                                                                                                  >  1:860173/1‑148 (MQ=255)
                  ttAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACt                                                                                                                                  <  2:1590297/148‑1 (MQ=255)
                    aaGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTgg                                                                                                                                <  1:1240094/148‑1 (MQ=255)
                      gCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGgca                                                                                                                              >  2:2953805/1‑148 (MQ=255)
                        tGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGgcagc                                                                                                                            >  2:2409019/1‑148 (MQ=255)
                            tCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCAcgg                                                                                                                        <  1:1058/148‑1 (MQ=255)
                             cTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCAcggc                                                                                                                       >  2:3179508/1‑148 (MQ=255)
                                  gcgGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCa                                                                                                                  <  1:949015/148‑1 (MQ=255)
                                   cgGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCAt                                                                                                                 >  1:1580698/1‑148 (MQ=255)
                                      aTGCCCAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCAtttt                                                                                                              >  2:916017/1‑148 (MQ=255)
                                      aTGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCAtttt                                                                                                              <  2:725846/148‑1 (MQ=255)
                                         ccAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGc                                                                                                           <  2:302530/148‑1 (MQ=255)
                                         ccAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGc                                                                                                           >  2:125749/1‑148 (MQ=255)
                                          cAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCt                                                                                                          <  1:2460877/148‑1 (MQ=255)
                                           aaaaGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCtt                                                                                                         >  1:37386/1‑148 (MQ=255)
                                              agaagaGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGAt                                                                                                      <  2:2196045/148‑1 (MQ=255)
                                               gaagaGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATg                                                                                                     >  2:1661869/1‑148 (MQ=255)
                                                aagaGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATgg                                                                                                    >  2:757742/1‑148 (MQ=255)
                                                 agagCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATggg                                                                                                   <  1:187697/148‑1 (MQ=255)
                                                 agagCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATggg                                                                                                   <  2:1242155/148‑1 (MQ=255)
                                                   agCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGtt                                                                                                 >  1:1044671/1‑148 (MQ=255)
                                                   agCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGtt                                                                                                 >  1:925485/1‑148 (MQ=255)
                                                   agCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGtt                                                                                                 >  2:1343362/1‑148 (MQ=255)
                                                       ggCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAg                                                                                             >  2:1413575/1‑148 (MQ=255)
                                                        gCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAgg                                                                                            >  1:1778412/1‑148 (MQ=255)
                                                         cGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGa                                                                                           <  1:330474/148‑1 (MQ=255)
                                                          ggaggaAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGAGATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGaa                                                                                          <  2:1572802/148‑1 (MQ=255)
                                                              gaAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTgg                                                                                      <  1:1490921/148‑1 (MQ=255)
                                                                 gagTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATg                                                                                   >  1:479434/1‑148 (MQ=255)
                                                                 gagTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATg                                                                                   >  1:591231/1‑148 (MQ=255)
                                                                    tCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGACCTGGATgcgc                                                                                >  2:505727/1‑148 (MQ=255)
                                                                    tCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATgcgc                                                                                <  1:1567494/148‑1 (MQ=255)
                                                                          cgcTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATg                                                                          >  1:1528627/1‑148 (MQ=255)
                                                                          cgcTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATg                                                                          <  2:129487/148‑1 (MQ=255)
                                                                               ggCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCaa                                                                     >  1:3067597/1‑148 (MQ=255)
                                                                               ggCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCaa                                                                     >  1:3068320/1‑148 (MQ=255)
                                                                                   gagaTCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGc                                                                 >  1:1634462/1‑148 (MQ=255)
                                                                                   gagaTCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGc                                                                 >  2:1219194/1‑148 (MQ=255)
                                                                                   gagaTCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGc                                                                 >  1:1314305/1‑148 (MQ=255)
                                                                                   gagaTCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGc                                                                 <  2:1268395/148‑1 (MQ=255)
                                                                                    agaTCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCt                                                                >  1:2698154/1‑148 (MQ=255)
                                                                                    agaTCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCt                                                                >  2:769032/1‑148 (MQ=255)
                                                                                     gaTCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTc                                                               <  2:2398511/148‑1 (MQ=255)
                                                                                     gaTCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTc                                                               <  2:739578/148‑1 (MQ=255)
                                                                                       tcatcGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAg                                                             >  1:913924/1‑148 (MQ=255)
                                                                                          tcGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTaa                                                          <  2:1528635/148‑1 (MQ=255)
                                                                                           cGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAAt                                                         <  2:599651/148‑1 (MQ=255)
                                                                                           cGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAAt                                                         <  1:785810/148‑1 (MQ=255)
                                                                                           cGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAAt                                                         >  1:2899767/1‑148 (MQ=255)
                                                                                             aaaTGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAAttt                                                       >  1:2209965/1‑148 (MQ=255)
                                                                                             aaaTGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAAttt                                                       >  1:1198975/1‑148 (MQ=255)
                                                                                               aTGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTaa                                                     >  2:1967759/1‑148 (MQ=255)
                                                                                                    gCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTg                                                <  2:1778413/148‑1 (MQ=255)
                                                                                                         aCCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAAt                                           >  2:3063643/1‑148 (MQ=255)
                                                                                                         aCCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAAt                                           >  1:2362896/1‑148 (MQ=255)
                                                                                                         aCCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAAt                                           >  2:2454832/1‑148 (MQ=255)
                                                                                                         aCCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAAt                                           <  1:813393/148‑1 (MQ=255)
                                                                                                          ccTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATg                                          >  1:1579014/1‑148 (MQ=255)
                                                                                                          ccTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATg                                          <  1:304809/148‑1 (MQ=255)
                                                                                                          ccTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATg                                          <  2:1889873/148‑1 (MQ=255)
                                                                                                           cTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGt                                         <  2:2596551/148‑1 (MQ=255)
                                                                                                             tGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCg                                       >  2:2438175/1‑148 (MQ=255)
                                                                                                              gAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCgg                                      <  2:394120/148‑1 (MQ=255)
                                                                                                                 cAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGgtgt                                   >  1:812822/1‑148 (MQ=255)
                                                                                                                     ccTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGttctt                               <  1:2326750/148‑1 (MQ=255)
                                                                                                                     ccTCCGAGATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGttctt                               <  1:667534/148‑1 (MQ=255)
                                                                                                                       tCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTtcttct                             <  1:1152736/148‑1 (MQ=255)
                                                                                                                        ccGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTAGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTtctc                            <  1:431585/148‑1 (MQ=255)
                                                                                                                         cGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTctct                           <  1:2171677/148‑1 (MQ=255)
                                                                                                                            tatCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGc                        >  1:1653978/1‑148 (MQ=255)
                                                                                                                             atCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCg                       <  1:2041391/148‑1 (MQ=255)
                                                                                                                                   gAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCa                 >  1:1181967/1‑148 (MQ=255)
                                                                                                                                    aGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCAc                >  2:1142274/1‑148 (MQ=255)
                                                                                                                                    aGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCAc                <  2:1799908/148‑1 (MQ=255)
                                                                                                                                        cATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAgc            >  1:425973/1‑148 (MQ=255)
                                                                                                                                         aTGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAgcg           <  2:208059/148‑1 (MQ=255)
                                                                                                                                                gCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCtcgtcg    <  1:1294265/148‑1 (MQ=255)
                                                                                                                                                gCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCtcgtcg    <  1:424553/148‑1 (MQ=255)
                                                                                                                                                 cGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTcgtcgt   >  2:551981/1‑148 (MQ=255)
                                                                                                                                                 cGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTcgtcgt   <  2:2297641/148‑1 (MQ=255)
                                                                                                                                                  gAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGtt  <  1:3179508/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
ACGTTCTCTACACCGCCATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTT  >  CP009273/4218782‑4219075

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: