Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,884,185 G→T 13.4% T358K (ACG→AAG)  cysJ ← sulfite reductase, alpha subunit, flavoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,884,1850GT13.4% 1040.8 / 91.2 321T358K (ACG→AAG) cysJsulfite reductase, alpha subunit, flavoprotein
Reads supporting (aligned to +/- strand):  ref base G (131/146);  new base T (19/24);  total (150/171)
Fisher's exact test for biased strand distribution p-value = 7.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

Only 100 of 333 total aligned reads displayed.

CACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCA  >  CP009273/2884037‑2884332
                                                                                                                                                    |                                                                                                                                                   
cacaccaACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGcg                                                                                                                                                     <  2:3132998/149‑1 (MQ=255)
  caccaACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGcgtc                                                                                                                                                   <  2:1972295/149‑1 (MQ=255)
    ccaACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGcgtcgt                                                                                                                                                 <  2:2697775/149‑1 (MQ=255)
     caACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCgtcgtc                                                                                                                                                >  2:211942/1‑149 (MQ=255)
       aCGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTcgc                                                                                                                                              <  2:1790045/149‑1 (MQ=255)
         ggTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCgcgg                                                                                                                                            >  1:1000291/1‑149 (MQ=255)
          gTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCgcggc                                                                                                                                           >  2:2283554/1‑149 (MQ=255)
              cGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTaa                                                                                                                                       <  2:2082027/149‑1 (MQ=255)
                   aCTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTg                                                                                                                                  <  2:2891377/149‑1 (MQ=255)
                     ttCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTa                                                                                                                                <  1:2973273/149‑1 (MQ=255)
                      tCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTaa                                                                                                                               >  1:4716741/1‑149 (MQ=255)
                      tCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTaa                                                                                                                               <  2:184090/149‑1 (MQ=255)
                       cGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAAc                                                                                                                              <  1:2754409/149‑1 (MQ=255)
                        gTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACt                                                                                                                             <  1:4640404/149‑1 (MQ=255)
                           ctcGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGAAACTTCg                                                                                                                          >  1:4093570/1‑149 (MQ=255)
                            tcGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGc                                                                                                                         <  2:1135100/149‑1 (MQ=255)
                              gACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCtt                                                                                                                       <  2:1479774/149‑1 (MQ=255)
                               aCTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCttt                                                                                                                      >  1:187096/1‑149 (MQ=255)
                                cgtCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCTTCGTCGCGGCGTAATGCTGTAACTTCGCTTTa                                                                                                                     <  1:1543041/147‑1 (MQ=255)
                                cTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCTTCGTCGCGGCGTAATGCTGTAACTTCGCTTTa                                                                                                                     >  2:4256293/1‑149 (MQ=255)
                                cTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTa                                                                                                                     >  2:2876270/1‑149 (MQ=255)
                                 ttCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTAt                                                                                                                    <  1:2234821/149‑1 (MQ=255)
                                  tCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCTTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATc                                                                                                                   >  2:437493/1‑149 (MQ=255)
                                    cGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGc                                                                                                                 >  1:3754676/1‑149 (MQ=255)
                                     gCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGcc                                                                                                                <  1:1452826/149‑1 (MQ=255)
                                      ccTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGccc                                                                                                               <  2:2211625/149‑1 (MQ=255)
                                       cTGAGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCTAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCTTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCcca                                                                                                              <  1:2660528/149‑1 (MQ=255)
                                          cGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCccacca                                                                                                           <  1:48517/149‑1 (MQ=255)
                                           gaggagGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCTTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAg                                                                                                          <  2:1281965/149‑1 (MQ=255)
                                            aggaggCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGc                                                                                                         <  2:4097/149‑1 (MQ=255)
                                            aggaggCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGc                                                                                                         <  2:3911/149‑1 (MQ=255)
                                              gaggCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCgg                                                                                                       >  2:3447336/1‑149 (MQ=255)
                                                    aTGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGcagcag                                                                                                 >  1:1667513/1‑149 (MQ=255)
                                                     tGGAATACAGTCGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTTACTTCGCTTTATCGCCCACCAGCGGCAGCAgt                                                                                                <  2:4180821/149‑1 (MQ=255)
                                                     tGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAgt                                                                                                >  2:629349/1‑149 (MQ=255)
                                                      ggAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCAGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAgtg                                                                                               >  1:2808011/1‑149 (MQ=255)
                                                       gAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAgtgt                                                                                              >  1:1149448/1‑149 (MQ=255)
                                                        aaTACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGtt                                                                                             >  2:1726796/1‑149 (MQ=255)
                                                           aCAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCa                                                                                          <  1:4206900/149‑1 (MQ=255)
                                                            cAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCAc                                                                                         <  1:3361720/149‑1 (MQ=255)
                                                            cAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCAc                                                                                         >  2:1111226/1‑149 (MQ=255)
                                                            cAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCGTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCAc                                                                                         <  2:552284/149‑1 (MQ=255)
                                                             aGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACt                                                                                        <  2:3088943/149‑1 (MQ=255)
                                                                   ggcgTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGt                                                                                  <  2:4227604/149‑1 (MQ=255)
                                                                    gcgTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTa                                                                                 >  1:1479542/1‑149 (MQ=255)
                                                                       tCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGc                                                                              >  2:3706877/1‑149 (MQ=255)
                                                                        cAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCACCAAGCTGTGCCGGGGAGAAACGCACCATGACAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCg                                                                             <  2:3365422/149‑1 (MQ=255)
                                                                          gCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCACGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTg                                                                           <  1:1256773/149‑1 (MQ=255)
                                                                          gCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCTTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTg                                                                           >  2:4140701/1‑149 (MQ=255)
                                                                            gggCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCTTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGc                                                                         <  2:1446844/149‑1 (MQ=255)
                                                                              gcgcAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGt                                                                       >  2:2464117/1‑149 (MQ=255)
                                                                                gcagcaGATTACTTCGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTaa                                                                     <  1:4521177/149‑1 (MQ=255)
                                                                                 cagcagATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAAt                                                                    <  2:3476970/149‑1 (MQ=255)
                                                                                  agcagATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAAtt                                                                   <  1:3617910/149‑1 (MQ=255)
                                                                                   gcagATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATtc                                                                  <  2:1252926/149‑1 (MQ=255)
                                                                                   gcagATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATtc                                                                  <  1:3057188/149‑1 (MQ=255)
                                                                                   gcagATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATtc                                                                  <  1:3578412/149‑1 (MQ=255)
                                                                                     agATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATtctc                                                                <  1:248237/149‑1 (MQ=255)
                                                                                       attaattaGAGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTcaa                                                              <  1:2135517/149‑1 (MQ=255)
                                                                                       attaagtaGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTcaa                                                              <  2:344491/149‑1 (MQ=255)
                                                                                        ttaattaGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCTTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTcaac                                                             <  1:3171462/149‑1 (MQ=255)
                                                                                        ttaattaGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTcaac                                                             <  2:685054/149‑1 (MQ=255)
                                                                                          aattaGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTCATTCTcaacaa                                                           >  2:2615096/1‑149 (MQ=255)
                                                                                          aatgaGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTcaacaa                                                           <  2:3632767/149‑1 (MQ=255)
                                                                                           attaGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAAt                                                          >  2:4696986/1‑149 (MQ=255)
                                                                                           attaGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAAt                                                          <  1:826917/149‑1 (MQ=255)
                                                                                           attaGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAAt                                                          <  2:2973369/149‑1 (MQ=255)
                                                                                            ttaGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCTTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATg                                                         >  2:32341/1‑149 (MQ=255)
                                                                                                     cGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTATCAACAATGTTGGCGgtg                                                <  2:1552514/149‑1 (MQ=255)
                                                                                                      ggCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCTTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGgtgt                                               >  1:2606378/1‑149 (MQ=255)
                                                                                                      ggCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGgtgt                                               <  1:4438749/149‑1 (MQ=255)
                                                                                                         aTCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGa                                            >  2:1345560/1‑149 (MQ=255)
                                                                                                          tCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCTTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGAc                                           <  2:3028740/149‑1 (MQ=255)
                                                                                                          tCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGAc                                           <  2:3304022/149‑1 (MQ=255)
                                                                                                          tCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGAc                                           >  2:4440588/1‑149 (MQ=255)
                                                                                                            aaGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACgg                                         <  1:3882810/149‑1 (MQ=255)
                                                                                                              gCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTc                                       <  1:4506521/149‑1 (MQ=255)
                                                                                                              gCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTc                                       <  2:4500781/149‑1 (MQ=255)
                                                                                                              gCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTc                                       >  2:3651854/1‑149 (MQ=255)
                                                                                                               cTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCTTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCa                                      >  1:724302/1‑149 (MQ=255)
                                                                                                                 gtgCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGt                                    >  2:888365/1‑149 (MQ=255)
                                                                                                                 gtgCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGt                                    <  2:2580419/149‑1 (MQ=255)
                                                                                                                 gtgCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGt                                    <  2:187722/149‑1 (MQ=255)
                                                                                                                   gCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTc                                  <  2:1000291/149‑1 (MQ=255)
                                                                                                                          agaAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGc                           <  1:4066151/149‑1 (MQ=255)
                                                                                                                           gaAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGcc                          >  2:3775337/1‑149 (MQ=255)
                                                                                                                             aaCGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCAc                        <  2:4669781/149‑1 (MQ=255)
                                                                                                                                gCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGt                     <  1:3482620/149‑1 (MQ=255)
                                                                                                                                 cACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTa                    <  2:2517375/149‑1 (MQ=255)
                                                                                                                                  aCCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAg                   >  2:4317191/1‑149 (MQ=255)
                                                                                                                                     aTGTCAACAATCGGCTTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAgcgc                <  2:2372779/149‑1 (MQ=255)
                                                                                                                                       gTCATCAATCGGCTTCGTCGCGGCGTAATGCTGTATCTTCGCTTTATCGCCCACCAGCGGCAGAAGTGTTTCTCTGCTGGTAAGCGTGTCGTAATTCTCAACAATGTTTGCGGTGTTGACGGTCTGTTCGAAGTGCCACTGTAGCGCtt              <  1:1678147/149‑1 (MQ=255)
                                                                                                                                       gTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCtt              <  1:2396156/149‑1 (MQ=255)
                                                                                                                                       gTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCtt              >  1:1047939/1‑149 (MQ=255)
                                                                                                                                        tCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTc             <  2:3502240/149‑1 (MQ=255)
                                                                                                                                         caacaaTCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCg            >  2:1943213/1‑149 (MQ=255)
                                                                                                                                          aacaaTCGGCTTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGt           <  1:3501117/149‑1 (MQ=255)
                                                                                                                                          aacaaTCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGt           <  1:682152/149‑1 (MQ=255)
                                                                                                                                           acaaTCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCGCCAGTGTTAAGATGGGGGTGAGCTTGGCATATGTGTAGACAATGTTGGCGGTGATGACTGTAAGTTCGAAGTGCAAGTGTAGCGAGTCGtt          >  2:3816284/1‑149 (MQ=255)
                                                                                                                                              aTCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCag       >  1:974289/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCA  >  CP009273/2884037‑2884332

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: