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,3850GT96.4% 589.3 / 11.8 192V291V (GTC→GTAcysJsulfite reductase, alpha subunit, flavoprotein
Reads supporting (aligned to +/- strand):  ref base G (5/2);  new base T (84/101);  total (89/103)
Fisher's exact test for biased strand distribution p-value = 2.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 192 total aligned reads displayed.

TTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAA  >  CP009273/2884241‑2884530
                                                                                                                                                |                                                                                                                                                 
ttCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTAcc                                                                                                                                                >  2:2027271/1‑147 (MQ=255)
ttCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACa                                                                                                                                                >  2:1161730/1‑148 (MQ=255)
ttCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACa                                                                                                                                                >  2:1779563/1‑148 (MQ=255)
ttCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACa                                                                                                                                                <  2:638262/148‑1 (MQ=255)
 tCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACaa                                                                                                                                               <  2:2858472/148‑1 (MQ=255)
   aCTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGt                                                                                                                                             <  2:2094191/148‑1 (MQ=255)
   aCTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGt                                                                                                                                             >  1:844752/1‑148 (MQ=255)
    cTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGtt                                                                                                                                            <  1:431929/148‑1 (MQ=255)
    cTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGtt                                                                                                                                            >  1:1268002/1‑148 (MQ=255)
    cTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGtt                                                                                                                                            >  2:1538681/1‑148 (MQ=255)
    cTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCAACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGtt                                                                                                                                            <  2:2278092/148‑1 (MQ=255)
     tGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTc                                                                                                                                           <  1:2321815/148‑1 (MQ=255)
      gCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCt                                                                                                                                          <  1:616816/148‑1 (MQ=255)
       cGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCtt                                                                                                                                         <  1:2075191/148‑1 (MQ=255)
        gggTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCttt                                                                                                                                        <  2:1902418/148‑1 (MQ=255)
        gggTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCttt                                                                                                                                        <  1:1564178/148‑1 (MQ=255)
          gTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCa                                                                                                                                      <  1:1749916/148‑1 (MQ=255)
           tAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCAc                                                                                                                                     <  1:2162434/148‑1 (MQ=255)
            aaGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCAcc                                                                                                                                    <  1:2698886/148‑1 (MQ=255)
             aGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCa                                                                                                                                   >  1:1143364/1‑148 (MQ=255)
                gTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTg                                                                                                                                <  1:3318762/148‑1 (MQ=255)
                     gTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGa                                                                                                                           <  2:1060699/148‑1 (MQ=255)
                        aTTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCg                                                                                                                        >  1:1625245/1‑148 (MQ=255)
                        aTTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCg                                                                                                                        <  1:555983/148‑1 (MQ=255)
                          tctcAACAATGTTGGCGGTGTTGACTGTCAGTTCGAAGTGAGACTGTAGCGCTTCGTTCAGAGGCATCGTTTTGCCCTCGACGGTGACAGGTTAATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCTGTTCCGGATCGtt                                                                                                                      <  1:2180201/148‑1 (MQ=255)
                           ctcAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCTTTc                                                                                                                     >  2:2951200/1‑148 (MQ=255)
                           ctcAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTc                                                                                                                     >  2:221683/1‑148 (MQ=255)
                             caacaaTGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTg                                                                                                                   <  1:3370914/148‑1 (MQ=255)
                              aacaaTGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGa                                                                                                                  <  2:1705227/148‑1 (MQ=255)
                                caaTGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATa                                                                                                                >  1:3539437/1‑148 (MQ=255)
                                    gTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAg                                                                                                            >  1:1864559/1‑148 (MQ=255)
                                     ttGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGa                                                                                                           <  1:722802/148‑1 (MQ=255)
                                       ggCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACg                                                                                                         >  2:1465166/1‑148 (MQ=255)
                                          ggTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGccc                                                                                                      <  1:1392416/148‑1 (MQ=255)
                                              ttGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAgcg                                                                                                  >  2:2965411/1‑148 (MQ=255)
                                               tGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAgcgc                                                                                                 >  2:629857/1‑148 (MQ=255)
                                                     tCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCAcc                                                                                           >  2:2677855/1‑148 (MQ=255)
                                                         ttCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGc                                                                                       <  1:525466/148‑1 (MQ=255)
                                                          tCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCt                                                                                      <  2:2646148/148‑1 (MQ=255)
                                                              aGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTa                                                                                  <  2:1986456/148‑1 (MQ=255)
                                                               gTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTTACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTaa                                                                                 >  1:3352750/1‑148 (MQ=255)
                                                               gTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTaa                                                                                 >  2:3121961/1‑148 (MQ=255)
                                                                  ccACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGc                                                                              <  2:67790/148‑1 (MQ=255)
                                                                       gTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGcc                                                                         >  1:2551463/1‑148 (MQ=255)
                                                                        tAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGccc                                                                        >  1:3170861/1‑148 (MQ=255)
                                                                          gcgcTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGa                                                                      >  1:375291/1‑148 (MQ=255)
                                                                           cgcTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAg                                                                     <  1:3082831/148‑1 (MQ=255)
                                                                            gcTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGt                                                                    <  2:1888365/148‑1 (MQ=255)
                                                                            gcTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGt                                                                    <  2:440459/148‑1 (MQ=255)
                                                                              ttCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCa                                                                  >  1:1320140/1‑148 (MQ=255)
                                                                              ttCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCa                                                                  >  2:382238/1‑148 (MQ=255)
                                                                                  ttCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTa                                                              >  1:3094811/1‑148 (MQ=255)
                                                                                    cAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAg                                                            >  2:3382545/1‑148 (MQ=255)
                                                                                    cAGAGGCAACGTTTTGCCCTAGACGGTGACAGGTTCATCGCCTTTGAGACACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTTATAACAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAg                                                            <  2:822513/148‑1 (MQ=255)
                                                                                     agagGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGt                                                           >  2:1127093/1‑148 (MQ=255)
                                                                                         gCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAAt                                                       >  2:2233504/1‑148 (MQ=255)
                                                                                         gCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAAt                                                       <  2:1250521/148‑1 (MQ=255)
                                                                                         gCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAAt                                                       <  2:1706652/148‑1 (MQ=255)
                                                                                         gCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAAt                                                       <  2:1706632/148‑1 (MQ=255)
                                                                                             cGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCg                                                   <  2:3566266/148‑1 (MQ=255)
                                                                                              gTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGa                                                  <  2:3201604/148‑1 (MQ=255)
                                                                                               ttttGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGat                                                 <  2:2610114/148‑1 (MQ=255)
                                                                                                 ttGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGatat                                               <  1:604134/148‑1 (MQ=255)
                                                                                                         gACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACg                                       >  1:1298891/1‑148 (MQ=255)
                                                                                                         gACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACg                                       >  1:3073261/1‑148 (MQ=255)
                                                                                                          aCGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGt                                      >  2:205113/1‑148 (MQ=255)
                                                                                                            ggTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCt                                    <  1:2272540/148‑1 (MQ=255)
                                                                                                             gTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCtt                                   <  1:2090267/148‑1 (MQ=255)
                                                                                                             gTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCtt                                   <  1:2091027/148‑1 (MQ=255)
                                                                                                              tGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCttt                                  <  1:1524260/148‑1 (MQ=255)
                                                                                                              tGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCttt                                  >  2:1524260/1‑148 (MQ=255)
                                                                                                              tGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCttt                                  <  1:327796/148‑1 (MQ=255)
                                                                                                              tGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCttt                                  >  2:327796/1‑148 (MQ=255)
                                                                                                               gACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCtttt                                 >  1:2331595/1‑148 (MQ=255)
                                                                                                                aCAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCttttt                                >  1:1300840/1‑148 (MQ=255)
                                                                                                                   ggTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCag                             <  2:1822780/148‑1 (MQ=255)
                                                                                                                   ggTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCag                             >  2:1117292/1‑148 (MQ=255)
                                                                                                                      tCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGt                          <  2:1159512/148‑1 (MQ=255)
                                                                                                                        aTCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTa                        <  1:2382080/148‑1 (MQ=255)
                                                                                                                         tCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTAc                       >  1:1096529/1‑148 (MQ=255)
                                                                                                                         tCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTAc                       >  1:1096464/1‑148 (MQ=255)
                                                                                                                         tCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTATTTCAGAGTTAc                       >  2:2787012/1‑148 (MQ=255)
                                                                                                                            ccTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGcc                    <  2:2407202/148‑1 (MQ=255)
                                                                                                                               ttCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCgg                 >  1:2453800/1‑148 (MQ=255)
                                                                                                                                 cAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTa               <  1:205113/148‑1 (MQ=255)
                                                                                                                                 cAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTa               <  1:629857/148‑1 (MQ=255)
                                                                                                                                   gCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAAt             <  2:2794551/148‑1 (MQ=255)
                                                                                                                                    ccACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAAtt            <  2:155900/148‑1 (MQ=255)
                                                                                                                                    ccACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAAtt            >  1:1730837/1‑148 (MQ=255)
                                                                                                                                      acaGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAAtttt          >  2:1309870/1‑148 (MQ=255)
                                                                                                                                      acaGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAAtttt          >  2:1310209/1‑148 (MQ=255)
                                                                                                                                      acaGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAAtttt          <  2:3073788/148‑1 (MQ=255)
                                                                                                                                       cagcagTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTc         <  1:1608386/148‑1 (MQ=255)
                                                                                                                                        agcagTTCTACAAGGTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCt        <  2:3459413/148‑1 (MQ=255)
                                                                                                                                         gcagTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTg       <  1:791838/148‑1 (MQ=255)
                                                                                                                                         gcagTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTg       >  1:862259/1‑148 (MQ=255)
                                                                                                                                             ttCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTAACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTa   <  1:3028200/148‑1 (MQ=255)
                                                                                                                                             ttCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTAACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTa   <  1:3027843/148‑1 (MQ=255)
                                                                                                                                              tCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTaa  >  1:2433131/1‑148 (MQ=255)
                                                                                                                                              tCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTaa  >  1:2433127/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
TTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAA  >  CP009273/2884241‑2884530

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: