Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,884,385 G→T 88.4% 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,3850GT88.4% 383.1 / 47.0 155V291V (GTC→GTAcysJsulfite reductase, alpha subunit, flavoprotein
Reads supporting (aligned to +/- strand):  ref base G (10/8);  new base T (67/70);  total (77/78)
Fisher's exact test for biased strand distribution p-value = 6.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 155 total aligned reads displayed.

GCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTA  >  CP009273/2884247‑2884529
                                                                                                                                          |                                                                                                                                                
gCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCt                                                                                                                                         <  1:2741924/148‑1 (MQ=255)
gCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCt                                                                                                                                         <  2:292147/148‑1 (MQ=255)
  gggTAAGCGTGGCGTAATTCTCAACCATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCttt                                                                                                                                       <  2:1715555/148‑1 (MQ=255)
  gggTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCttt                                                                                                                                       <  2:2557729/148‑1 (MQ=255)
  gggTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCttt                                                                                                                                       <  1:697606/148‑1 (MQ=255)
  gggTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCttt                                                                                                                                       <  1:2359472/148‑1 (MQ=255)
         cGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGt                                                                                                                                >  2:1883199/1‑148 (MQ=255)
           tGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGc                                                                                                                              >  2:176841/1‑148 (MQ=255)
           tGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGc                                                                                                                              <  2:208533/148‑1 (MQ=255)
           tGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGc                                                                                                                              <  2:3106589/148‑1 (MQ=255)
           tGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGc                                                                                                                              >  1:1597808/1‑148 (MQ=255)
           tGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGc                                                                                                                              >  1:2160372/1‑148 (MQ=255)
             gCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCg                                                                                                                            >  1:1714140/1‑148 (MQ=255)
                   ttCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGt                                                                                                                      >  2:3268542/1‑148 (MQ=255)
                    tctcAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGtt                                                                                                                     <  1:177469/148‑1 (MQ=255)
                    tctcAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGtt                                                                                                                     <  1:1845633/148‑1 (MQ=255)
                    tctcAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGtt                                                                                                                     <  1:1054840/148‑1 (MQ=255)
                     ctcAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTc                                                                                                                    >  2:2087446/1‑148 (MQ=255)
                     ctcAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTc                                                                                                                    <  1:2052076/148‑1 (MQ=255)
                         acaaTGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGAt                                                                                                                >  2:1164176/1‑148 (MQ=255)
                             tGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCa                                                                                                            >  1:1195560/1‑148 (MQ=255)
                              gTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAg                                                                                                           >  2:1888897/1‑148 (MQ=255)
                                tGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGAc                                                                                                         >  1:965573/1‑148 (MQ=255)
                                  gCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGc                                                                                                       >  2:139777/1‑148 (MQ=255)
                                  gCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGc                                                                                                       >  1:308628/1‑148 (MQ=255)
                                   cGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGcc                                                                                                      >  2:1857277/1‑148 (MQ=255)
                                      tgtTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAg                                                                                                   >  1:384876/1‑148 (MQ=255)
                                      tgtTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAg                                                                                                   >  1:714666/1‑148 (MQ=255)
                                           aCGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGt                                                                                              >  2:2454612/1‑148 (MQ=255)
                                             ggTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCa                                                                                            <  1:2271639/148‑1 (MQ=255)
                                               tCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCAcc                                                                                          >  1:1579076/1‑148 (MQ=255)
                                                cAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCAccc                                                                                         >  1:417763/1‑148 (MQ=255)
                                                  gTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCgg                                                                                       <  1:2295708/148‑1 (MQ=255)
                                                  gTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCgg                                                                                       <  2:227907/148‑1 (MQ=255)
                                                  gTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCgg                                                                                       >  1:1587474/1‑148 (MQ=255)
                                                      gAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTgg                                                                                   >  2:41682/1‑148 (MQ=255)
                                                       aaGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGt                                                                                  >  1:9517/1‑148 (MQ=255)
                                                       aaGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGt                                                                                  >  1:8412/1‑148 (MQ=255)
                                                        aGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTa                                                                                 >  1:1051594/1‑148 (MQ=255)
                                                        aGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTa                                                                                 >  1:1051113/1‑148 (MQ=255)
                                                         gTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTaa                                                                                <  1:1703270/148‑1 (MQ=255)
                                                         gTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGATGGTaa                                                                                >  2:2746063/1‑148 (MQ=255)
                                                          tGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAAc                                                                               >  1:79145/1‑148 (MQ=255)
                                                             cACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCa                                                                            <  2:916511/148‑1 (MQ=255)
                                                             cACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCa                                                                            <  2:1281612/148‑1 (MQ=255)
                                                               cTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATg                                                                          <  1:315888/148‑1 (MQ=255)
                                                               cTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATg                                                                          <  1:2168580/148‑1 (MQ=255)
                                                                     cgcTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAg                                                                    >  2:2380984/1‑148 (MQ=255)
                                                                      gcTTCGTTCAGAGGCAACGTTTTGCCCTCGACTGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACTAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGt                                                                   <  1:3081995/148‑1 (MQ=255)
                                                                      gcTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGt                                                                   <  2:1250554/148‑1 (MQ=255)
                                                                      gcTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGt                                                                   >  1:2386823/1‑148 (MQ=255)
                                                                        ttCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCa                                                                 >  1:533307/1‑148 (MQ=255)
                                                                         tCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCAGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCAc                                                                <  2:2648353/148‑1 (MQ=255)
                                                                         tCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCAc                                                                >  1:373084/1‑148 (MQ=255)
                                                                          cGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCAcc                                                               >  2:2989061/1‑148 (MQ=255)
                                                                            ttCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTa                                                             <  1:3113306/148‑1 (MQ=255)
                                                                             tCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTaa                                                            <  1:2293761/148‑1 (MQ=255)
                                                                               agagGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGt                                                          <  2:3170555/148‑1 (MQ=255)
                                                                                gagGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTc                                                         <  2:278582/148‑1 (MQ=255)
                                                                                    cAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAAtt                                                     <  1:2116590/148‑1 (MQ=255)
                                                                                    cAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAAtt                                                     <  1:2116473/148‑1 (MQ=255)
                                                                                    cAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAAtg                                                     >  2:1562296/1‑147 (MQ=255)
                                                                                     aaCGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAAttt                                                    >  1:164780/1‑148 (MQ=255)
                                                                                       cGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCg                                                  <  2:2160372/148‑1 (MQ=255)
                                                                                        gTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGa                                                 <  2:1516820/148‑1 (MQ=255)
                                                                                          tttGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGata                                               >  1:3117218/1‑148 (MQ=255)
                                                                                           tttcccTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGatat                                              <  2:1201792/145‑1 (MQ=255)
                                                                                             gCCCTCGACGGTGACAGGTTCCTCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTGTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATgg                                            <  2:2008845/148‑1 (MQ=255)
                                                                                                 tCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGaa                                        <  1:1883199/148‑1 (MQ=255)
                                                                                                    aCGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGt                                     >  2:450901/1‑148 (MQ=255)
                                                                                                    aCGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGt                                     >  2:450906/1‑148 (MQ=255)
                                                                                                     cGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTc                                    >  1:106954/1‑148 (MQ=255)
                                                                                                     cGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTc                                    >  1:1780695/1‑148 (MQ=255)
                                                                                                      ggTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCt                                   <  1:1857277/148‑1 (MQ=255)
                                                                                                         gACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCtttt                                <  2:1191263/148‑1 (MQ=255)
                                                                                                          aCAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCttttt                               <  1:3273780/148‑1 (MQ=255)
                                                                                                          aCAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCttttt                               >  1:318442/1‑148 (MQ=255)
                                                                                                          aCAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCttttt                               <  1:1248061/148‑1 (MQ=255)
                                                                                                          aCAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCttttt                               <  1:2087449/148‑1 (MQ=255)
                                                                                                          aCAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCttttt                               <  1:2087446/148‑1 (MQ=255)
                                                                                                              ggtCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGGGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCaga                           <  2:811443/146‑1 (MQ=255)
                                                                                                              gTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCaga                           <  1:41682/148‑1 (MQ=255)
                                                                                                               ttCATCGCCTTTCAGCCACATCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCagag                          <  2:2253510/148‑1 (MQ=255)
                                                                                                                tCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGt                         >  1:447726/1‑148 (MQ=255)
                                                                                                                tCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGt                         <  2:373084/148‑1 (MQ=255)
                                                                                                                 cATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGtt                        <  1:2856277/148‑1 (MQ=255)
                                                                                                                    cGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACg                     >  1:716659/1‑148 (MQ=255)
                                                                                                                     gCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGc                    >  1:3217805/1‑148 (MQ=255)
                                                                                                                     gCCTTTCAGCCACAGCAGGTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGc                    <  2:2549932/148‑1 (MQ=255)
                                                                                                                          tCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGt               <  2:1216318/148‑1 (MQ=255)
                                                                                                                          tCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGt               <  1:280084/148‑1 (MQ=255)
                                                                                                                            aGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTaa             <  2:2066006/148‑1 (MQ=255)
                                                                                                                            aGCCACAGCAGTTCGACAAGTTCTGTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTaa             <  2:2176012/148‑1 (MQ=255)
                                                                                                                                acaGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAAtttt         >  2:856277/1‑148 (MQ=255)
                                                                                                                                 cagcagTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTc        >  1:2743363/1‑148 (MQ=255)
                                                                                                                                  agcagTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCt       >  1:1687731/1‑148 (MQ=255)
                                                                                                                                    cagTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTgg     >  2:1590478/1‑148 (MQ=255)
                                                                                                                                    cagTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTgg     >  2:1589536/1‑148 (MQ=255)
                                                                                                                                       ttCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTa  >  1:2791504/1‑148 (MQ=255)
                                                                                                                                       ttCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTa  >  1:2039436/1‑148 (MQ=255)
                                                                                                                                          |                                                                                                                                                
GCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTA  >  CP009273/2884247‑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: