Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,884,385 G→T 100% V291V (GTC→GTA)  cysJ ← sulfite reductase, alpha subunit, flavoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,884,3850GT97.1% 533.9 / 6.7 172V291V (GTC→GTA) cysJsulfite reductase, alpha subunit, flavoprotein
Reads supporting (aligned to +/- strand):  ref base G (2/3);  new base T (70/97);  total (72/100)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 172 total aligned reads displayed.

TTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTA  >  CP009273/2884240‑2884529
                                                                                                                                                 |                                                                                                                                                
    aCTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGt                                                                                                                                            <  2:2100683/148‑1 (MQ=255)
       gCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCt                                                                                                                                         <  1:2411377/148‑1 (MQ=255)
       gCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCt                                                                                                                                         <  1:2311562/148‑1 (MQ=255)
          ggTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTc                                                                                                                                      <  2:3026605/148‑1 (MQ=255)
            tAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCAc                                                                                                                                    >  1:3422388/1‑148 (MQ=255)
            tAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCAc                                                                                                                                    <  2:112260/148‑1 (MQ=255)
            tAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCAc                                                                                                                                    >  2:2577924/1‑148 (MQ=255)
            tAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCAc                                                                                                                                    >  2:659941/1‑148 (MQ=255)
             aaGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCAcc                                                                                                                                   >  2:781065/1‑148 (MQ=255)
              aGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCa                                                                                                                                  <  2:987250/148‑1 (MQ=255)
               gCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAg                                                                                                                                 >  1:1374795/1‑148 (MQ=255)
                 gTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTg                                                                                                                               <  1:1270072/148‑1 (MQ=255)
                 gTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTg                                                                                                                               <  2:644481/148‑1 (MQ=255)
                  tGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGc                                                                                                                              <  1:769518/148‑1 (MQ=255)
                   ggCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGcc                                                                                                                             >  1:613890/1‑148 (MQ=255)
                     cGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCgg                                                                                                                           <  2:1338421/148‑1 (MQ=255)
                       tAATTCTCAACCATGGTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGAt                                                                                                                         <  2:1959572/148‑1 (MQ=255)
                       tAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGAt                                                                                                                         <  2:3603046/148‑1 (MQ=255)
                         aTTCTCAACCATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCg                                                                                                                       <  2:3274842/148‑1 (MQ=255)
                         aTTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCg                                                                                                                       >  1:856783/1‑148 (MQ=255)
                          ttCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGt                                                                                                                      >  2:1939656/1‑148 (MQ=255)
                           tctcAACACTGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGtt                                                                                                                     <  2:682611/148‑1 (MQ=255)
                              caaccaTGTGGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTg                                                                                                                  <  2:491044/148‑1 (MQ=255)
                              caacaaTGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTg                                                                                                                  >  2:3038453/1‑148 (MQ=255)
                                 caaTGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATa                                                                                                               >  1:3627737/1‑148 (MQ=255)
                                     gTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAg                                                                                                           >  2:264635/1‑148 (MQ=255)
                                     gTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAg                                                                                                           <  1:1435605/148‑1 (MQ=255)
                                          cGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGcc                                                                                                      >  1:2275768/1‑148 (MQ=255)
                                                tGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAgcgc                                                                                                >  2:1197732/1‑148 (MQ=255)
                                                  aCGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGt                                                                                              <  2:910443/148‑1 (MQ=255)
                                                     gTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCAc                                                                                           <  1:1556973/148‑1 (MQ=255)
                                                      tCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCAcc                                                                                          >  2:1221689/1‑148 (MQ=255)
                                                          ttCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGc                                                                                      >  1:1396616/1‑148 (MQ=255)
                                                          ttCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGc                                                                                      >  2:1713490/1‑148 (MQ=255)
                                                          ttCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGc                                                                                      >  1:2982126/1‑148 (MQ=255)
                                                          ttCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGc                                                                                      <  1:395980/148‑1 (MQ=255)
                                                          ttCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGc                                                                                      >  2:533904/1‑148 (MQ=255)
                                                            cGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTg                                                                                    <  1:2513992/148‑1 (MQ=255)
                                                            cGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTg                                                                                    <  1:2513985/148‑1 (MQ=255)
                                                                 tGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAAc                                                                               >  1:2203811/1‑148 (MQ=255)
                                                                 tGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAAc                                                                               >  1:2200199/1‑148 (MQ=255)
                                                                  gCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACg                                                                              >  1:1557705/1‑148 (MQ=255)
                                                                    cACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCa                                                                            >  1:3798982/1‑148 (MQ=255)
                                                                    cACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCa                                                                            >  1:1951507/1‑148 (MQ=255)
                                                                      cTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATg                                                                          <  2:2459401/148‑1 (MQ=255)
                                                                            cgcTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAg                                                                    <  2:3582034/148‑1 (MQ=255)
                                                                              cTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTc                                                                  >  2:359959/1‑148 (MQ=255)
                                                                                tCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCAc                                                                <  2:716741/148‑1 (MQ=255)
                                                                                tCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCAc                                                                <  1:1984429/148‑1 (MQ=255)
                                                                                tCGTTCAGAGGCAACGTTGTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCAc                                                                <  2:2418711/148‑1 (MQ=255)
                                                                                 cGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCAcc                                                               <  2:2779437/148‑1 (MQ=255)
                                                                                   ttCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTa                                                             <  2:3591828/148‑1 (MQ=255)
                                                                                   ttCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTa                                                             <  1:414205/148‑1 (MQ=255)
                                                                                   ttCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTa                                                             <  1:411810/148‑1 (MQ=255)
                                                                                     cAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAg                                                           >  2:3307034/1‑148 (MQ=255)
                                                                                      agagGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGt                                                          >  2:1814111/1‑148 (MQ=255)
                                                                                      agagGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGt                                                          <  1:3378219/148‑1 (MQ=255)
                                                                                        agGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCa                                                        <  1:2789489/148‑1 (MQ=255)
                                                                                          gCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTACCGCATGCCCGAGTCACCTAAGTCAAt                                                      >  2:134921/1‑148 (MQ=255)
                                                                                          gCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAAt                                                      >  2:134915/1‑148 (MQ=255)
                                                                                          gCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAAt                                                      <  2:3422388/148‑1 (MQ=255)
                                                                                            aaCGTTTTGCCCTCGACGGGGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGGTCTTGCCCCAGTGCCGGATCGGTCTGATACCCGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAAttt                                                    <  2:2687392/148‑1 (MQ=255)
                                                                                               gTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGa                                                 <  2:3039250/148‑1 (MQ=255)
                                                                                                ttttGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGat                                                <  1:3508042/148‑1 (MQ=255)
                                                                                                ttttGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGat                                                <  1:3302435/148‑1 (MQ=255)
                                                                                                 tttGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGata                                               <  2:1970767/148‑1 (MQ=255)
                                                                                                 tttGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGata                                               <  1:593901/148‑1 (MQ=255)
                                                                                                  ttGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGatat                                              <  2:279802/148‑1 (MQ=255)
                                                                                                   tGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATg                                             <  1:905529/148‑1 (MQ=255)
                                                                                                       cTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGa                                         <  1:3384602/148‑1 (MQ=255)
                                                                                                       cTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGa                                         <  1:2249486/148‑1 (MQ=255)
                                                                                                         cGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAAc                                       >  2:3159251/1‑148 (MQ=255)
                                                                                                          gACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACg                                      <  2:3627737/148‑1 (MQ=255)
                                                                                                          gACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACg                                      <  2:478754/148‑1 (MQ=255)
                                                                                                           aCGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGt                                     >  1:1802474/1‑148 (MQ=255)
                                                                                                            cGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTc                                    <  2:3283256/148‑1 (MQ=255)
                                                                                                            cGGTGACAGGGTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTc                                    <  2:2277753/148‑1 (MQ=255)
                                                                                                             ggTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCt                                   <  2:1255943/148‑1 (MQ=255)
                                                                                                                    ggTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGAGACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCag                            <  2:1005839/148‑1 (MQ=255)
                                                                                                                      ttCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCagag                          >  1:1868316/1‑148 (MQ=255)
                                                                                                                        cATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGtt                        >  2:1391005/1‑148 (MQ=255)
                                                                                                                         aTCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTa                       <  1:1291819/148‑1 (MQ=255)
                                                                                                                             ccTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGcc                   <  1:2943854/148‑1 (MQ=255)
                                                                                                                              cTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGccc                  >  1:2840491/1‑148 (MQ=255)
                                                                                                                              cTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGccc                  >  1:2497742/1‑148 (MQ=255)
                                                                                                                                ttCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCgg                >  2:3480409/1‑148 (MQ=255)
                                                                                                                                ttCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCgg                <  1:1197733/148‑1 (MQ=255)
                                                                                                                                ttCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCgg                <  1:1197732/148‑1 (MQ=255)
                                                                                                                                  cAGCCACAGCAGTTCGACACTTTCTTTCCCCAGTGGCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCAGGGGGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGACCGTCTTTTTCAGAGTTACGCCCGGTa              <  2:2461580/148‑1 (MQ=255)
                                                                                                                                   aGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTaa             >  2:2360221/1‑148 (MQ=255)
                                                                                                                                   aGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTaa             <  2:2214499/148‑1 (MQ=255)
                                                                                                                                    gCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAAt            <  1:264635/148‑1 (MQ=255)
                                                                                                                                     ccACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAAtt           <  2:1557705/148‑1 (MQ=255)
                                                                                                                                        cagcagTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTc        <  1:3734538/148‑1 (MQ=255)
                                                                                                                                           cagTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTgg     <  1:533904/148‑1 (MQ=255)
                                                                                                                                           cagTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTgg     <  1:1713490/148‑1 (MQ=255)
                                                                                                                                              ttCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTACTTTTCTGGTTa  >  2:717247/1‑148 (MQ=255)
                                                                                                                                              ttCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTa  >  1:2577448/1‑148 (MQ=255)
                                                                                                                                              ttCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTa  >  2:495856/1‑148 (MQ=255)
                                                                                                                                              ttCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTa  >  1:2561115/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTA  >  CP009273/2884240‑2884529

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