Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,533,556 G→T 17.2% H1098N (CAC→AAC)  narZ ← nitrate reductase 2 (NRZ), alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,533,5560GT17.2% 963.4 / 128.2 314H1098N (CAC→AAC) narZnitrate reductase 2 (NRZ), alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (134/126);  new base T (25/29);  total (159/155)
Fisher's exact test for biased strand distribution p-value = 5.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 319 total aligned reads displayed.

GCGGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGCCACCAGGCTTTCACCAAAAGCACGCATCCACGGAT  >  CP009273/1533407‑1533705
                                                                                                                                                     |                                                                                                                                                     
   ggCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTTAAt                                                                                                                                                    <  1:1680078/150‑1 (MQ=255)
   ggCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTTAAt                                                                                                                                                    <  1:3412183/150‑1 (MQ=255)
    gCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTTAATg                                                                                                                                                   <  1:1291678/150‑1 (MQ=255)
      aGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACGGTAGGTTGAGTTAATGcc                                                                                                                                                 >  1:189717/1‑150 (MQ=255)
        tCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGcccc                                                                                                                                               <  1:929203/150‑1 (MQ=255)
          aGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCAt                                                                                                                                             >  2:2844198/1‑150 (MQ=255)
           gcgcgcCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTTAATGCCCCAtt                                                                                                                                            <  2:4250422/150‑1 (MQ=255)
             gcgcCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTTAATGCCCCATTTc                                                                                                                                          <  2:3642427/150‑1 (MQ=255)
               gcCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTTAATGCCCCATTTCTg                                                                                                                                        <  1:4491717/150‑1 (MQ=255)
               gcCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTg                                                                                                                                        <  2:4143291/150‑1 (MQ=255)
                ccGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCAACTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTgg                                                                                                                                       >  1:2125805/1‑150 (MQ=255)
                   ttGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTgc                                                                                                                                    <  2:786620/150‑1 (MQ=255)
                   ttGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTgc                                                                                                                                    >  1:259401/1‑150 (MQ=255)
                    tGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTgcg                                                                                                                                   >  2:1683986/1‑150 (MQ=255)
                       cGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTgcggcg                                                                                                                                <  1:1972904/150‑1 (MQ=255)
                       cGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTgcggcg                                                                                                                                >  1:2485421/1‑150 (MQ=255)
                       cGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTgcggcg                                                                                                                                >  1:327345/1‑150 (MQ=255)
                        gTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGt                                                                                                                               >  2:2661612/1‑150 (MQ=255)
                         ttGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTc                                                                                                                              <  1:810907/150‑1 (MQ=255)
                         ttGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTc                                                                                                                              >  1:2831534/1‑150 (MQ=255)
                          tGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCa                                                                                                                             >  2:1616085/1‑150 (MQ=255)
                           gAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAg                                                                                                                            >  2:3684339/1‑150 (MQ=255)
                            aaTACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAgg                                                                                                                           <  2:1111721/150‑1 (MQ=255)
                             aTACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTAGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTATGGTGCGGCGTCAGGa                                                                                                                          >  1:4350731/1‑150 (MQ=255)
                              tACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGaa                                                                                                                         >  2:2609647/1‑150 (MQ=255)
                               aCTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAg                                                                                                                        <  1:769363/150‑1 (MQ=255)
                                 ttCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGtt                                                                                                                      <  2:3140296/150‑1 (MQ=255)
                                   ccACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTTAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTaa                                                                                                                    <  1:3921678/150‑1 (MQ=255)
                                        cAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCt                                                                                                               <  1:2987276/150‑1 (MQ=255)
                                        cAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCt                                                                                                               >  1:3100518/1‑150 (MQ=255)
                                        cAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCt                                                                                                               <  1:1154482/150‑1 (MQ=255)
                                            cGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCttttt                                                                                                           >  1:2555380/1‑150 (MQ=255)
                                            cGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCttttt                                                                                                           >  2:445153/1‑150 (MQ=255)
                                             gTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTc                                                                                                          <  2:2354434/150‑1 (MQ=255)
                                               tATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCg                                                                                                        <  2:3932622/150‑1 (MQ=255)
                                               tATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCg                                                                                                        <  1:1083661/150‑1 (MQ=255)
                                                 tCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCggg                                                                                                      <  2:938872/150‑1 (MQ=255)
                                                  cGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTTAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGa                                                                                                     >  2:187359/1‑150 (MQ=255)
                                                  cGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGa                                                                                                     >  2:1464014/1‑150 (MQ=255)
                                                    aCAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTTAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAg                                                                                                   <  1:487624/150‑1 (MQ=255)
                                                    aCAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAg                                                                                                   <  1:3806581/150‑1 (MQ=255)
                                                    aCAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAg                                                                                                   >  1:829225/1‑150 (MQ=255)
                                                        tGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGt                                                                                               >  1:1747515/1‑150 (MQ=255)
                                                         ggTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGtt                                                                                              >  2:3414506/1‑150 (MQ=255)
                                                          gTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTTAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGttt                                                                                             >  2:1119237/1‑150 (MQ=255)
                                                          gTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGttt                                                                                             >  1:70212/1‑150 (MQ=255)
                                                          gTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGttt                                                                                             <  1:3297615/150‑1 (MQ=255)
                                                              gTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGGCg                                                                                         >  2:1321056/1‑150 (MQ=255)
                                                                  tcGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGGCGGTAt                                                                                     >  1:2010961/1‑150 (MQ=255)
                                                                   cGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGGCGGTATc                                                                                    <  2:2928825/150‑1 (MQ=255)
                                                                     ggCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGGCGGTATCTg                                                                                  <  2:3391288/150‑1 (MQ=255)
                                                                      gCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGGCGGTATCTgg                                                                                 >  2:1339596/1‑150 (MQ=255)
                                                                        aTCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGGCGGTATCTGgcg                                                                               <  2:3079646/150‑1 (MQ=255)
                                                                             tttCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGGCGGTATCTGGCGCAtct                                                                          >  2:3773462/1‑150 (MQ=255)
                                                                                cGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCAc                                                                       <  2:2796199/150‑1 (MQ=255)
                                                                                 gCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTTAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACt                                                                      <  2:1263517/150‑1 (MQ=255)
                                                                                  cTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTg                                                                     <  2:668177/150‑1 (MQ=255)
                                                                                   tGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGa                                                                    >  1:4582337/1‑150 (MQ=255)
                                                                                   tGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGa                                                                    >  1:2262029/1‑150 (MQ=255)
                                                                                       ccAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGCt                                                                <  1:4323964/150‑1 (MQ=255)
                                                                                        cAGACAATCGGTCCACCGCGAGAGAGCGTCNGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGCTa                                                               >  1:4282371/1‑150 (MQ=255)
                                                                                        cAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGCTa                                                               >  1:4116559/1‑150 (MQ=255)
                                                                                         aGACAATCGGTCCACCGCTAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGCTAc                                                              <  2:2019135/150‑1 (MQ=255)
                                                                                            cAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGCTACggg                                                           >  2:494502/1‑150 (MQ=255)
                                                                                            cAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGCTACggg                                                           <  1:1337521/150‑1 (MQ=255)
                                                                                             aaTCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGCTACGGgt                                                          <  2:3518497/150‑1 (MQ=255)
                                                                                             aaTCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACAGACGCTACGGgt                                                          >  2:1611280/1‑150 (MQ=255)
                                                                                               tCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGCTACGGgtgt                                                        >  1:1533572/1‑150 (MQ=255)
                                                                                                cGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTTAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTc                                                       >  2:3022013/1‑150 (MQ=255)
                                                                                                   tCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGAt                                                    >  1:2892146/1‑150 (MQ=255)
                                                                                                   tCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGAt                                                    <  1:2633757/150‑1 (MQ=255)
                                                                                                      aCCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAgg                                                 >  2:3356400/1‑150 (MQ=255)
                                                                                                      aCCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAgg                                                 <  1:1451520/150‑1 (MQ=255)
                                                                                                        cgcgAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAggcg                                               >  1:1222202/1‑150 (MQ=255)
                                                                                                               gagCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTTAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATa                                        >  1:542301/1‑150 (MQ=255)
                                                                                                               gagCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATa                                        <  2:874694/150‑1 (MQ=255)
                                                                                                                   gTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATGTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGTGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGcc                                    <  2:1900512/150‑1 (MQ=255)
                                                                                                                    tCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGcca                                   <  1:4013644/150‑1 (MQ=255)
                                                                                                                     cAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGccac                                  >  2:2214410/1‑150 (MQ=255)
                                                                                                                      aGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGccacc                                 >  2:1943343/1‑150 (MQ=255)
                                                                                                                      aGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGccacc                                 <  2:4116686/150‑1 (MQ=255)
                                                                                                                      aGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGccacc                                 >  2:3796282/1‑150 (MQ=255)
                                                                                                                       gCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGccacca                                <  2:4317021/150‑1 (MQ=255)
                                                                                                                            aGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGCCACCAGGCtt                           <  2:2517202/150‑1 (MQ=255)
                                                                                                                              cAGGTTTTCACTGTAGGTTGAGTTAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGCCACCAGGCTTTc                         <  2:3202091/150‑1 (MQ=255)
                                                                                                                              cAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGCCACCAGGCTTTc                         >  2:4387536/1‑150 (MQ=255)
                                                                                                                              cAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGCCACCAGGCTTTc                         <  2:1136023/150‑1 (MQ=255)
                                                                                                                              cAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGCCACCAGGCTTTc                         >  2:2898000/1‑150 (MQ=255)
                                                                                                                               aGGTTTTCACTGTAGGTTGAGTTAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGCCACCAGGCTTTCa                        >  1:1632242/1‑150 (MQ=255)
                                                                                                                               aGGTTTTCACTGTAGGTTGAGTTAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGCCACCAGGCTTTCa                        <  1:990018/150‑1 (MQ=255)
                                                                                                                                      cACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGCCACCAGGCTTTCACCAAAAg                 <  2:3548889/150‑1 (MQ=255)
                                                                                                                                        cTGTAGGTTGAGTTAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGCCACCAGGCTTTCACCAAAAGCa               >  1:3141188/1‑150 (MQ=255)
                                                                                                                                          gTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGCCACCAGGCTTTCACCAAAAGCACg             <  1:3547975/150‑1 (MQ=255)
                                                                                                                                                 gAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGCCACCAGGCTTTCACCAAAAGCACGCATCCAc      <  1:2359912/150‑1 (MQ=255)
                                                                                                                                                 gAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGCCACCAGGCTTTCACCAAAAGCACGCATCCAc      <  1:1509766/150‑1 (MQ=255)
                                                                                                                                                 gAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGCCACCAGGCTTTCACCAAAAGCACGCATCCAc      >  2:331569/1‑150 (MQ=255)
                                                                                                                                                 gAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGCCACCAGGCTTTCACCAAAAGCACGCATCCAc      <  2:3000526/150‑1 (MQ=255)
                                                                                                                                                 gAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGCCACCAGGCTTTCACCAAAAGCACGCATCCAc      >  1:3358980/1‑150 (MQ=255)
                                                                                                                                                     gAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGCCACCAGGCTTTCACCAAAAGCACGCATCCACGgat  >  2:1815557/1‑150 (MQ=255)
                                                                                                                                                     gAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGCCACCAGGCTTTCACCAAAAGCACGCATCCACGgat  <  1:100819/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCGGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGCCACCAGGCTTTCACCAAAAGCACGCATCCACGGAT  >  CP009273/1533407‑1533705

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: