Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 172,107 G→T 9.6% E172* (GAA→TAA)  clcA → H(+)/Cl(‑) exchange transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273172,1070GT 9.6% 635.2 / 26.9 181E172* (GAA→TAA) clcAH(+)/Cl(‑) exchange transporter
Reads supporting (aligned to +/- strand):  ref base G (76/87);  new base T (4/14);  total (80/101)
Fisher's exact test for biased strand distribution p-value = 7.77e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 181 total aligned reads displayed.

CTGGTGGCGTGTATTGCCGGTGAAGTTCTTTGGCGGGCTGGGGACACTCGGCGGAGGCATGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATACGTTAATTTCGATTAAAGCGGTATTTATT  >  CP009273/171962‑172253
                                                                                                                                                 |                                                                                                                                                  
cTGGTGGCGTGTATTGCCGGTGAAGTTCTTTGGCGGGCTGGGGACACTCGGCGGAGGCATGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGaa                                                                                                                                                  <  2:3283871/148‑1 (MQ=255)
 tggtggCGTGTATTGCCGGTGAAGTTCTTTGGCGGGCTGGGGACACTCGGCGGAGGCATGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACTAAg                                                                                                                                                 <  1:1000896/148‑1 (MQ=255)
 tggtggCGTGTATTGCCGGTGAAGTTCTTTGGCGGGCTGGGGACACTCGGCGGAGGCATGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAg                                                                                                                                                 <  1:2587049/148‑1 (MQ=255)
     ggCGTGTATTGCCGGTGAAGTTCTTTGGCGGGCTGGGGACACTCGGCGGAGGCATGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCg                                                                                                                                             >  2:983707/1‑148 (MQ=255)
     ggCGTGTATTGCCGGTGAAGTTCTTTGGCGGGCTGGGGACACTCGGCGGAGGCATGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCg                                                                                                                                             >  2:557275/1‑148 (MQ=255)
        gtgtATTGCCGGTGAAGTTCTTTGGCGGGCTGGGGACACTCGGCGGAGGCATGGTGTTGGGGCGCGAAGGGCCACCCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCa                                                                                                                                          <  1:429687/148‑1 (MQ=255)
         tgtATTGCCGGTGAAGTTCTTTTGCGGGCTGGGGACACTCGGCGGAGGCATGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCAt                                                                                                                                         >  1:3100514/1‑148 (MQ=255)
         tgtATTGCCGGTGAAGTTCTTTGGCGGGCTGGGGACACTCGGCGGAGGCATGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCAt                                                                                                                                         >  2:1582038/1‑148 (MQ=255)
          gtATTGCCGGTGAAGTTCTTTGGCGGGCTGGGGACACTCGGCGGAGGCATGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATa                                                                                                                                        >  1:1044593/1‑148 (MQ=255)
            aTTGCCGGTGAAGTTCTTTGGCGGGCTGGGGACACTCGGCGGAGGCATGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACTAAGCTCGCCATACg                                                                                                                                      <  2:236462/148‑1 (MQ=255)
            aTTGCCGGTGAAGTTCTTTGGCGGGCTGGGGACACTCGGCGGAGGCATGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACg                                                                                                                                      >  1:1313699/1‑148 (MQ=255)
              tGCCGGTGAAGTTCTTTGGCGGGCTGGGGACACTCGGCGGAGGCATGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACgct                                                                                                                                    >  1:790579/1‑148 (MQ=255)
                   gTGAAGTTCTTTGGCGGGCTGGGGACACTCGGCGGAGGCATGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTgg                                                                                                                               >  2:1476691/1‑148 (MQ=255)
                    tGAAGTTCTTTGGCGGGCTGGGGACACTCGGCGGAGGCATGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGc                                                                                                                              >  1:2892998/1‑148 (MQ=255)
                      aaGGTCTTGGGCGGGCTGGGGACACTCGGCGGAGGCATGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCaa                                                                                                                            <  2:1812469/148‑1 (MQ=255)
                         ttCTTTGGCGGGCTGGGGACACTCGGCGGAGGCATGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCg                                                                                                                         >  2:1976287/1‑148 (MQ=255)
                         ttCTTTGGCGGGCTGGGGACACTCGGCGGAGGCATGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCg                                                                                                                         >  1:998481/1‑148 (MQ=255)
                         ttCTTTGGCGGGCTGGGGACACTCGGCGGAGGCATGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCg                                                                                                                         >  1:617445/1‑148 (MQ=255)
                         ttCTTTGGCGGGCTGGGGACACTCGGCGGAGGCATGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCg                                                                                                                         >  1:2653378/1‑148 (MQ=255)
                            tttGGCGGGCTGGGGACACTCGGCGGAGGCATGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACTAAGCTCGCCATACGCTGCTGGCAACCGGtg                                                                                                                      <  2:2985123/148‑1 (MQ=255)
                               ggCGGGCTGGGGACACTCGGCGGAGGCATGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGtgctg                                                                                                                   <  1:3083784/148‑1 (MQ=255)
                                  gggCTGGGGACACTCGGCGGAGGCATGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGtgctgctg                                                                                                                <  1:15020/148‑1 (MQ=255)
                                   ggCTGGGGACACTCGGCGGAGGCATGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGtgctgctgc                                                                                                               >  1:1558007/1‑148 (MQ=255)
                                        gggACACTCGGCGGAGGCATGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGggc                                                                                                          >  1:321746/1‑148 (MQ=255)
                                          gACACTCGGCGGAGGCATGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGggctg                                                                                                        <  2:1090410/148‑1 (MQ=255)
                                                    ggAGGCATGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCt                                                                                              <  1:2615791/148‑1 (MQ=255)
                                                       ggCATGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTa                                                                                           >  1:2358979/1‑148 (MQ=255)
                                                       ggCATGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTa                                                                                           <  1:659468/148‑1 (MQ=255)
                                                        gCATGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTaa                                                                                          >  1:1255116/1‑148 (MQ=255)
                                                            ggTGTTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAAcgcg                                                                                      <  2:95828/148‑1 (MQ=255)
                                                             gtgtTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACgcgc                                                                                     <  2:783947/148‑1 (MQ=255)
                                                               gtTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGcgccg                                                                                   <  1:936286/148‑1 (MQ=255)
                                                                   gggCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTgg                                                                               <  1:2540660/148‑1 (MQ=255)
                                                                     gcgcGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCg                                                                             >  2:1957989/1‑148 (MQ=255)
                                                                         gAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTa                                                                         >  2:15844/1‑148 (MQ=255)
                                                                         gAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTa                                                                         <  1:2176932/148‑1 (MQ=255)
                                                                          aaGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTAt                                                                        <  1:983707/148‑1 (MQ=255)
                                                                          aaGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTAt                                                                        <  1:557275/148‑1 (MQ=255)
                                                                          aaGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTAt                                                                        <  2:1930962/148‑1 (MQ=255)
                                                                           aGGGCCAACCGTGCAGATCGGCGGTAACATTGTCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACTAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTAtt                                                                       <  1:316030/148‑1 (MQ=255)
                                                                            gggCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTAttt                                                                      <  1:2332951/148‑1 (MQ=255)
                                                                             cgCCAACCGTGCAGAACGGCGCTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTAtttt                                                                     <  2:2964036/147‑1 (MQ=255)
                                                                              gCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTg                                                                    <  1:1368836/148‑1 (MQ=255)
                                                                              gCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTg                                                                    <  1:780975/148‑1 (MQ=255)
                                                                                 aaCCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGttt                                                                 >  1:1884231/1‑148 (MQ=255)
                                                                                  aCCGTGCAGATCGGCGGTAACCTTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACTAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTtta                                                                <  2:1791152/148‑1 (MQ=255)
                                                                                   ccGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTttat                                                               >  1:2555104/1‑148 (MQ=255)
                                                                                     gTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTttatta                                                             <  2:2630506/148‑1 (MQ=255)
                                                                                      tGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTtattat                                                            <  2:997831/148‑1 (MQ=255)
                                                                                      tGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTtattat                                                            <  1:2046414/148‑1 (MQ=255)
                                                                                         aGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGa                                                         <  2:2772144/148‑1 (MQ=255)
                                                                                            tCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAaga                                                      >  1:1496463/1‑148 (MQ=255)
                                                                                             cggcggTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAagag                                                     <  2:822673/148‑1 (MQ=255)
                                                                                             cggcggTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAagag                                                     <  2:306110/148‑1 (MQ=255)
                                                                                               gcggTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGAt                                                   <  2:38824/148‑1 (MQ=255)
                                                                                                   tAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGt                                               >  2:3217664/1‑148 (MQ=255)
                                                                                                       aTTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGc                                           <  1:3248173/148‑1 (MQ=255)
                                                                                                        ttGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCa                                          >  2:308388/1‑148 (MQ=255)
                                                                                                         tGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAg                                         >  1:181019/1‑148 (MQ=255)
                                                                                                         tGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAg                                         >  2:1649703/1‑148 (MQ=255)
                                                                                                         tGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAg                                         >  2:2717916/1‑148 (MQ=255)
                                                                                                          ggCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGt                                        <  1:379165/148‑1 (MQ=255)
                                                                                                           gCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGtt                                       <  2:2648453/148‑1 (MQ=255)
                                                                                                           gCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGtt                                       <  1:1362607/148‑1 (MQ=255)
                                                                                                           gCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGtt                                       <  2:254719/148‑1 (MQ=255)
                                                                                                           gCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGtt                                       <  1:479097/148‑1 (MQ=255)
                                                                                                            ccGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGttt                                      <  2:271093/148‑1 (MQ=255)
                                                                                                             cGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTc                                     >  2:348534/1‑148 (MQ=255)
                                                                                                              gTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCg                                    >  2:938791/1‑148 (MQ=255)
                                                                                                                aTGGTGCTTGATATTTTCCGCCTGAAAGGTGACTAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCt                                  >  1:1548589/1‑148 (MQ=255)
                                                                                                                 tGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCta                                 <  2:198470/148‑1 (MQ=255)
                                                                                                                  ggTGCTTGATATTTTCCGCCTGAAAGGTGACTAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCtat                                <  2:2079415/148‑1 (MQ=255)
                                                                                                                  ggTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCtat                                <  2:1146772/148‑1 (MQ=255)
                                                                                                                  ggTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCtat                                <  2:471755/148‑1 (MQ=255)
                                                                                                                     gCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATACg                             <  1:1975943/148‑1 (MQ=255)
                                                                                                                        tGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATACGTTa                          <  1:1723753/148‑1 (MQ=255)
                                                                                                                        tGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATACGTTa                          >  2:881344/1‑148 (MQ=255)
                                                                                                                         gATATTTTCCGCCTGAAAGGTGACTAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATACGTTaa                         <  2:2401663/148‑1 (MQ=255)
                                                                                                                         gATATTTTCCGCCTGAAAGGTGACTAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATACGTTaa                         <  2:742826/148‑1 (MQ=255)
                                                                                                                         gATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATACGTTaa                         <  1:1428233/148‑1 (MQ=255)
                                                                                                                          atatTTTCCGCCTGAAAGGTGACTAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATACGTTAAt                        <  2:1411647/148‑1 (MQ=255)
                                                                                                                          atatTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATACGTTAAt                        <  1:3057784/148‑1 (MQ=255)
                                                                                                                            atTTTCCGCCTGAAAGGTGACTAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATACGTTAAttt                      >  1:149140/1‑148 (MQ=255)
                                                                                                                            atTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATACGTTAAttt                      <  2:36044/148‑1 (MQ=255)
                                                                                                                             ttttCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATACGTTAATTTc                     >  2:641418/1‑148 (MQ=255)
                                                                                                                             ttttCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATACGTTAATTTc                     >  2:641489/1‑148 (MQ=255)
                                                                                                                             ttttCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATACGTTAATTTc                     >  1:2079669/1‑148 (MQ=255)
                                                                                                                              tttCCGCCTGAAAGGTGACTAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATACGTTAATTTCg                    <  1:115060/148‑1 (MQ=255)
                                                                                                                              tttCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATACGTTAATTTCg                    <  2:2447258/148‑1 (MQ=255)
                                                                                                                                tCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATACGTTAATTTCGAt                  <  2:1313699/148‑1 (MQ=255)
                                                                                                                                      tGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATACGTTAATTTCGATTAAAGc            <  2:1542533/148‑1 (MQ=255)
                                                                                                                                       gAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATACGTTAATTTCGATTAAAGCg           >  2:884905/1‑148 (MQ=255)
                                                                                                                                       gAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGACGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATACGTTAATTTCGATTAAAGCg           >  1:3281172/1‑148 (MQ=255)
                                                                                                                                        aaaGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATACGTTAATTTCGATTAAAGCgg          <  2:2103768/148‑1 (MQ=255)
                                                                                                                                        aaaGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATACGTTAATTTCGATTAAAGCgg          >  1:2031426/1‑148 (MQ=255)
                                                                                                                                           ggTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATACGTTAATTTCGATTAAAGCGGtat       <  2:1255116/148‑1 (MQ=255)
                                                                                                                                             tGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATACGTTAATTTCGATTAAAGCGGtattt     <  1:348534/148‑1 (MQ=255)
                                                                                                                                              gACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATACGTTAATTTCGATTAAAGCGGtattta    <  2:1884231/148‑1 (MQ=255)
                                                                                                                                              gACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATACGTTAATTTCGATTAAAGCGGtattta    >  1:2948866/1‑148 (MQ=255)
                                                                                                                                                cGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATACGTTAATTTCGATTAAAGCGGtatttatt  <  1:849043/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
CTGGTGGCGTGTATTGCCGGTGAAGTTCTTTGGCGGGCTGGGGACACTCGGCGGAGGCATGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATACGTTAATTTCGATTAAAGCGGTATTTATT  >  CP009273/171962‑172253

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: ‑