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,3850GT99.0% 678.4 / ‑1.6 206V291V (GTC→GTA) cysJsulfite reductase, alpha subunit, flavoprotein
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base T (98/106);  total (98/108)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 206 total aligned reads displayed.

TTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTA  >  CP009273/2884240‑2884529
                                                                                                                                                 |                                                                                                                                                
tttCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTAc                                                                                                                                                >  1:1061871/1‑148 (MQ=255)
   cACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAg                                                                                                                                             <  1:1089850/148‑1 (MQ=255)
    aCTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGt                                                                                                                                            <  1:569832/148‑1 (MQ=255)
     cTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGtt                                                                                                                                           <  1:714789/148‑1 (MQ=255)
      tGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTc                                                                                                                                          >  2:2446033/1‑148 (MQ=255)
      tGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTc                                                                                                                                          >  2:681844/1‑148 (MQ=255)
                cGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGt                                                                                                                                <  2:875842/148‑1 (MQ=255)
                   ggCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGATGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGcc                                                                                                                             <  2:2480287/148‑1 (MQ=255)
                     cGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCgg                                                                                                                           >  2:1959865/1‑148 (MQ=255)
                     cGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCgg                                                                                                                           >  2:1604346/1‑148 (MQ=255)
                     cGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCgg                                                                                                                           >  2:115985/1‑148 (MQ=255)
                       tAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGAt                                                                                                                         >  2:64295/1‑148 (MQ=255)
                       tAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGAt                                                                                                                         <  2:3345908/148‑1 (MQ=255)
                         aTTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCg                                                                                                                       >  2:2826274/1‑148 (MQ=255)
                           tctcAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGtt                                                                                                                     <  1:705090/148‑1 (MQ=255)
                            ctcAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTc                                                                                                                    <  2:690798/148‑1 (MQ=255)
                             tcAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCt                                                                                                                   >  1:2464005/1‑148 (MQ=255)
                              caacaaTGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTg                                                                                                                  >  1:498958/1‑148 (MQ=255)
                               aacaaTGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGa                                                                                                                 <  1:3155965/148‑1 (MQ=255)
                               aacaaTGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGa                                                                                                                 <  1:3368084/148‑1 (MQ=255)
                                acaaTGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGAt                                                                                                                >  2:2057848/1‑148 (MQ=255)
                                 caaTGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATa                                                                                                               >  1:365343/1‑148 (MQ=255)
                                  aaTGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATAc                                                                                                              >  2:925731/1‑148 (MQ=255)
                                   aTGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATAcc                                                                                                             >  1:1926811/1‑148 (MQ=255)
                                   aTGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATAcc                                                                                                             >  2:3383803/1‑148 (MQ=255)
                                     gTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGTTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAg                                                                                                           >  2:782296/1‑148 (MQ=255)
                                     gTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAg                                                                                                           <  2:889386/148‑1 (MQ=255)
                                     gTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAg                                                                                                           >  2:1605026/1‑148 (MQ=255)
                                     gTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAg                                                                                                           <  2:2338208/148‑1 (MQ=255)
                                     gTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAg                                                                                                           >  2:2543089/1‑148 (MQ=255)
                                     gTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAg                                                                                                           >  1:143976/1‑148 (MQ=255)
                                     gTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAg                                                                                                           >  1:3050025/1‑148 (MQ=255)
                                      ttGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGa                                                                                                          >  2:517751/1‑148 (MQ=255)
                                      ttGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGa                                                                                                          >  2:2283493/1‑148 (MQ=255)
                                         gCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGc                                                                                                       >  1:2333252/1‑148 (MQ=255)
                                         gCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGc                                                                                                       >  2:1116309/1‑148 (MQ=255)
                                          cGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGcc                                                                                                      <  1:1758826/148‑1 (MQ=255)
                                          cGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGcc                                                                                                      <  2:1505630/148‑1 (MQ=255)
                                             tgtTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAg                                                                                                   >  2:1655881/1‑148 (MQ=255)
                                             tgtTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAg                                                                                                   <  1:600587/148‑1 (MQ=255)
                                             tgtTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAg                                                                                                   >  2:2403493/1‑148 (MQ=255)
                                                  aCGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGt                                                                                              >  2:2152476/1‑148 (MQ=255)
                                                  aCGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGt                                                                                              <  2:1802415/148‑1 (MQ=255)
                                                   cGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTc                                                                                             >  1:2663586/1‑148 (MQ=255)
                                                    ggTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCa                                                                                            <  2:2781361/148‑1 (MQ=255)
                                                       cAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCAccc                                                                                         >  1:1151874/1‑148 (MQ=255)
                                                       cAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCAccc                                                                                         >  1:2421339/1‑148 (MQ=255)
                                                              aaGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGt                                                                                  >  2:819191/1‑148 (MQ=255)
                                                                gTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTaa                                                                                >  2:2101901/1‑148 (MQ=255)
                                                                  gCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACg                                                                              >  2:3491235/1‑148 (MQ=255)
                                                                   ccACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGc                                                                             <  2:1580348/148‑1 (MQ=255)
                                                                        gTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGcc                                                                        >  1:984914/1‑148 (MQ=255)
                                                                             gcTTCGTTCAGAGGCATCGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGNTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGt                                                                   <  1:1780005/148‑1 (MQ=255)
                                                                              cTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTc                                                                  <  1:2489891/148‑1 (MQ=255)
                                                                               ttCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCa                                                                 <  2:2695647/148‑1 (MQ=255)
                                                                                tCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCAc                                                                <  2:1751121/148‑1 (MQ=255)
                                                                                   ttCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTa                                                             >  2:2372948/1‑148 (MQ=255)
                                                                                   ttCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTa                                                             <  2:3446975/148‑1 (MQ=255)
                                                                                     cAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAg                                                           <  2:143976/148‑1 (MQ=255)
                                                                                      agagGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGt                                                          >  1:1198483/1‑148 (MQ=255)
                                                                                       gagGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTc                                                         >  2:1517134/1‑148 (MQ=255)
                                                                                            accGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAAttt                                                    >  2:3155382/3‑148 (MQ=255)
                                                                                            aaCGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAAttt                                                    <  1:673968/148‑1 (MQ=255)
                                                                                            aaCGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAAttt                                                    <  1:671519/148‑1 (MQ=255)
                                                                                            aaCGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAAttt                                                    >  2:761830/1‑148 (MQ=255)
                                                                                                ttttGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGat                                                <  2:231909/148‑1 (MQ=255)
                                                                                                ttttGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGat                                                >  1:2400616/1‑148 (MQ=255)
                                                                                                  ttGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGatat                                              >  1:1110484/1‑148 (MQ=255)
                                                                                                    gCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATgg                                            <  2:1018442/148‑1 (MQ=255)
                                                                                                    gCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATgg                                            <  1:1601202/148‑1 (MQ=255)
                                                                                                    gCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATgg                                            <  1:2786579/148‑1 (MQ=255)
                                                                                                       cTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGa                                         <  1:347484/148‑1 (MQ=255)
                                                                                                         cGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAAc                                       >  2:3042611/1‑148 (MQ=255)
                                                                                                          gACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACg                                      >  1:1377265/1‑148 (MQ=255)
                                                                                                             ggTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCt                                   <  1:2607160/148‑1 (MQ=255)
                                                                                                             ggTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCt                                   <  2:2595041/148‑1 (MQ=255)
                                                                                                             ggTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCt                                   <  2:365343/148‑1 (MQ=255)
                                                                                                                    ggTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCag                            <  2:2743597/148‑1 (MQ=255)
                                                                                                                    ggTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCag                            <  2:3049857/148‑1 (MQ=255)
                                                                                                                       tCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGt                         <  1:1835192/148‑1 (MQ=255)
                                                                                                                         aTCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTa                       <  2:1075754/148‑1 (MQ=255)
                                                                                                                         aTCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTa                       >  1:2953180/1‑148 (MQ=255)
                                                                                                                           cGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACg                     >  2:2257474/1‑148 (MQ=255)
                                                                                                                              cTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGccc                  >  2:931511/1‑148 (MQ=255)
                                                                                                                               tttCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCg                 >  2:1416136/1‑148 (MQ=255)
                                                                                                                                 tCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGt               >  1:2366091/1‑148 (MQ=255)
                                                                                                                                  cAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTa              <  1:199667/148‑1 (MQ=255)
                                                                                                                                  cAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTa              <  1:1631226/148‑1 (MQ=255)
                                                                                                                                       acaGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGCTATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAAtttt         <  1:1258292/148‑1 (MQ=255)
                                                                                                                                       acaGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAAtttt         >  1:3210075/1‑148 (MQ=255)
                                                                                                                                        cagcagTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTc        <  2:1232565/148‑1 (MQ=255)
                                                                                                                                         agcagTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCt       <  2:1151759/148‑1 (MQ=255)
                                                                                                                                          gcagTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTg      >  1:2993316/1‑148 (MQ=255)
                                                                                                                                           cagTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTgg     >  2:2992607/1‑148 (MQ=255)
                                                                                                                                           cagTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTgg     <  1:2999212/148‑1 (MQ=255)
                                                                                                                                            agTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGt    <  1:651599/148‑1 (MQ=255)
                                                                                                                                            agTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGt    <  1:1116309/148‑1 (MQ=255)
                                                                                                                                             gTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGtt   >  2:1301769/1‑148 (MQ=255)
                                                                                                                                              ttCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTCATTTTCTGGTTa  >  2:2762598/1‑148 (MQ=255)
                                                                                                                                              ttCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTa  <  1:886988/148‑1 (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: ‑