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.2% 822.4 / ‑2.2 250V291V (GTC→GTAcysJsulfite reductase, alpha subunit, flavoprotein
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base T (105/143);  total (106/144)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01
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 250 total aligned reads displayed.

TGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAAC  >  CP009273/2884238‑2884531
                                                                                                                                                   |                                                                                                                                                  
 gtTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACa                                                                                                                                                 <  1:1532757/150‑1 (MQ=255)
     cACTGCGTGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGtt                                                                                                                                             <  2:3463202/150‑1 (MQ=255)
      aCTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTc                                                                                                                                            >  2:3834688/1‑150 (MQ=255)
       cTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCt                                                                                                                                           <  2:4000369/150‑1 (MQ=255)
        tGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCtt                                                                                                                                          <  1:283771/150‑1 (MQ=255)
              tAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCa                                                                                                                                    <  1:3269749/150‑1 (MQ=255)
              tAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCa                                                                                                                                    >  1:3926002/1‑150 (MQ=255)
                 gCGTGGCGTAATTCTCAACACTGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTg                                                                                                                                 <  2:2879738/150‑1 (MQ=255)
                  cGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCAGCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGc                                                                                                                                >  2:854801/1‑150 (MQ=255)
                   gTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGcc                                                                                                                               <  1:4718099/150‑1 (MQ=255)
                    tGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCg                                                                                                                              <  1:293418/150‑1 (MQ=255)
                      gCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGa                                                                                                                            <  2:2415676/150‑1 (MQ=255)
                        gTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATc                                                                                                                          <  1:4424292/150‑1 (MQ=255)
                        gTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATc                                                                                                                          >  1:3475929/1‑150 (MQ=255)
                         tAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCg                                                                                                                         <  1:573809/150‑1 (MQ=255)
                          aaTTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGt                                                                                                                        >  1:2113318/1‑150 (MQ=255)
                          aaTTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGt                                                                                                                        >  1:3724805/1‑150 (MQ=255)
                           aTTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCTGTTCTACAAGTTCTTTCACCAGTGCCGGATCGtt                                                                                                                       <  1:3375330/150‑1 (MQ=255)
                           aTTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGtt                                                                                                                       >  2:2443293/1‑150 (MQ=255)
                            ttCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTc                                                                                                                      <  1:2959793/150‑1 (MQ=255)
                              ctcAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTg                                                                                                                    >  1:5264937/1‑150 (MQ=255)
                               tcAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGa                                                                                                                   >  1:1643503/1‑150 (MQ=255)
                                 aacaaTGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATa                                                                                                                 >  2:1588909/1‑150 (MQ=255)
                                   caaTGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATAcc                                                                                                               <  1:4296526/150‑1 (MQ=255)
                                   caaTGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATAcc                                                                                                               <  1:3425793/150‑1 (MQ=255)
                                    aaTGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCa                                                                                                              >  1:1476120/1‑150 (MQ=255)
                                           gCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGccc                                                                                                       <  2:3299822/150‑1 (MQ=255)
                                            cGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCTCTTCGTTCAGAGGCAACGTTTTGACCTAGACGGTGACAGGTTCATCTCCTTTCAGCCACATAAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTTATAACAGACTCCCa                                                                                                      <  2:2249768/150‑1 (MQ=255)
                                            cGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCa                                                                                                      <  1:3822294/150‑1 (MQ=255)
                                             ggTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAg                                                                                                     <  2:5190575/150‑1 (MQ=255)
                                             ggTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAg                                                                                                     >  2:5004135/1‑150 (MQ=255)
                                             ggTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAg                                                                                                     <  2:4897325/150‑1 (MQ=255)
                                             ggTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAg                                                                                                     >  1:1396060/1‑150 (MQ=255)
                                              gtgtTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAgc                                                                                                    <  1:2345878/150‑1 (MQ=255)
                                               tgtTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAgcg                                                                                                   >  1:452589/1‑150 (MQ=255)
                                                  tGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGt                                                                                                <  2:2617246/150‑1 (MQ=255)
                                                   gACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTc                                                                                               <  2:3606996/150‑1 (MQ=255)
                                                   gACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTc                                                                                               >  1:673148/1‑150 (MQ=255)
                                                      ggTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCAcc                                                                                            >  2:1092342/1‑150 (MQ=255)
                                                       gTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCAccc                                                                                           >  2:2568955/1‑150 (MQ=255)
                                                          aGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGc                                                                                        >  1:1093474/1‑150 (MQ=255)
                                                          aGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGAAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGc                                                                                        >  1:1948271/1‑150 (MQ=255)
                                                              cGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGt                                                                                    >  1:3544520/1‑150 (MQ=255)
                                                              cGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGt                                                                                    <  1:1258742/150‑1 (MQ=255)
                                                              cGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGt                                                                                    <  2:206957/150‑1 (MQ=255)
                                                                 aGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAAc                                                                                 <  1:2953828/150‑1 (MQ=255)
                                                                         tGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGccc                                                                         <  1:4055318/150‑1 (MQ=255)
                                                                         tGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGccc                                                                         <  1:4507834/150‑1 (MQ=255)
                                                                          gTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCg                                                                        <  1:2534289/150‑1 (MQ=255)
                                                                            aGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAg                                                                      >  1:4202987/1‑150 (MQ=255)
                                                                             gcgcTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGt                                                                     <  1:1784564/150‑1 (MQ=255)
                                                                               gcTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCa                                                                   <  2:1250034/150‑1 (MQ=255)
                                                                                  tCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACCGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCt                                                                <  2:2994669/150‑1 (MQ=255)
                                                                                   cGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTa                                                               <  2:2367660/150‑1 (MQ=255)
                                                                                     ttCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAg                                                             <  1:4578687/150‑1 (MQ=255)
                                                                                      tCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCGCCCGGCTGGTAACGCATGCCCGAGTCACCTAAGt                                                            <  2:682994/150‑1 (MQ=255)
                                                                                        agagGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCa                                                          >  2:195663/1‑150 (MQ=255)
                                                                                           ggCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAAtt                                                       >  2:827108/1‑150 (MQ=255)
                                                                                           ggCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAAtt                                                       >  1:3172652/1‑150 (MQ=255)
                                                                                            gCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAAttt                                                      <  1:4342552/150‑1 (MQ=255)
                                                                                            gCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAAttt                                                      <  1:1563214/150‑1 (MQ=255)
                                                                                             cAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTc                                                     <  1:792878/150‑1 (MQ=255)
                                                                                             cAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTc                                                     <  2:3164968/150‑1 (MQ=255)
                                                                                             cAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTc                                                     >  1:2983750/1‑150 (MQ=255)
                                                                                                cGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCgtt                                                  <  1:2239149/150‑3 (MQ=255)
                                                                                                   tttGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATg                                               <  1:2907779/150‑1 (MQ=255)
                                                                                                   tttGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATg                                               >  2:4799046/1‑150 (MQ=255)
                                                                                                       cccTCGACGGTGACAGGTTCATCGCATTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGa                                           <  2:1862914/150‑1 (MQ=255)
                                                                                                        ccTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGaa                                          >  2:1587458/1‑150 (MQ=255)
                                                                                                        ccTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGaa                                          <  2:3475929/150‑1 (MQ=255)
                                                                                                            gACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTc                                      >  2:2712376/1‑150 (MQ=255)
                                                                                                            gACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTc                                      >  1:3344117/1‑150 (MQ=255)
                                                                                                            gACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTc                                      >  2:3192593/1‑150 (MQ=255)
                                                                                                               ggTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCttt                                   <  2:947170/150‑1 (MQ=255)
                                                                                                               ggTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCttt                                   >  2:4343959/1‑150 (MQ=255)
                                                                                                               ggTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCttt                                   <  2:3121192/150‑1 (MQ=255)
                                                                                                                gTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCtttt                                  >  2:5310869/1‑150 (MQ=255)
                                                                                                                gTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCtttt                                  <  1:903211/150‑1 (MQ=255)
                                                                                                                gTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCtttt                                  <  2:4574588/150‑1 (MQ=255)
                                                                                                                  gACAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTc                                <  2:3691955/150‑1 (MQ=255)
                                                                                                                   aCAGGTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCa                               <  1:1593057/150‑1 (MQ=255)
                                                                                                                       gTTCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGt                           <  1:4059386/150‑1 (MQ=255)
                                                                                                                        ttCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGtt                          <  1:1864694/150‑1 (MQ=255)
                                                                                                                        ttCATCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGtt                          >  2:4769465/1‑150 (MQ=255)
                                                                                                                            tCGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGc                      >  1:4883354/1‑150 (MQ=255)
                                                                                                                             cGCCTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGcc                     <  2:3352380/150‑1 (MQ=255)
                                                                                                                                cTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCgg                  >  2:1740154/1‑150 (MQ=255)
                                                                                                                                cTTTCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCgg                  >  1:187011/1‑150 (MQ=255)
                                                                                                                                   tCAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTaa               <  1:5193893/150‑1 (MQ=255)
                                                                                                                                    cAGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGAGACCAGACGCCCAGCGCGTCACCCGGCTGGGAACGCAGGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAAt              <  2:1133397/150‑1 (MQ=255)
                                                                                                                                     aGCCACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAAtt             <  2:4625454/150‑1 (MQ=255)
                                                                                                                                       ccACAGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAAtttt           >  1:176431/1‑150 (MQ=255)
                                                                                                                                        cacaGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTc          >  2:3577858/1‑150 (MQ=255)
                                                                                                                                         acaGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCt         >  2:5201505/1‑150 (MQ=255)
                                                                                                                                         acaGCAGTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCt         <  1:5004135/150‑1 (MQ=255)
                                                                                                                                              agTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTa    >  2:1339079/1‑150 (MQ=255)
                                                                                                                                              agTTCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTa    >  2:917979/1‑150 (MQ=255)
                                                                                                                                               tttCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTaa   >  1:3658260/2‑150 (MQ=255)
                                                                                                                                                ttCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAAc  <  2:4097704/150‑1 (MQ=255)
                                                                                                                                                ttCTACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGATAATTTTCTGGTTAAc  <  2:4411184/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAAC  >  CP009273/2884238‑2884531

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: