Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 172,107 G→T 5.4% 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 5.4% 631.9 / 9.4 166E172* (GAA→TAA) clcAH(+)/Cl(‑) exchange transporter
Reads supporting (aligned to +/- strand):  ref base G (87/70);  new base T (4/5);  total (91/75)
Fisher's exact test for biased strand distribution p-value = 7.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 168 total aligned reads displayed.

CGCTGGTGGCGTGTATTGCCGGTGAAGTTCTTTGGCGGGCTGGGGACACTCGGCGGAGGCATGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATACGTTAATTTCGATTAAAGCGGTATTTATT  >  CP009273/171960‑172253
                                                                                                                                                   |                                                                                                                                                  
cGCTGGTGGCGTGTATTGCCGGTGAAGTTCTTTGGCGGGCTGGGGACACTCGGCGGAGGCATGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACg                                                                                                                                                    >  2:1777820/1‑148 (MQ=255)
cGCTGGTGGCGTGTATTGCCGGTGAAGTTCTTTGGCGGGCTGGGGACACTCGGCGGAGGCATGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACg                                                                                                                                                    >  2:1778123/1‑148 (MQ=255)
  cTGGTGGCGTGTATTGCCGGTGAAGTTCTTTGGCGGGCTGGGGACACTCGGCGGAGGCATGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGaa                                                                                                                                                  >  1:2892149/1‑148 (MQ=255)
     gtggCGTGTATTGCCGGTGAAGTTCTTTGGCGGGCTGGGGACACTCGGCGGAGGCATGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCt                                                                                                                                               >  1:1007139/1‑148 (MQ=255)
     gtggCGTGTATTGCCGGTGAAGTTCTTTGGCGGGCTGGGGACACTCGGCGGAGGCATGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTAGGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCt                                                                                                                                               >  1:3038869/1‑148 (MQ=255)
      tggCGTGTATTGCCGGTGAAGTTCTTTGGCGGGCTGGGGACACTCGGCGGAGGCATGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTc                                                                                                                                              >  1:1948542/1‑148 (MQ=255)
            gtATTGCCGGTGAAGTTCTTTGGCGGGCTGGGGACACTCGGCGGAGGCATGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATa                                                                                                                                        <  2:138374/148‑1 (MQ=255)
              aTTGCCGGTGAAGTTCTTTGGCGGGCTGGGGACACTCGGCGGAGGCATGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACg                                                                                                                                      >  2:567536/1‑148 (MQ=255)
               ttGCCGGTGAAGTTCTTTGGCGGGCTGGGGACACTCGGCGGAGGCATGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTCTGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACgc                                                                                                                                     >  1:2156973/1‑148 (MQ=255)
               ttGCCGGTGAAGTTCTTTGGCGGGCTGGGGACACTCGGCGGAGGCATGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACgc                                                                                                                                     >  2:2438843/1‑148 (MQ=255)
                tGCCGGTGAAGTTCTTTGGCGGGCTGGGGACACTCGGCGGAGGCATGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACgct                                                                                                                                    >  1:2481048/1‑148 (MQ=255)
                         aGTTCTTTGGCGGGCTGGGGACACTCGGCGGAGGCATGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAAc                                                                                                                           >  1:829215/1‑148 (MQ=255)
                          gTTCTTTGGCGGGCTGGGGACACTCGGCGGAGGCATGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAAcc                                                                                                                          >  1:727152/1‑148 (MQ=255)
                           ttCTTTGGCGGGCTGGGGACACTCGGCGGAGGCATGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCg                                                                                                                         >  2:1976112/1‑148 (MQ=255)
                           ttCTTTGGCGGGCTGGGGACACTCGGCGGAGGCATGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCg                                                                                                                         >  2:3136907/1‑148 (MQ=255)
                           ttCTTTGGCGGGCTGGGGACACTCGGCGGAGGCATGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCg                                                                                                                         >  1:2266229/1‑148 (MQ=255)
                            tCTTTGGCGGGCTGGGGACACTCGGCGGAGGCATGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCgg                                                                                                                        >  1:1157739/1‑148 (MQ=255)
                             cTTTGGCGGGCTGGGGACACTCGGCGGAGGCATGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGt                                                                                                                       >  2:918911/1‑148 (MQ=255)
                              tttGGCGGGCTGGGGACACTCGGCGGAGGCATGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGtg                                                                                                                      <  2:183469/148‑1 (MQ=255)
                               ttGGCGGGCTGGGGACACTCGGCGGAGGCATGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGtgc                                                                                                                     <  2:3268375/148‑1 (MQ=255)
                                     ggCTGGGGACACTCGGCGGAGGCCTGGGGAGGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGagctgctgc                                                                                                               >  1:2387075/1‑148 (MQ=255)
                                     ggCTGGGGACACTCGGCGGAGGCATGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGtgctgctgc                                                                                                               >  1:993529/1‑148 (MQ=255)
                                      gCTGGGGACACTCGGCGGAGGCATGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCg                                                                                                              >  2:1390530/1‑148 (MQ=255)
                                          gggACACTCGGCGGAGGCATGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGggc                                                                                                          <  1:2378119/148‑1 (MQ=255)
                                          gggACACTCGGCGGAGGCATGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGggc                                                                                                          <  1:2818634/148‑1 (MQ=255)
                                          gggACACTCGGCGGAGGCATGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGggc                                                                                                          >  1:2711213/1‑148 (MQ=255)
                                             acacTCGGCGGAGGCATGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGggctgg                                                                                                       <  2:28467/148‑1 (MQ=255)
                                               acTCGGCGGAGGCCTGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACTAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGggctggct                                                                                                     >  1:3338223/1‑148 (MQ=255)
                                                cTCGGCGGAGGCATGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGgctggctg                                                                                                    >  2:2294610/1‑148 (MQ=255)
                                                       gAGGCATGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCtt                                                                                             <  1:3218634/148‑1 (MQ=255)
                                                          gCATGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTaa                                                                                          <  1:318054/148‑1 (MQ=255)
                                                             tGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAAcgc                                                                                       <  1:3303472/148‑1 (MQ=255)
                                                               gtgtTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACTAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACgcgc                                                                                     >  1:154498/1‑148 (MQ=255)
                                                               gtgtTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACgcgc                                                                                     <  2:2432050/148‑1 (MQ=255)
                                                                     gggCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTgg                                                                               >  1:2607657/1‑148 (MQ=255)
                                                                       gcgcGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCg                                                                             <  1:2781115/148‑1 (MQ=255)
                                                                          cgAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGt                                                                          <  1:596544/148‑1 (MQ=255)
                                                                            aaGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTAt                                                                        <  1:3549244/148‑1 (MQ=255)
                                                                             aGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTAtt                                                                       <  1:1460090/148‑1 (MQ=255)
                                                                             aGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTAtt                                                                       <  2:2679999/148‑1 (MQ=255)
                                                                             aGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTAtt                                                                       <  1:2359004/148‑1 (MQ=255)
                                                                             aGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTAtt                                                                       <  2:2701314/148‑1 (MQ=255)
                                                                              gggCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTCttt                                                                      >  2:621243/1‑148 (MQ=255)
                                                                              gggCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTAttt                                                                      >  1:1194296/1‑148 (MQ=255)
                                                                              gggCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTAttt                                                                      >  2:1677115/1‑148 (MQ=255)
                                                                              gggCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTAttt                                                                      >  1:1194319/1‑148 (MQ=255)
                                                                              gggCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTAttt                                                                      >  1:2754511/1‑148 (MQ=255)
                                                                                gCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTg                                                                    <  2:2556911/148‑1 (MQ=255)
                                                                                    aCCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACTAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTtta                                                                <  2:1883540/148‑1 (MQ=255)
                                                                                     ccGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTttat                                                               <  1:325642/148‑1 (MQ=255)
                                                                                     ccGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTttat                                                               >  1:1592725/1‑148 (MQ=255)
                                                                                      cGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTttatt                                                              <  2:2921301/148‑1 (MQ=255)
                                                                                       gTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTttatta                                                             <  1:2343638/148‑1 (MQ=255)
                                                                                       gTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTttatta                                                             <  2:2757631/148‑1 (MQ=255)
                                                                                          cAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCg                                                          <  1:1861488/148‑1 (MQ=255)
                                                                                          cAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCg                                                          <  1:1859145/148‑1 (MQ=255)
                                                                                             aTCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAag                                                       >  2:2143827/1‑148 (MQ=255)
                                                                                             aTCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAag                                                       >  1:429082/1‑148 (MQ=255)
                                                                                             aTCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAag                                                       <  1:713505/148‑1 (MQ=255)
                                                                                             aTCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAag                                                       >  1:2057139/1‑148 (MQ=255)
                                                                                              tCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAaga                                                      >  1:630119/1‑148 (MQ=255)
                                                                                               cggcggTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAagag                                                     >  2:1201983/1‑148 (MQ=255)
                                                                                                ggcggTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAgaga                                                    <  1:3555869/148‑1 (MQ=255)
                                                                                                   ggTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGc                                                 <  1:3314872/148‑1 (MQ=255)
                                                                                                     tAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGt                                               <  2:1532634/148‑1 (MQ=255)
                                                                                                     tAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGt                                               <  2:3402582/148‑1 (MQ=255)
                                                                                                      aaCATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTc                                              >  2:2553954/1‑148 (MQ=255)
                                                                                                      aaCATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTc                                              <  1:791273/148‑1 (MQ=255)
                                                                                                         aTTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGc                                           >  1:351532/1‑148 (MQ=255)
                                                                                                           tGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACTAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAg                                         >  2:3411359/1‑148 (MQ=255)
                                                                                                                 tATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGc                                   >  1:993243/1‑148 (MQ=255)
                                                                                                                  aTGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCt                                  >  1:1741271/1‑148 (MQ=255)
                                                                                                                    ggTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCtat                                <  1:1685329/148‑1 (MQ=255)
                                                                                                                    ggTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCtat                                >  2:385652/1‑148 (MQ=255)
                                                                                                                        cTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATACGt                            <  2:1336166/148‑1 (MQ=255)
                                                                                                                           gATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATACGTTaa                         <  1:1390530/148‑1 (MQ=255)
                                                                                                                             tatTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTCTTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATACGTTAAtt                       >  2:3004228/1‑148 (MQ=255)
                                                                                                                              atTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATACGTTAAttt                      >  2:675547/1‑148 (MQ=255)
                                                                                                                                tttCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCTACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATACGTTAATTTCg                    <  2:1978794/148‑1 (MQ=255)
                                                                                                                                tttCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATACGTTAATTTCg                    >  1:1321814/1‑148 (MQ=255)
                                                                                                                                 ttCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATACGTTAATTTCGa                   <  1:3156971/148‑1 (MQ=255)
                                                                                                                                  tCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATACGTTAATTTCGAt                  <  2:2481048/148‑1 (MQ=255)
                                                                                                                                    cGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATACGTTAATTTCGATTa                <  2:1105066/148‑1 (MQ=255)
                                                                                                                                      ccTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATACGTTAATTTCGATTaaa              <  1:1976112/148‑1 (MQ=255)
                                                                                                                                      ccTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATACGTTAATTTCGATTaaa              >  2:2970026/1‑148 (MQ=255)
                                                                                                                                      ccTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATACGTTAATTTCGATTaaa              >  1:2777480/1‑148 (MQ=255)
                                                                                                                                        tGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATACGTTAATTTCGATTAAAGc            >  2:2210119/1‑148 (MQ=255)
                                                                                                                                        tGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATACGTTAATTTCGATTAAAGc            >  2:2210108/1‑148 (MQ=255)
                                                                                                                                         gAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATACGTTAATTTCGATTAAAGCg           >  1:2664386/1‑148 (MQ=255)
                                                                                                                                          aaaGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATACGTTAATTTCGATTAAAGCgg          >  1:3137578/1‑148 (MQ=255)
                                                                                                                                           aaGGTGACTAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATACGTTAATTTCGATTAAAGCGGt         <  1:1944139/148‑1 (MQ=255)
                                                                                                                                           aaGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATACGTTAATTTCGATTAAAGCGGt         >  2:3383465/1‑148 (MQ=255)
                                                                                                                                            aGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATACGTTAATTTCGATTAAAGCGGta        <  1:1048002/148‑1 (MQ=255)
                                                                                                                                            aGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATACGTTAATTTCGATTAAAGCGGta        <  1:3438410/148‑1 (MQ=255)
                                                                                                                                            aGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATACGTTAATTTCGATTAAAGCGGta        <  1:3438608/148‑1 (MQ=255)
                                                                                                                                             ggTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATACGTTAATTTCGATTAAAGCGGtat       >  1:218186/1‑148 (MQ=255)
                                                                                                                                               tGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATACGTTAATTTCGATTAAAGCGGtattt     <  1:918911/148‑1 (MQ=255)
                                                                                                                                                gACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATACGTTAATTTCGATTAAAGCGGtattta    >  1:3314294/1‑148 (MQ=255)
                                                                                                                                                 aCGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATACGTTAATTTCGATTAAAGCGGtatttat   >  1:3293895/1‑148 (MQ=255)
                                                                                                                                                  cGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATACGTTAATTTCGATTAAAGCGGtatttatt  <  1:2451441/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CGCTGGTGGCGTGTATTGCCGGTGAAGTTCTTTGGCGGGCTGGGGACACTCGGCGGAGGCATGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATACGTTAATTTCGATTAAAGCGGTATTTATT  >  CP009273/171960‑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: