Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,936,836 C→A 100% A47A (GCC→GCA csdA → cysteine sulfinate desulfinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,936,8360CA99.6% 872.7 / ‑4.4 263A47A (GCC→GCAcsdAcysteine sulfinate desulfinase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (119/143);  total (119/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.96e-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 263 total aligned reads displayed.

GTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCCGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTGGACGCGCGGCACCACTGAATCCATCAACA  >  CP009273/2936691‑2936981
                                                                                                                                                 |                                                                                                                                                 
gTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGa                                                                                                                                                <  2:3250042/149‑1 (MQ=255)
gTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGa                                                                                                                                                >  1:2435626/1‑149 (MQ=255)
  aCCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGaaa                                                                                                                                              >  1:446001/1‑149 (MQ=255)
  aCCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGaaa                                                                                                                                              <  1:3475258/149‑1 (MQ=255)
  aCCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGaaa                                                                                                                                              >  1:3182173/1‑149 (MQ=255)
  aCCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGaaa                                                                                                                                              <  1:2832195/149‑1 (MQ=255)
    cATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTGACAGTCTGAGCGCAGGAAACg                                                                                                                                            >  2:120967/1‑149 (MQ=255)
       gAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTcc                                                                                                                                         >  1:92495/1‑149 (MQ=255)
        aaCGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCa                                                                                                                                        <  2:2131744/149‑1 (MQ=255)
            tttttAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGc                                                                                                                                    >  2:1652133/1‑149 (MQ=255)
               ttAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGc                                                                                                                                 <  2:3773281/149‑1 (MQ=255)
               ttAATCCCGCGCAGTATCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCAATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGc                                                                                                                                 >  2:1960933/1‑149 (MQ=255)
                   tCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGt                                                                                                                             <  2:3117248/149‑1 (MQ=255)
                   tCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGt                                                                                                                             >  1:1924096/1‑149 (MQ=255)
                    cccGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGGGCAGGAAACGTCCATCGCAGCCAGtt                                                                                                                            <  1:3829053/149‑1 (MQ=255)
                    cccGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGtt                                                                                                                            >  2:2093962/1‑149 (MQ=255)
                     ccGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGttt                                                                                                                           <  2:2133289/149‑1 (MQ=255)
                     ccGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGttt                                                                                                                           <  2:2653864/149‑1 (MQ=255)
                      cgcgCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTg                                                                                                                          >  2:2403541/1‑149 (MQ=255)
                       gcgcAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGc                                                                                                                         <  1:2976931/149‑1 (MQ=255)
                           aGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGaa                                                                                                                     <  2:3379248/149‑1 (MQ=255)
                             tttCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGc                                                                                                                   >  2:675287/1‑149 (MQ=255)
                             tttCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGc                                                                                                                   <  2:1538101/149‑1 (MQ=255)
                                 gcgcCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCaa                                                                                                               <  1:2698159/149‑1 (MQ=255)
                                 gcgcCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCaa                                                                                                               <  2:4009248/149‑1 (MQ=255)
                                  cgcCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAAc                                                                                                              <  1:349937/149‑1 (MQ=255)
                                  cgcCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAAc                                                                                                              <  1:2951056/149‑1 (MQ=255)
                                   gcCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACg                                                                                                             >  2:2280770/1‑149 (MQ=255)
                                     ccAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGcc                                                                                                           <  1:4053945/149‑1 (MQ=255)
                                       aGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTg                                                                                                         <  2:3914049/149‑1 (MQ=255)
                                       aGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTg                                                                                                         <  2:1930173/149‑1 (MQ=255)
                                        gTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGa                                                                                                        <  1:2010176/149‑1 (MQ=255)
                                         tttCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGAc                                                                                                       <  1:3740874/149‑1 (MQ=255)
                                          ttCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGAcc                                                                                                      >  2:1326121/1‑149 (MQ=255)
                                                 aCTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGtt                                                                                               <  2:319304/149‑1 (MQ=255)
                                                 aCTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGtt                                                                                               <  1:633597/149‑1 (MQ=255)
                                                  cTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTa                                                                                              >  2:3143433/1‑149 (MQ=255)
                                                  cTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTa                                                                                              <  1:2284483/149‑1 (MQ=255)
                                                   tACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTAt                                                                                             <  1:1519105/149‑1 (MQ=255)
                                                      aGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGaa                                                                                          >  1:3771989/1‑149 (MQ=255)
                                                       ggATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAg                                                                                         >  2:3000470/1‑149 (MQ=255)
                                                       ggATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAg                                                                                         <  1:3004070/149‑1 (MQ=255)
                                                          tGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTg                                                                                      <  2:2747252/149‑1 (MQ=255)
                                                           gCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGc                                                                                     <  2:446001/149‑1 (MQ=255)
                                                               gCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACga                                                                                 <  2:220395/149‑1 (MQ=255)
                                                               gCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACga                                                                                 <  1:2857745/149‑1 (MQ=255)
                                                                cGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACgag                                                                                <  1:2336142/149‑1 (MQ=255)
                                                                cGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACgag                                                                                <  2:183242/149‑1 (MQ=255)
                                                                  tCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACgagag                                                                              >  2:3498414/1‑149 (MQ=255)
                                                                  tCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACgagag                                                                              >  2:2635166/1‑149 (MQ=255)
                                                                  tCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACgagag                                                                              >  1:3316966/1‑149 (MQ=255)
                                                                       ctcGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGt                                                                         <  2:2970622/149‑1 (MQ=255)
                                                                            cAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGgcgc                                                                    <  1:1573145/149‑1 (MQ=255)
                                                                            cAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGgcgc                                                                    <  1:1207225/149‑1 (MQ=255)
                                                                              gcgcCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCaa                                                                  >  2:3831704/1‑149 (MQ=255)
                                                                              gcgcCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCaa                                                                  >  1:2311534/1‑149 (MQ=255)
                                                                              gcgcCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCaa                                                                  >  1:2175928/1‑149 (MQ=255)
                                                                               cgccgcGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAAt                                                                 <  2:491867/149‑1 (MQ=255)
                                                                                   gcgACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACt                                                             >  2:2896179/1‑149 (MQ=255)
                                                                                      aCCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGaa                                                          >  1:2599928/1‑149 (MQ=255)
                                                                                       ccGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAAt                                                         >  1:2267914/1‑149 (MQ=255)
                                                                                           gcTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCAc                                                     <  1:1082445/149‑1 (MQ=255)
                                                                                             ttAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCg                                                   <  2:910795/149‑1 (MQ=255)
                                                                                               aaaCCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGga                                                 <  1:1453463/149‑1 (MQ=255)
                                                                                               aaaCCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGga                                                 <  1:835642/149‑1 (MQ=255)
                                                                                               aaaCCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGga                                                 <  1:3585855/149‑1 (MQ=255)
                                                                                                aaCCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGgat                                                <  1:57619/149‑1 (MQ=255)
                                                                                                 aCCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGgatg                                               <  1:1977990/149‑1 (MQ=255)
                                                                                                 aCCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGgatg                                               >  1:3037630/1‑149 (MQ=255)
                                                                                                  ccTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGgatga                                              <  1:2320216/149‑1 (MQ=255)
                                                                                                   cTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGgatgat                                             <  2:784014/149‑1 (MQ=255)
                                                                                                   cTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGgatgat                                             <  1:2093962/149‑1 (MQ=255)
                                                                                                   cTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGgatgat                                             <  1:2116233/149‑1 (MQ=255)
                                                                                                     gAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATaa                                           >  2:1488212/1‑149 (MQ=255)
                                                                                                     gAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATaa                                           >  2:1034820/1‑149 (MQ=255)
                                                                                                      aaGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATaaa                                          <  2:2932839/149‑1 (MQ=255)
                                                                                                       aGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATaaaa                                         <  1:325927/149‑1 (MQ=255)
                                                                                                        gCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAAc                                        >  2:1384192/1‑149 (MQ=255)
                                                                                                        gCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAAc                                        >  2:3700167/1‑149 (MQ=255)
                                                                                                          cGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTa                                      <  2:2937150/149‑1 (MQ=255)
                                                                                                            tGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATc                                    <  2:2440954/149‑1 (MQ=255)
                                                                                                             ggTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCg                                   >  1:1999531/1‑149 (MQ=255)
                                                                                                               ttGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTc                                 <  2:1924096/149‑1 (MQ=255)
                                                                                                                tGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCt                                <  1:3405932/149‑1 (MQ=255)
                                                                                                                tGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCt                                <  2:2818033/149‑1 (MQ=255)
                                                                                                                tGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCt                                <  2:436804/149‑1 (MQ=255)
                                                                                                                 gAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTg                               >  1:1681037/1‑149 (MQ=255)
                                                                                                                   aGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTGGa                             <  2:1599854/149‑1 (MQ=255)
                                                                                                                       aCCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTGGAcgcg                         >  2:2181063/1‑149 (MQ=255)
                                                                                                                       aCCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTAGAcgcg                         >  2:2181045/1‑149 (MQ=255)
                                                                                                                        cccAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCGTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTGGAcgcgc                        <  1:1479849/149‑1 (MQ=255)
                                                                                                                                ttttACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTGGACGCGCGGcaccac                >  1:2788828/1‑149 (MQ=255)
                                                                                                                                ttttACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTGGACGCGCGGcaccac                >  2:2502112/1‑149 (MQ=255)
                                                                                                                                ttttACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTGGACGCGCGGcaccac                >  1:726645/1‑149 (MQ=255)
                                                                                                                                   tACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTGGACGCGCGGCACCACTGa             >  1:3079758/1‑149 (MQ=255)
                                                                                                                                   tACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTGGACGCGCGGCACCACTGa             >  1:1023056/1‑149 (MQ=255)
                                                                                                                                    aCAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTGGACGCGCGGCACCACTGaa            >  2:145227/1‑149 (MQ=255)
                                                                                                                                          tGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTGGACGCGCGGCACCACTGAATCCATc      >  2:788973/1‑149 (MQ=255)
                                                                                                                                            aGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTGGACGCGCGGCACCACTGAATCCATCaa    >  1:3348662/1‑149 (MQ=255)
                                                                                                                                             gcgcAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTGGACGCGCGGCACCACTGAATCCATCAAc   >  2:3266592/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
GTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCCGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTGGACGCGCGGCACCACTGAATCCATCAACA  >  CP009273/2936691‑2936981

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: