Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,659,882 C→A 93.4% S104* (TCG→TAG)  clcB → H(+)/Cl(‑) exchange transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,659,8820CA93.4% 558.4 / 30.6 198S104* (TCG→TAG) clcBH(+)/Cl(‑) exchange transporter
Reads supporting (aligned to +/- strand):  ref base C (9/4);  new base A (101/83);  total (110/88)
Fisher's exact test for biased strand distribution p-value = 3.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

Only 100 of 198 total aligned reads displayed.

ACTCCGGCGCTCGGCGGACTGGCGGCGGGTTTGTTGCTGATGGGCTGGCAGAAATTTACCCAACAACGCCCTCATGCGCCGACCGATTATATGGAAGCGTTGCAAACCGATGGACAGTTCGATTACGCAGCAAGCCTGGTTAAATCGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCGATGATTCTTTTAGCTGCCCTTGCCGCCTCCTGTTTTGCCCAACGTTTTACGCCACGCCAGGAGTGGAAATTATGGATCGCCTGTGGG  >  CP009273/1659737‑1660024
                                                                                                                                                 |                                                                                                                                              
 cTCCGGCGCTCGGCGGACTGGCGGCGGGTTTGTTGCTGATGGGCTGGCAGAAATTTACCCAACAACGCCCTCATGCGCCGACCGATTATATGGAAGCGTTGCAAACCGATGGACAGTTCGATTACGCAGCAAGCCTGGTTAAATAGCt                                                                                                                                             >  1:1458677/1‑148 (MQ=255)
     ggCGCTCGGCGGACTGGCGGCGGGTTTGTTGCTGATGGGCTGGCAGAAATTTACCCAACAACGCCCTCATGCGCCGACCGATTATATGGAAGCGTTGCAAACCGATGGACAGTTCGATTACGCAGCAAGCCTGGTTAAATAGCTTGcc                                                                                                                                         <  1:585670/148‑1 (MQ=255)
      gcgcTCGGCGGACTGGCGGCGGGTTTGTTGCTGATGGGCTGGCAGAAATTTACCCAACAACGCCCTCATGCGCCGACCGATTATATGGAAGCGTTGCAAACCGATGGACAGTTCGATTACGCAGCAAGCCTGGTTAAATAGCTTGCct                                                                                                                                        <  2:940378/148‑1 (MQ=255)
      gcgcTCGGCGGACTGGCGGCGGGTTTGTTGCTGATGGGCTGGCAGAAATTTACCCAACAACGCCCTCATGCGCCGACCGATTATATGGAAGCGTTGCAAACCGATGGACAGTTCGATTACGCAGCAAGCCTGGTTAAATAGCTTGCct                                                                                                                                        >  2:2831466/1‑148 (MQ=255)
        gcTCGGCGGACTGGCGGCGGGTTTGTTGCTGATGGGCTGGCAGAAATTTACCCAACAACGCCCTCATGCGCCGACCGATTATATGGAAGCGTTGCAAACCGATGGACAGTTCGATTACGCAGCAAGCCTGGTTAAATAGCTTGCctct                                                                                                                                      <  2:840434/148‑1 (MQ=255)
         cTCGGCGGACTGGCGGCGGGTTTGTTGCTGATGGGCTGGCAGAAATTTACCCAACAACGCCCTCATGCGCCGACCGATTATATGGAAGCGTTGCAAACCGATGGACAGTTCGATTACGCAGCAAGCCTGGTTAAATCGCTTGCctctc                                                                                                                                     <  2:1741880/148‑1 (MQ=255)
         cTCGGCGGACTGGCGGCGGGTTTGTTGCTGATGGGCTGGCAGAAATTTACCCAACAACGCCCTCATGCGCCGACCGATTATATGGAAGCGTTGCAAACCGATGGACAGTTCGATTACGCAGCAAGCCTGGTTAAATAGCTTGCctctc                                                                                                                                     <  1:1821994/148‑1 (MQ=255)
           cggcggACTGGCGGCGGGTTTGTTGCTGATGGGCTGGCAGAAATTTACCCAACAACGCCCTCATGCGCCGACCGATTATATGGAAGCGTTGCAAACCGATGGACAGTTCGATTACGCAGCAAGCCTGGTTAAATAGCTTGCCTCTctg                                                                                                                                   <  1:2210748/148‑1 (MQ=255)
            ggcggACTGGCGGCGGGTTTGTTGCTGATGGGCTGGCAGAAATTTACCCAACAACGCCCTCATGCGCCGACCGATTATATGGAAGCGTTGCAAACCGATGGACAGTTCGATTACGCAGCAAGCCTGGTTAAATAGCTTGCCTCTctgc                                                                                                                                  <  2:1836995/148‑1 (MQ=255)
               ggACTGGCGGCGGGTTTGTTGCTGATGGGCTGGCAGAAATTTACCCAACAACGCCCTCATGCGCCGACCGATTATATGGAAGCGTTGCAAACCGATGGACAGTTCGATTACGCAGCAAGCCTGGTTAAATAGCTTGCCTCTCTGCTgg                                                                                                                               <  1:1254410/148‑1 (MQ=255)
                      cggcggGTTTGTTGCTGATGGGCTGGCAGAAATTTACCCAACAACGCCCTCATGCGCCGACCGATTATATGGAAGCGTTGCAAACCGATGGACAGTTCGATTACGCAGCAAGCCTGGTTAAATCGCTTGCCTCTCTGCTGGTAGTAAc                                                                                                                        <  1:895743/148‑1 (MQ=255)
                            gTTTGTTGCTGATGGGCTGGCAGAAATTTACCCAACAACGCCCTCATGCGCCGACCGATTATATGGAAGCGTTGCAAACCGATGGACAGTTCGATTACGCAGCAAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCgg                                                                                                                  >  2:1595155/1‑148 (MQ=255)
                             tttGTTGCTGATGGGCTGGCAGAAATTTACCCAACAACGCCCTCATGCGCCGACCGATTATATGGAAGCGTTGCAAACCGATGGACAGTTCGATTACGCAGCAAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGc                                                                                                                 >  1:2848775/1‑148 (MQ=255)
                             tttGTTGCTGATGGGCTGGCAGAAATTTACCCAACAACGCCCTCATGCGCCGACCGATTATATGGAAGCGTTGCAAACCGATGGACAGTTCGATTACGCAGCAAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGc                                                                                                                 >  1:1097088/1‑148 (MQ=255)
                             tttGTTGCTGATGGGCTGGCAGAAATTTACCCAACAACGCCCTCATGCGCCGACCGATTATATGGAAGCGTTGCAAACCGATGGACAGTTCGATTACGCAGCAAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGc                                                                                                                 >  1:198132/1‑148 (MQ=255)
                                gttgCTGATGGGCTGGCAGAAATTTACCCAACAACGCCCTCATGCGCCGACCGATTATATGGAAGCGTTGCAAACCGATGGACAGTTCGATTACGCAGCAAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGt                                                                                                              >  2:1709338/1‑148 (MQ=255)
                                 ttgCTGATGGGCTGGCAGAAATTTACCCAACAACGCCCTCATGCGCCGACCGATTATATGGAAGCGTTGCAAACCGATGGACAGTTCGATTACGCAGCAAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTg                                                                                                             >  1:1786344/1‑148 (MQ=255)
                                      gATGGGCTGGCAGAAATTTACCCAACAACGCCCTCATGCGCCGACCGATTATATGGAAGCGTTGCAAACCGATGGACAGTTCGATTACGCAGCAAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAAtt                                                                                                        <  2:219061/148‑1 (MQ=255)
                                        tGGGCTGGCAGAAATTTACCCAACAACGCCCTCATGCGCCGACCGATTATATGGAAGCGTTGCAAACCGATGGACAGTTCGATTACGCAGCAAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTgg                                                                                                      <  1:3108640/148‑1 (MQ=255)
                                          ggCTGGCAGAAATTTACCCAACAACGCCCTCATGCGCCGACCGATTATATGGAAGCGTTGCAAACCGATGGACAGTTCGATTACGCAGCAAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTc                                                                                                    >  1:1575912/1‑148 (MQ=255)
                                              ggCAGAAATTTACCCAACAACGCCCTCATGCGCCGACCGATTATATGGAAGCGTTGCAAACCGATGGACAGTTCGATTACGCAGCAAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGa                                                                                                >  1:2808500/1‑148 (MQ=255)
                                                cAGAAATTTACCCAACAACGCCCTCATGCGCCGACCGATTATATGGAAGCGTTGCAAACCGATGGACAGTTCGATTACGCAGCAAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAg                                                                                              >  2:240390/1‑148 (MQ=255)
                                                   aaaTTTACCCAACAACGCCCTCATGCGCCGACCGATTATATGGAAGCGTTGCAAACCGATGGACAGTTCGATTACGCAGCAAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTg                                                                                           <  1:287391/148‑1 (MQ=255)
                                                      tttACCCAACAACGCCCTCATGCGCCGACCGATTATATGGAAGCGTTGCAAACCGATGGACAGTTCGATTACGCAGCAAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCga                                                                                        >  1:3487109/1‑148 (MQ=255)
                                                         aCCCAACAACGCCCTCATGCGCCGACCGATTATATGGAAGCGTTGCAAACCGATGGACAGTTCGATTACGCAGCAAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCgatga                                                                                     >  1:3009846/1‑148 (MQ=255)
                                                            caacaaCGCCCTCATGCGCCGACCGATTATATGGAAGCGTTGCAAACCGATGGACAGTTCGATTACGCAGCAAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCGATGATTc                                                                                  >  1:684691/1‑148 (MQ=255)
                                                              acaacGCCCTCATGCGCCGACCGATTATATGGAAGCGTTGCAAACCGATGGACAGTTCGATTACGCAGCAAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCGATGATTCtt                                                                                <  2:1786344/148‑1 (MQ=255)
                                                              acaacGCCCTCATGCGCCGACCGATTATATGGAAGCGTTGCAAACCGATGGACAGTTCGATTACGCAGCAAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCGATGATTCtt                                                                                <  2:422770/148‑1 (MQ=255)
                                                                 acGCCCTCATGCGCCGACCGATTATATGGAAGCGTTGCAAACCGATGGACAGTTCGATTACGCAGCAAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCGATGATTCTTTTa                                                                             >  1:1011254/1‑148 (MQ=255)
                                                                 acGCCCTCATGCGCCGACCGATTATATGGAAGCGTTGCAAACCGATGGACAGTTCGATTACGCAGCAAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCGATGATTCTTTTa                                                                             <  2:3234283/148‑1 (MQ=255)
                                                                  cGCCCTCATGCGCCGACCGATTATATGGAAGCGTTGCAAACCGATGGACAGTTCGATTACGCAGCAAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCGATGATTCTTTTAg                                                                            >  1:225686/1‑148 (MQ=255)
                                                                  cGCCCTCATGCGCCGACCGATTATATGGAAGCGTTGCAAACCGATGGACAGTTCGATTACGCAGCAAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCGATGATTCTTTTAg                                                                            <  1:2545644/148‑1 (MQ=255)
                                                                   gCCCTCATGCGCCGACCGATTATATGGAAGCGTTGCAAACCGATGGACAGTTCGATTACGCAGCAAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCGATGATTCTTTTAGc                                                                           <  1:3439155/148‑1 (MQ=255)
                                                                    cccTCATGCGCCGACCGATTATATGGAAGCGTTGCAAACCGATGGACAGTTCGATTACGCAGCAAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCGATGATTCTTTTAGCt                                                                          >  2:2533383/1‑148 (MQ=255)
                                                                     ccTCATGCGCCGACCGATTATATGGAAGCGTTGCAAACCGATGGACAGTTCGATTACGCAGCAAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAATGTGCGATGATTCTTTTAGCTg                                                                         >  1:177662/1‑148 (MQ=255)
                                                                     ccTCATGCGCCGACCGATTATATGGAAGCGTTGCAAACCGATGGACAGTTCGATTACGCAGCAAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCGATGATTCTTTTAGCtt                                                                         >  2:2506046/1‑147 (MQ=255)
                                                                     ccTCATGCGCCGACCGATTATATGGAAGCGTTGCAAACCGATGGACAGTTCGATTACGCAGCAAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCGATGATTCTTTTAGCTg                                                                         >  1:3301992/1‑148 (MQ=255)
                                                                      cTCATGCGCCGACCGATTATATGGAAGCGTTGCAAACCGATGGACAGTTCGATTACGCAGCAAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCGATGATTCTTTTAGCTGc                                                                        <  2:1069367/148‑1 (MQ=255)
                                                                      cTCATGCGCCGACCGATTATATGGAAGCGTTGCAAACCGATGGACAGTTCGATTACGCAGCAAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCGATGATTCTTTTAGCTGc                                                                        >  2:3468458/1‑148 (MQ=255)
                                                                        cATGCGCCGACCGATTATATGGAAGCGTTGCAAACCGATGGACAGTTCGATTACGCAGCAAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCGATGATTCTTTTAGCTGccc                                                                      <  1:3075926/148‑1 (MQ=255)
                                                                          tGCGCCGACCGATTATATGGAAGCGTTGCAAACCGATGGACAGTTCGATTACGCAGCAAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCGATGATTCTTTTAGCTGCCCtt                                                                    >  2:1970711/1‑148 (MQ=255)
                                                                           gcgcCGACCGATTATATGGAAGCGTTGCAAACCGATGGACAGTTCGATTACGCAGCAAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCGATGATTCTTTTAGCTGCCCTTg                                                                   <  2:751778/148‑1 (MQ=255)
                                                                            cgcCGACCGATTATATGGAAGCGTTGCAAACCGATGGACAGTTCGATTACGCAGCAAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCGATGATTCTTTTAGCTGCCCTTgc                                                                  >  2:275352/1‑148 (MQ=255)
                                                                            cgcCGACCGATTATATGGAAGCGTTGCAAACCGATGGACAGTTCGATTACGCAGCAAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCGATGATTCTTTTAGCTGCCCTTgc                                                                  >  2:275924/1‑148 (MQ=255)
                                                                                 accgaTTATATGGAAGCGTTGCAAACCGATGGACAGTTCGATTACGCAGCAAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCGATGATTCTTTTAGCTGCCCTTGCCGcct                                                             >  2:2523921/1‑148 (MQ=255)
                                                                                  ccgaTTATATGGAAGCGTTGCAAACCGATGGACAGTTCGATTACGCAGCAAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCGATGATTCTTTTAGCTGCCCTTGCCGcctc                                                            >  2:1869985/1‑148 (MQ=255)
                                                                                   cgaTTATCTGGAAGCGTTGCAAACGGATGGACAGTTCGATTACGCAGCAAGCCTGGTTATATTGCTTGCCTCTCTGCTGGTAGTCACCAGCGGCAGTGCTATTGGTCGCGTAGGTGCGATGATTCTTTTAGCTGCCCTTGCCGcctcc                                                           <  1:2927672/148‑1 (MQ=255)
                                                                                      ttATATGGAAGCGTTGCAAACCGATGGACAGTTCGATTACGCAGCAAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCGATGATTCTTTTAGCTGCCCTTGCCGCCTCCTGt                                                        >  1:2738152/1‑148 (MQ=255)
                                                                                      ttATATGGAAGCGTTGCAAACCGATGGACAGTTCGATTACGCAGCAAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCGATGATTCTTTTAGCTGCCCTTGCCGCCTCCTGt                                                        >  1:2737898/1‑148 (MQ=255)
                                                                                        atatGGAAGCGTTGCAAACCGATGGACAGTTCGATTACGCAGCAAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCGATGATTCTTTTAGCTGCCCTTGCCGCCTCCTGttt                                                      <  1:1416974/148‑1 (MQ=255)
                                                                                          atGGAAGCGTTGCAAACCGATGGACAGTTCGATTACGCAGCAAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCGATGATTCTTTTAGCTGCCCTTGCCGCCTCCTGTTTTg                                                    >  2:3438873/1‑148 (MQ=255)
                                                                                           tGGAAGCGTTGCAAACCGATGGACAGTTCGATTACGCAGCAAGCCTGGTTAAATCGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCGATGATTCTTTTAGCTGCCCTTGCCGCCTCCTGTTTTGc                                                   >  1:1514562/1‑148 (MQ=255)
                                                                                            ggAAGCGTTGCAAACCGATGGACAGTTCGATTACGCAGCAAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCGATGATTCTTTTAGCTGCCCTTGCCGCCTCCTGTTTTGcc                                                  >  1:1593160/1‑148 (MQ=255)
                                                                                              aaGCGTTGCAAACCGATGGACCGTTCGATTACGCAGCAAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCACTTGGTCGCGAAGGTGCGATGATTCGTTTAGCTGCCCTTGCCGCGTCGTGTTTTGgcct                                                >  2:2105459/1‑147 (MQ=255)
                                                                                                 cGTTGCAAACCGATGGACAGTTCGATTACGCAGCAAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCGATGATTCTTTTAGCTGCCCTTGCCGCCTCCTGTTTTGCCCAACg                                             >  1:2705310/1‑148 (MQ=255)
                                                                                                   ttGCAAACCGATGGACAGTTCGATTACGCAGCAAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCGATGATTCTTTTAGCTGCCCTTGCCGCCTCCTGTTTTGCCCAACGtt                                           <  2:1972870/148‑1 (MQ=255)
                                                                                                   ttGCAAACCGATGGACAGTTCGATTACGCAGCAAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCGATGATTCTTTTAGCTGCCCTTGCCGCCTCCTGTTTTGCCCAACGtt                                           <  2:1458677/148‑1 (MQ=255)
                                                                                                    tGCAAACCGATGGACAGTTCGATTACGCAGCAAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCGATGATTCTTTTAGCTGCCCTTGCCGCCTCCTGTTTTGCCCAACGttt                                          >  2:2568855/1‑148 (MQ=255)
                                                                                                    tGCAAACCGATGGACAGTTCGATTACGCAGCAAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCGATGATTCTTTTAGCTGCCCTTGCCGCCTCCTGTTTTGCCCAACGttt                                          <  1:2352457/148‑1 (MQ=255)
                                                                                                     gCAAACCGATGGACAGTTCGATTACGCAGCAAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCGATGATTCTTTTAGCTGCCCTTGCCGCCTCCTGTTTTGCCCAACGtttt                                         <  1:1670003/148‑1 (MQ=255)
                                                                                                      cAAACCGATGGACAGTTCGATTACGCAGCAAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCGATGATTCTTTTAGCTGCCCTTGCCGCCTCCTGTTTTGCCCAACGTTTTa                                        >  1:1708019/1‑148 (MQ=255)
                                                                                                      cAAACCGATGGACAGTTCGATTACGCAGCAAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCGATGATTCTTTTAGCTGCCCTTGCCGCCTCCTGTTTTGCCCAACGTTTTa                                        <  2:2237407/148‑1 (MQ=255)
                                                                                                        aaCCGATGGACAGTTCGATTACGCAGCAAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCGATGATTCTTTTAGCTGCCCTTGCCGCCTCCTGTTTTGCCCAACGTTTTacg                                      >  2:2090163/1‑148 (MQ=255)
                                                                                                         aCCGATGGACAGTTCGATTACGCAGCAAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCGATGATTCTTTTAGCTGCCCTTGCCGCCTCCTGTTTTGCCCAACGTTTTacgc                                     >  2:1960373/1‑148 (MQ=255)
                                                                                                           cGATGGACAGTTCGATTACGCAGCAAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCGATGATTCTTTTAGCTGCCCTTGCCGCCTCCTGTTTTGCCCAACGTTTTacgcca                                   <  2:621905/148‑1 (MQ=255)
                                                                                                             aTGGACAGTTCGATTACGCAGCAAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCGATGATTCTTTTAGCTGCCCTTGCCGCCTCCTGTTTTGCCCAACGTTTTacgccacg                                 <  1:101391/148‑1 (MQ=255)
                                                                                                              tGGACAGTTCGATTACGCAGCAAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCGATGATTCTTTTAGCTGCCCTTGCCGCCTCCTGTTTTGCCCAACGTTTTacgccacgc                                >  1:2735222/1‑148 (MQ=255)
                                                                                                               ggACAGTTCGATTACGCAGCAAGCCTGGTTAAATCGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCGATGATTCTTTTAGCTGCCCTTGCCGCCTCCTGTTTTGCCCAACGTTTTacgccacgcc                               <  1:3043686/148‑1 (MQ=255)
                                                                                                                 aCAGTTCGATTACGCAGCAAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCGATGATTCTTTTAGCTGCCCTTGCCGCCTCCTGTTTTGCCCAACGTTTTACGCCACGCCAg                             >  1:695810/1‑148 (MQ=255)
                                                                                                                  cAGTTCGATTACGCAGCAAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCGATGATTCTTTTAGCTGCCCTTGCCGCCTCCTGTTTTGCCCAACGTTTTACGCCACGCCAgg                            >  2:3345943/1‑148 (MQ=255)
                                                                                                                    gTTCGATTACGCAGCAAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCGATGATTCTTTTAGCTGCCCTTGCCGCCTCCTGTTTTGCCCAACGTTTTACGCCACGCCAGGAg                          <  2:1097121/148‑1 (MQ=255)
                                                                                                                    gTTCGATTACGCAGCAAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCGATGATTCTTTTAGCTGCCCTTGCCGCCTCCTGTTTTGCCCAACGTTTTACGCCACGCCAGGAg                          <  2:2848802/148‑1 (MQ=255)
                                                                                                                     ttCGATTACGCAGCAAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCGATGATTCTTTTAGCTGCCCTTGCCGCCTCCTGTTTTGCCCAACGTTTTACGCCACGCCAGGAGt                         >  2:1265167/1‑148 (MQ=255)
                                                                                                                     ttCGATTACGCAGCAAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCGATGATTCTTTTAGCTGCCCTTGCCGCCTCCTGTTTTGCCCAACGTTTTACGCCACGCCAGGAGt                         >  1:3453119/1‑148 (MQ=255)
                                                                                                                     ttCGATTACGCAGCAAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCGATGATTCTTTTAGCTGCCCTTGCCGCCTCCTGTTTTGCCCAACGTTTTACGCCACGCCAGGAGt                         <  1:2021319/148‑1 (MQ=255)
                                                                                                                        gATGACGCAGCAAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCGATGATTCTTTTAGCTGCCCGTGCCGCCTCCTGTTTTGCCCAACGTTTTACGCCACGCCAGGAGTGGa                      <  2:1803625/148‑1 (MQ=255)
                                                                                                                         aTTACGCAGCAAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCGATGATTCTTTTAGCTGCCCTTGCCGCCTCCTGTTTTGCCCAACGTTTTACGCCACGCCAGGAGTGGaa                     <  1:397362/148‑1 (MQ=255)
                                                                                                                          ttACGCAGCAAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCGATGATTCTTTTAGCTGCCCTTGCCGCCTCCTGTTTTGCCCAACGTTTTACGCCACGCCAGGAGTGGaaa                    >  1:631668/1‑148 (MQ=255)
                                                                                                                          ttACGCAGCAAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCGATGATTCTTTTAGCTGCCCTTGCCGCCTCCTGTTTTGCCCAACGTTTTACGCCACGCCAGGAGTGGaaa                    >  1:307123/1‑148 (MQ=255)
                                                                                                                            aCGCAGCAAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCGATGATTCTTTTAGCTGCCCTTGCCGCCTCCTGTTTTGCCCAACGTTTTACGCCACGCCAGGAGTGGAAAtt                  <  2:160001/148‑1 (MQ=255)
                                                                                                                            aCGCAGCAAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCGATGATTCTTTTAGCTGCCCTTGCCGCCTCCTGTTTTGCCCAACGTTTTACGCCACGCCAGGAGTGGAAAtt                  <  2:814577/148‑1 (MQ=255)
                                                                                                                             cGCAGCAAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCGATGATTCTTTTAGCTGCCCTTGCCGCCTCCTGTTTTGCCCAACGTTTTACGCCACGCCAGGAGTGGAAATTa                 <  1:999079/148‑1 (MQ=255)
                                                                                                                             cGCAGCAAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCGATGATTCTTTTAGCTGCCCTTGCCGCCTCCTGTTTTGCCCAACGTTTTACGCCACGCCAGGAGTGGAAATTa                 >  1:1123357/1‑148 (MQ=255)
                                                                                                                               cagcaAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCGATGATTCTTTTAGCTGCCCTTGCCGCCTCCTGTTTTGCCCAACGTTTTACGCCACGCCAGGAGTGGAAATTATg               >  2:1515942/1‑148 (MQ=255)
                                                                                                                               cagcaAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCGATGATTCTTTTAGCTGCCCTTGCCGCCTCCTGTTTTGCCCAACGTTTTACGCCACGCCAGGAGTGGAAATTATg               >  2:2631859/1‑148 (MQ=255)
                                                                                                                                agcaAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCGATGATTCTTTTAGCTGCCCTTGCCGCCTCCTGTTTTGCCCAACGTTTTACGCCACGCCAGGAGTGGAAATTATgg              <  1:3362630/148‑1 (MQ=255)
                                                                                                                                 gcaAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCGATGATTCTTTTAGCTGCCCTTGCCGCCTCCTGTTTTGCCCAACGTTTTACGCCACGCCAGGAGTGGAAATTATGGa             <  1:1143185/148‑1 (MQ=255)
                                                                                                                                 gcaAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCGATGATTCTTTTAGCTGCCCTTGCCGCCTCCTGTTTTGCCCAACGTTTTACGCCACGCCAGGAGTGGAAATTATGGa             <  1:1143199/148‑1 (MQ=255)
                                                                                                                                 gcaAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCGATGATTCTTTTAGCTGCCCTTGCCGCCTCCTGTTTTGCCCAACGTTTTACGCCACGCCAGGAGTGGAAATTATGGa             <  1:2228575/148‑1 (MQ=255)
                                                                                                                                 gcaAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCGATGATTCTGTTAGCTGCCCTTGCCGCCTCCTGTTTTGCCCAACGGTTTACGCCACGCCAGGAGTGGAAATTCTGGa             >  2:992568/1‑148 (MQ=255)
                                                                                                                                  caAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCGATGATTCTTTTAGCTGCCCTTGCCGCCTCCTGTTTTGCCCAACGTTTTACGCCACGCCAGGAGTGGAAATTATGGAt            >  1:2181666/1‑148 (MQ=255)
                                                                                                                                  caAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCGATGATTCTTTTAGCTGCCCTTGCCGCCTCCTGTTTTGCCCAACGTTTTACGCCACGCCAGGAGTGGAAATTATGGAt            <  2:3021510/148‑1 (MQ=255)
                                                                                                                                  caAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCGATGATTCTTTTAGCTGCCCTTGCCGCCTCCTGTTTTGCCCAACGTTTTACGCCACGCCAGGAGTGGAAATTATGGAt            >  1:830305/1‑148 (MQ=255)
                                                                                                                                  caAGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCGATGATTCTTTTAGCTGCCCTTGCCGCCTCCTGTTTTGCCCAACGTTTTACGCCACGCCAGGAGTGGAAATTATGGAt            >  2:1876037/1‑148 (MQ=255)
                                                                                                                                    aGCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCGATGATTCTTTTAGCTGCCCTTGCCGCCTCCTGTTTTGCCCAACGTTTTACGCCACGCCAGGAGTGGACATTATGGATCg          >  2:2612704/1‑148 (MQ=255)
                                                                                                                                     gCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCGATGATTCTTTTAGCTGCCCTTGCCGCCTCCTGTTTTGCCCAACGTTTTACGCCACGCCAGGAGTGGAAATTATGGATCGc         >  1:2457690/1‑148 (MQ=255)
                                                                                                                                     gCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCGATGATTCTTTTAGCTGCCCTTGCCGCCTCCTGTTTTGCCCAACGTTTTACGCCACGCCAGGAGTGGAAATTATGGATCGc         >  1:1419753/1‑148 (MQ=255)
                                                                                                                                     gCCTGGTTAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCGATGATTCTTTTAGCTGCCCTTGCCGCCTCCTGTTTTGCCCAACGTTTTACGCCACGCCAGGAGTGGAAATTATGGATCGc         <  1:3249517/148‑1 (MQ=255)
                                                                                                                                            tAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCTATGATTCTTTTAGCTGCCCTTGCCGCCTCCTGTTTTGCCCAACGTTTTACGCCACGCCAGGAGTGGAAATTATGGATCGCCTGTggg  >  1:232133/1‑148 (MQ=255)
                                                                                                                                            tAAATAGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCTATGATTCTTTTAGCTGCCCTTGCCGCCTCCTGTTTTGCCCAACGTTTTACGCCACGCCAGGAGTGGAAATTATGGATCGCCTGTggg  >  1:232063/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
ACTCCGGCGCTCGGCGGACTGGCGGCGGGTTTGTTGCTGATGGGCTGGCAGAAATTTACCCAACAACGCCCTCATGCGCCGACCGATTATATGGAAGCGTTGCAAACCGATGGACAGTTCGATTACGCAGCAAGCCTGGTTAAATCGCTTGCCTCTCTGCTGGTAGTAACCAGCGGCAGTGCAATTGGTCGCGAAGGTGCGATGATTCTTTTAGCTGCCCTTGCCGCCTCCTGTTTTGCCCAACGTTTTACGCCACGCCAGGAGTGGAAATTATGGATCGCCTGTGGG  >  CP009273/1659737‑1660024

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: