Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,882,054 C→A 7.0% D469Y (GAT→TAT)  cysI ← sulfite reductase, beta subunit, NAD(P)‑binding, heme‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,882,0540CA 7.0% 637.8 / 16.5 171D469Y (GAT→TAT) cysIsulfite reductase, beta subunit, NAD(P)‑binding, heme‑binding
Reads supporting (aligned to +/- strand):  ref base C (75/84);  new base A (5/7);  total (80/91)
Fisher's exact test for biased strand distribution p-value = 7.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 172 total aligned reads displayed.

ATCCGTGGGATACGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAA  >  CP009273/2881909‑2882201
                                                                                                                                                 |                                                                                                                                                   
aTCCGTGGGATACGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGc                                                                                                                                                   <  2:2028878/148‑1 (MQ=255)
aTCCGTGGGATACGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGc                                                                                                                                                   <  1:824978/148‑1 (MQ=255)
 tCCGTGGGATACGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCt                                                                                                                                                  >  2:1920165/1‑148 (MQ=255)
 tCCGTGGGATACGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCt                                                                                                                                                  >  1:177446/1‑148 (MQ=255)
  ccGTGGGATACGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTg                                                                                                                                                 <  1:2318581/148‑1 (MQ=255)
   cGTGGGATACGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGa                                                                                                                                                <  1:3119972/148‑1 (MQ=255)
       ggATACGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACAcc                                                                                                                                            <  1:908433/148‑1 (MQ=255)
        gATACGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCa                                                                                                                                           >  2:1396810/1‑148 (MQ=255)
        gATACGTGTCCCAATGCGGTTGCCGCAAAGATGCAGGTTGGAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCTGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCa                                                                                                                                           <  2:2656735/148‑1 (MQ=255)
          tACGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATg                                                                                                                                         <  1:1650094/148‑1 (MQ=255)
            cGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGtt                                                                                                                                       <  2:3157252/148‑1 (MQ=255)
             gtgtCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGttt                                                                                                                                      <  2:344620/148‑1 (MQ=255)
              tgtCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTc                                                                                                                                     >  1:1722437/1‑148 (MQ=255)
                tCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGc                                                                                                                                   <  1:717939/148‑1 (MQ=255)
                     aTGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTa                                                                                                                              >  2:3342236/1‑148 (MQ=255)
                     aTGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTa                                                                                                                              >  2:3342251/1‑148 (MQ=255)
                         ggTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAAtt                                                                                                                          <  1:2216755/148‑1 (MQ=255)
                         ggTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAAtt                                                                                                                          <  2:3263699/148‑1 (MQ=255)
                           ttGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTAt                                                                                                                        <  1:2743688/148‑1 (MQ=255)
                              ccgccAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGa                                                                                                                     >  1:2005311/1‑148 (MQ=255)
                              ccgccAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATAGCTGACACCATGTTTCGCCATTAAATTATCGa                                                                                                                     <  2:2821673/148‑1 (MQ=255)
                                gccAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGAtg                                                                                                                   >  2:2986736/1‑148 (MQ=255)
                                gccAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGAtg                                                                                                                   >  2:2986651/1‑148 (MQ=255)
                                  cAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGAtgtt                                                                                                                 <  1:2304720/148‑1 (MQ=255)
                                   aaGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGAtgttg                                                                                                                <  2:352163/148‑1 (MQ=255)
                                    aGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGAtgttgt                                                                                                               <  1:226606/148‑1 (MQ=255)
                                     gATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTc                                                                                                              >  1:2227473/1‑148 (MQ=255)
                                         cAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATa                                                                                                          >  1:2907838/1‑148 (MQ=255)
                                         cAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATACTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATa                                                                                                          >  2:1893378/1‑148 (MQ=255)
                                         cAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATACTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATa                                                                                                          >  2:2135008/1‑148 (MQ=255)
                                           ggTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATaaa                                                                                                        <  1:1599234/148‑1 (MQ=255)
                                             ttGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAg                                                                                                      <  2:2137238/148‑1 (MQ=255)
                                              tGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGa                                                                                                     <  1:538496/148‑1 (MQ=255)
                                              tGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGa                                                                                                     <  2:1436180/148‑1 (MQ=255)
                                              tGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGa                                                                                                     <  2:1436115/148‑1 (MQ=255)
                                              tGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGa                                                                                                     <  1:1031459/148‑1 (MQ=255)
                                              tGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGa                                                                                                     <  1:2437387/148‑1 (MQ=255)
                                               gTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGAc                                                                                                    <  1:2234834/148‑1 (MQ=255)
                                                  gCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGc                                                                                                 >  1:3456866/1‑148 (MQ=255)
                                                   cGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCa                                                                                                >  2:1626271/1‑148 (MQ=255)
                                                     aCCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAgg                                                                                              >  2:3514875/1‑148 (MQ=255)
                                                      cccGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGa                                                                                             >  2:662307/1‑148 (MQ=255)
                                                       ccGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGaa                                                                                            <  2:1463641/148‑1 (MQ=255)
                                                        cGGCGCTTGACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGaaa                                                                                           <  2:1980595/148‑1 (MQ=255)
                                                            gcTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGct                                                                                       <  1:2753837/148‑1 (MQ=255)
                                                             cTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGctc                                                                                      >  1:66909/1‑148 (MQ=255)
                                                             cTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGctc                                                                                      >  2:1213669/1‑148 (MQ=255)
                                                              tttACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGctct                                                                                     >  1:123555/1‑148 (MQ=255)
                                                               ttACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTg                                                                                    <  2:1852196/148‑1 (MQ=255)
                                                                 aCCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCt                                                                                  >  1:920507/1‑148 (MQ=255)
                                                                    caccaGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTcc                                                                               >  1:1788712/1‑148 (MQ=255)
                                                                      ccaGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGc                                                                             >  2:2612644/1‑148 (MQ=255)
                                                                       caGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGcc                                                                            >  1:133717/1‑148 (MQ=255)
                                                                         ggCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCAt                                                                          <  1:631167/148‑1 (MQ=255)
                                                                         ggCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCAt                                                                          <  1:631179/148‑1 (MQ=255)
                                                                               cTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAg                                                                    <  2:2202463/148‑1 (MQ=255)
                                                                               cTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATAGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAg                                                                    <  1:673847/148‑1 (MQ=255)
                                                                                        gCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGt                                                           >  2:1194138/1‑148 (MQ=255)
                                                                                        gCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGt                                                           <  1:1507808/148‑1 (MQ=255)
                                                                                          aTCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCg                                                         <  2:2306057/148‑1 (MQ=255)
                                                                                          aTCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCg                                                         <  2:3100183/148‑1 (MQ=255)
                                                                                           tCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCgg                                                        <  2:771917/148‑1 (MQ=255)
                                                                                              cgcgACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGaa                                                     >  1:1813601/1‑148 (MQ=255)
                                                                                                 gACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGa                                                  <  1:3535987/148‑1 (MQ=255)
                                                                                                 gACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGa                                                  <  1:3535586/148‑1 (MQ=255)
                                                                                                     acaACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACg                                              <  1:1009102/148‑1 (MQ=255)
                                                                                                       aaCCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCa                                            <  1:1744385/148‑1 (MQ=255)
                                                                                                       aaCCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCa                                            >  1:2123759/1‑148 (MQ=255)
                                                                                                            ttCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCa                                       <  1:2986736/148‑1 (MQ=255)
                                                                                                            ttCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCa                                       <  1:2986651/148‑1 (MQ=255)
                                                                                                            ttCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCa                                       <  2:1942834/148‑1 (MQ=255)
                                                                                                                 gCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAg                                  >  2:3298823/1‑148 (MQ=255)
                                                                                                                 gCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAg                                  >  1:2393944/1‑148 (MQ=255)
                                                                                                                  cAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGt                                 >  1:3075839/1‑148 (MQ=255)
                                                                                                                  cAGCCTGTTACACGCATCACGATATGCTCATAGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGt                                 <  1:2577125/148‑1 (MQ=255)
                                                                                                                    gCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGttt                               >  2:2556161/1‑148 (MQ=255)
                                                                                                                       tGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCa                            >  1:2725635/1‑148 (MQ=255)
                                                                                                                       tGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCa                            >  2:1314046/1‑148 (MQ=255)
                                                                                                                        gTTACACGCATCACGATATGCTCATAGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCAc                           >  1:1381927/1‑148 (MQ=255)
                                                                                                                        gTTACACGCATCACGATATGCTCATAGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCAc                           >  1:1382208/1‑148 (MQ=255)
                                                                                                                          tACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGc                         >  2:1623308/1‑148 (MQ=255)
                                                                                                                            cacGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTg                       >  2:2974023/1‑148 (MQ=255)
                                                                                                                                 aTCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTgcggcg                  <  2:2481655/148‑1 (MQ=255)
                                                                                                                                 aTCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTgcggcg                  <  2:2481456/148‑1 (MQ=255)
                                                                                                                                   cACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTg                <  1:75350/148‑1 (MQ=255)
                                                                                                                                   cACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTg                <  1:3342251/148‑1 (MQ=255)
                                                                                                                                       atatGCTCATAGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACgg            >  2:1712330/1‑148 (MQ=255)
                                                                                                                                         atGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGca          >  2:3098303/1‑148 (MQ=255)
                                                                                                                                         atGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGca          <  2:2935204/148‑1 (MQ=255)
                                                                                                                                         atGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGca          <  2:2934189/148‑1 (MQ=255)
                                                                                                                                         atGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGca          <  2:3456866/148‑1 (MQ=255)
                                                                                                                                          tGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGcat         <  2:2509445/148‑1 (MQ=255)
                                                                                                                                          tGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGcat         <  2:1100554/148‑1 (MQ=255)
                                                                                                                                           gCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGcatt        >  2:1640070/1‑148 (MQ=255)
                                                                                                                                           gCTCATAGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGTAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGcatt        <  1:70328/148‑1 (MQ=255)
                                                                                                                                           gCTCATAGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGcatt        <  1:3246/148‑1 (MQ=255)
                                                                                                                                            cTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGTCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGcattc       <  2:2291479/148‑1 (MQ=255)
                                                                                                                                              cATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGcattcat     >  1:1686759/1‑148 (MQ=255)
                                                                                                                                               aTCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGcattcatt    >  1:2284698/1‑148 (MQ=255)
                                                                                                                                                 cGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTaa  >  1:1152889/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
ATCCGTGGGATACGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAA  >  CP009273/2881909‑2882201

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: