Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,045,302 G→T 100% intergenic (+86/‑103) yihI → / → hemN activator of Der GTPase/coproporphyrinogen III oxidase, SAM and NAD(P)H dependent, oxygen‑independent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,045,3020GT98.8% 1062.5 / 1.6 326intergenic (+86/‑103)yihI/hemNactivator of Der GTPase/coproporphyrinogen III oxidase, SAM and NAD(P)H dependent, oxygen‑independent
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base T (157/165);  total (159/167)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 333 total aligned reads displayed.

ACGAAGAAGAGGAAGAAGACGAGAAGCAAGAAGACATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCGGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATCCAGAA  >  CP009273/4045155‑4045448
                                                                                                                                                   |                                                                                                                                                  
acgaAGAAGAGGAAGAAGACGAGAAGCAAGAAGACATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGt                                                                                                                                                  <  1:3058078/150‑1 (MQ=255)
 cgaAGAAGGGGAAGAAGACGAGAAGCAAGAAGACATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGtt                                                                                                                                                 >  1:904752/1‑150 (MQ=255)
 cgaAGAAGAGGAAGAAGACGAGAAGCAAGAAGACATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGtt                                                                                                                                                 >  2:3619273/1‑150 (MQ=255)
  gaagaaGAGGAAGAAGACGAGAAGCAAGAAGACATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTg                                                                                                                                                >  2:2902451/1‑150 (MQ=255)
  gaagaaGAGGAAGAAGACGAGAAGCAAGAAGACATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTg                                                                                                                                                <  1:2908185/150‑1 (MQ=255)
  gaagaaGAGGAAGAAGACGAGAAGCAAGAAGACATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTg                                                                                                                                                <  1:2368200/150‑1 (MQ=255)
   aagaagAGGAAGAAGACGAGAAGCAAGAAGACATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCtgttgt                                                                                                                                               >  2:2130240/1‑150 (MQ=255)
      aagaGGAAGAAGACGAGAAGCAAGAAGACATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTcgc                                                                                                                                            >  1:3125183/1‑150 (MQ=255)
       agagGAAGAAGACGAGAAGCAAGAAGACATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTcgcg                                                                                                                                           <  1:2283923/150‑1 (MQ=255)
       agagGAAGAAGACGAGAAGCAAGAAGACATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTACGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTcgcg                                                                                                                                           <  1:2855377/150‑1 (MQ=255)
           gaagaaGACGAGAAGCAAGAAGACATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGa                                                                                                                                       >  1:5150000/1‑150 (MQ=255)
             agaagaCGAGAAGCAAGAAGACATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACa                                                                                                                                     >  1:3437699/1‑150 (MQ=255)
                agaCGAGAAGCAAGAAGACATGATGCGTCTGCTGCGGGTCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACaaaa                                                                                                                                  >  2:2711118/1‑150 (MQ=255)
                agaCGAGAAGCAAGAAGACATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACaaaa                                                                                                                                  >  2:2096017/1‑150 (MQ=255)
                agaCGAGAAGCAAGAAGACATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTAGCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACaaaa                                                                                                                                  >  2:5131949/1‑150 (MQ=255)
                     agaAGCAAGAAGACATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGc                                                                                                                             >  2:5011763/1‑150 (MQ=255)
                       aaGCAAGAAGACATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTg                                                                                                                           >  1:1540896/1‑150 (MQ=255)
                         gCAAGAAGACATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTgcg                                                                                                                         >  1:250403/1‑150 (MQ=255)
                           aagaagACATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCa                                                                                                                       <  2:2566358/150‑1 (MQ=255)
                             gaagaCATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAAc                                                                                                                     >  2:5147877/1‑150 (MQ=255)
                              aagaCATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAAcc                                                                                                                    <  1:26440/150‑1 (MQ=255)
                              aagaCATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAAcc                                                                                                                    >  1:4910492/1‑150 (MQ=255)
                               agaCATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCg                                                                                                                   >  1:4436240/1‑150 (MQ=255)
                                 aCATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGc                                                                                                                 >  2:370476/1‑150 (MQ=255)
                                  cATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCt                                                                                                                <  1:3669328/150‑1 (MQ=255)
                                   atgatgCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCtg                                                                                                               <  1:2129359/150‑1 (MQ=255)
                                          gTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGAccc                                                                                                        >  2:394632/1‑150 (MQ=255)
                                            ctgctgCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGa                                                                                                      >  1:4581072/1‑150 (MQ=255)
                                            ctgctgCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGa                                                                                                      <  1:3955923/150‑1 (MQ=255)
                                                tgcGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACg                                                                                                  <  1:2683183/150‑1 (MQ=255)
                                                   ggggCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGc                                                                                               >  2:174604/1‑150 (MQ=255)
                                                    gggCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCa                                                                                              <  2:2390074/150‑1 (MQ=255)
                                                      gCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATCTGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATc                                                                                            >  2:252707/1‑150 (MQ=255)
                                                       cAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAAGCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATcc                                                                                           <  2:111674/150‑1 (MQ=255)
                                                       cAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATcc                                                                                           <  1:4940737/150‑1 (MQ=255)
                                                       cAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATcc                                                                                           <  1:4469831/150‑1 (MQ=255)
                                                        aaCTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATccg                                                                                          <  1:1657533/150‑1 (MQ=255)
                                                           tAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATccgccg                                                                                       >  1:107287/1‑150 (MQ=255)
                                                             aCGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGt                                                                                     <  2:3259909/150‑1 (MQ=255)
                                                               ggATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCGTCCGCCGGTAc                                                                                   >  1:1165167/1‑150 (MQ=255)
                                                               ggATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTAc                                                                                   <  2:1498016/150‑1 (MQ=255)
                                                                 aTTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGc                                                                                 <  1:4262244/150‑1 (MQ=255)
                                                                  ttGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGcc                                                                                <  2:3692267/150‑1 (MQ=255)
                                                                          gTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAgccgc                                                                        >  1:4680620/1‑150 (MQ=255)
                                                                            gggCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCa                                                                      >  2:3834281/1‑150 (MQ=255)
                                                                             ggCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCag                                                                     <  2:4135177/150‑1 (MQ=255)
                                                                              gCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCaga                                                                    >  1:3489069/1‑150 (MQ=255)
                                                                              gCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCaga                                                                    >  1:1016052/1‑150 (MQ=255)
                                                                                    ccGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGcc                                                              <  1:1921189/150‑1 (MQ=255)
                                                                                       ttCTGCTTATAACCCTTCAGGTTATATGTTATCTGGTGTGTTTATTCTTTAAACTACAGCTGTTGTCGCGGAGACAAAAGTGGCTGCGCAACCGGATGATGACCCGAAAAGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATc                                                           <  1:784336/150‑1 (MQ=255)
                                                                                          tgctgATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGa                                                        <  2:4198949/150‑1 (MQ=255)
                                                                                          tgctgATAACCCTTCCGGTTATAGGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGa                                                        >  2:3056712/1‑150 (MQ=255)
                                                                                           gctgATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGaa                                                       <  1:412495/150‑1 (MQ=255)
                                                                                            ctgATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAg                                                      >  1:2366102/1‑150 (MQ=255)
                                                                                            ctgATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAg                                                      >  1:4350738/1‑150 (MQ=255)
                                                                                            ctgATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAg                                                      >  2:1922165/1‑150 (MQ=255)
                                                                                            ctgATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCGGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAg                                                      <  2:2069233/150‑1 (MQ=255)
                                                                                              gATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGa                                                    <  1:5011763/150‑1 (MQ=255)
                                                                                               aTAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAg                                                   >  1:58782/1‑150 (MQ=255)
                                                                                               aTAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAg                                                   >  2:1008213/1‑150 (MQ=255)
                                                                                                      ttCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCAt                                            <  1:663687/150‑1 (MQ=255)
                                                                                                      ttCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCAt                                            <  1:2086938/150‑1 (MQ=255)
                                                                                                      ttCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCAt                                            >  1:1438822/1‑150 (MQ=255)
                                                                                                        ccGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGt                                          >  1:2034488/1‑150 (MQ=255)
                                                                                                        ccGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGt                                          >  2:46533/1‑150 (MQ=255)
                                                                                                        ccGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCGGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGt                                          >  1:2997013/1‑150 (MQ=255)
                                                                                                              atatGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGGTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGACACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTAc                                    <  2:4805124/150‑1 (MQ=255)
                                                                                                                atGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAg                                  >  2:2593872/1‑150 (MQ=255)
                                                                                                                 tGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGc                                 >  2:3439582/1‑150 (MQ=255)
                                                                                                                 tGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGc                                 >  1:802906/1‑150 (MQ=255)
                                                                                                                   ttATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAATAGTGGCTGCGCAACCGGCTGATGACCCGAAACTTGCATCCGCCGGTACGCCGTAGCCGCCAGAGAAGCCATCGGAAGGAGTGAGAATGTCTGTACatcta                               <  2:626939/150‑5 (MQ=255)
                                                                                                                   ttATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCaa                               >  1:4784501/1‑150 (MQ=255)
                                                                                                                    tATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCaaa                              >  2:2461545/1‑150 (MQ=255)
                                                                                                                     aTCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAAt                             >  1:1767163/1‑150 (MQ=255)
                                                                                                                     aTCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAAt                             <  1:2802561/150‑1 (MQ=255)
                                                                                                                         ggTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGAc                         <  1:3258307/150‑1 (MQ=255)
                                                                                                                          gtgtGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACt                        <  2:2961987/150‑1 (MQ=255)
                                                                                                                           tgtgGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTg                       >  2:2600167/1‑150 (MQ=255)
                                                                                                                            gtgGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTgg                      >  1:4028558/1‑150 (MQ=255)
                                                                                                                            gtgGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTgg                      >  2:2040566/1‑150 (MQ=255)
                                                                                                                             tgGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTggg                     >  1:296387/1‑150 (MQ=255)
                                                                                                                             tgGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTggg                     <  1:2094798/150‑1 (MQ=255)
                                                                                                                                ttATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATc                  <  2:725464/150‑1 (MQ=255)
                                                                                                                                  aTTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTg                <  2:4918628/150‑1 (MQ=255)
                                                                                                                                     cGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGcc             <  2:4358200/150‑1 (MQ=255)
                                                                                                                                      gTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGccc            >  2:1287021/1‑150 (MQ=255)
                                                                                                                                        tAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTg          <  1:5125542/150‑1 (MQ=255)
                                                                                                                                        tAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTg          <  2:680688/150‑1 (MQ=255)
                                                                                                                                        tAAACGACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTg          <  2:2294038/150‑1 (MQ=255)
                                                                                                                                         aaaCTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGa         >  1:1192716/1‑150 (MQ=255)
                                                                                                                                         aaaCTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGa         >  2:1634476/1‑150 (MQ=255)
                                                                                                                                          aaCTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTTGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGAt        <  2:4750246/150‑1 (MQ=255)
                                                                                                                                          aaCTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAACTCGACTGGGATCTGGCCCTGAt        >  2:1818044/1‑150 (MQ=255)
                                                                                                                                          aaCTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGAt        >  2:4579100/1‑150 (MQ=255)
                                                                                                                                           aCTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATc       <  2:3807869/150‑1 (MQ=255)
                                                                                                                                           aCTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATc       <  1:3117349/150‑1 (MQ=255)
                                                                                                                                            cTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATcc      <  1:3573795/150‑1 (MQ=255)
                                                                                                                                             tACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATCCa     <  2:3388218/150‑1 (MQ=255)
                                                                                                                                             tACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATCCa     <  2:812427/150‑1 (MQ=255)
                                                                                                                                              aCAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATCCAg    >  1:333455/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
ACGAAGAAGAGGAAGAAGACGAGAAGCAAGAAGACATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCGGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATCCAGAA  >  CP009273/4045155‑4045448

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: