Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,017,285 G→T 100% D871Y (GAC→TAC)  xdhD → putative hypoxanthine oxidase, molybdopterin‑binding/Fe‑S binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,017,2850GT96.7% 895.0 / 16.1 272D871Y (GAC→TAC) xdhDputative hypoxanthine oxidase, molybdopterin‑binding/Fe‑S binding
Reads supporting (aligned to +/- strand):  ref base G (5/4);  new base T (146/117);  total (151/121)
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.94e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 274 total aligned reads displayed.

AACACGCGTACGGGTGAAATCCGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGT  >  CP009273/3017138‑3017431
                                                                                                                                                   |                                                                                                                                                  
aaCACGCGTACGGGTGAAATCCGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTAc                                                                                                                                                  >  1:1078501/1‑150 (MQ=255)
 acacGCGTACGGGTGAAATCCGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACg                                                                                                                                                 <  2:2987550/150‑1 (MQ=255)
   acGCGTACGGGTGAAATCCGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGcc                                                                                                                                               <  1:3189874/150‑1 (MQ=255)
     gcgTACGGGTGAAATCCGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGa                                                                                                                                             >  2:418258/1‑150 (MQ=255)
     gcgTACGGGTGAAATCCGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGa                                                                                                                                             <  1:3878067/150‑1 (MQ=255)
     gcgTACGGGTGAAATCCGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGa                                                                                                                                             <  1:557899/150‑1 (MQ=255)
     gcgTACGGGTGAAATCCGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGa                                                                                                                                             <  2:3400344/150‑1 (MQ=255)
      cgTACGGGTGAAATCCGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGaa                                                                                                                                            <  2:2911670/150‑1 (MQ=255)
           gggTGAAATCCGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCa                                                                                                                                       >  1:319377/1‑150 (MQ=255)
             gTGAAATCCGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCAcc                                                                                                                                     >  1:3759749/1‑150 (MQ=255)
              tGAAATCCGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCAccc                                                                                                                                    <  2:4436884/150‑1 (MQ=255)
              tGAAATCCGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCAccc                                                                                                                                    <  1:2059153/150‑1 (MQ=255)
                 aaTCCGCCTTGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGtt                                                                                                                                 <  2:417400/150‑1 (MQ=255)
                       ccTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAAcgcg                                                                                                                           >  2:3723966/1‑150 (MQ=255)
                        cTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGt                                                                                                                          >  2:2158589/1‑150 (MQ=255)
                           gATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGAt                                                                                                                       >  2:3990576/1‑150 (MQ=255)
                             tAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGAttt                                                                                                                     >  2:3715417/1‑150 (MQ=255)
                              aaaTTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTa                                                                                                                    >  2:3431410/1‑150 (MQ=255)
                                 ttCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGc                                                                                                                 <  1:1027447/150‑1 (MQ=255)
                                   cTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGACACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAg                                                                                                               >  2:762699/1‑150 (MQ=255)
                                    tACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGt                                                                                                              >  2:3513923/1‑150 (MQ=255)
                                     aCGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGtt                                                                                                             >  1:291993/1‑150 (MQ=255)
                                     aCGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGtt                                                                                                             >  1:1539738/1‑150 (MQ=255)
                                      cGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTa                                                                                                            >  2:2021099/1‑150 (MQ=255)
                                      cGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTa                                                                                                            <  1:2673792/150‑1 (MQ=255)
                                       gCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTAc                                                                                                           <  2:1354154/150‑1 (MQ=255)
                                        ccTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACg                                                                                                          >  2:1728264/1‑150 (MQ=255)
                                        ccTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACg                                                                                                          <  2:3436135/150‑1 (MQ=255)
                                        ccTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACg                                                                                                          >  2:368645/1‑150 (MQ=255)
                                          ttGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGgc                                                                                                        >  1:4292573/1‑150 (MQ=255)
                                             cTGGACTGCGGTACACCGGTAAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCa                                                                                                     >  2:2120271/1‑150 (MQ=255)
                                              tGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCAc                                                                                                    >  2:871265/1‑150 (MQ=255)
                                               ggACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCAcc                                                                                                   >  1:4035259/1‑150 (MQ=255)
                                               ggACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCAcc                                                                                                   >  1:1133003/1‑150 (MQ=255)
                                                gACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCg                                                                                                  >  2:1172262/1‑150 (MQ=255)
                                                gACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCg                                                                                                  <  1:8175/150‑1 (MQ=255)
                                                gACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTCGGACAACTCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGACGAGATCATTTATTACGCCGAAGGTCACCCGTTAACCCGTGATTTACGCAGTTACGGCGCACCg                                                                                                  <  2:1934343/150‑1 (MQ=255)
                                                gACAGCGGTACCCCGGTCACACCAGAGATAGCGTTGGGCCAAAGCTACGCTGCCACCCTGCGAGCTATCGGCCACAGTATGAGCGAACAGATCATTTATGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCg                                                                                                  <  2:1811203/150‑1 (MQ=255)
                                                 aCTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGa                                                                                                 >  1:846348/1‑150 (MQ=255)
                                                     cGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAAt                                                                                             >  2:3274307/1‑150 (MQ=255)
                                                      ggTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAAtt                                                                                            >  1:3821264/1‑150 (MQ=255)
                                                      ggTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAAtt                                                                                            >  1:3821374/1‑150 (MQ=255)
                                                      ggTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAAtt                                                                                            >  1:525496/1‑150 (MQ=255)
                                                       gTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTg                                                                                           <  1:2002485/150‑1 (MQ=255)
                                                         acacCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGt                                                                                         >  2:1381051/1‑150 (MQ=255)
                                                         acacCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGt                                                                                         <  1:4212370/150‑1 (MQ=255)
                                                         acacCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGt                                                                                         >  2:2288709/1‑150 (MQ=255)
                                                                       agagTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGa                                                                           >  2:536579/1‑150 (MQ=255)
                                                                          gTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGAttt                                                                        >  1:606148/1‑150 (MQ=255)
                                                                           ttAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTc                                                                       >  2:2699240/1‑150 (MQ=255)
                                                                           ttAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTc                                                                       >  2:4500632/1‑150 (MQ=255)
                                                                            tAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTcc                                                                      <  1:3472166/150‑1 (MQ=255)
                                                                             aGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCcg                                                                     <  1:3114854/150‑1 (MQ=255)
                                                                               cGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCcgcg                                                                   <  2:2922398/150‑1 (MQ=255)
                                                                                  tGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCtg                                                                <  2:4420748/150‑1 (MQ=255)
                                                                                    ggACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCtgtg                                                              >  1:1272088/1‑150 (MQ=255)
                                                                                       cAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTg                                                           >  1:273421/1‑150 (MQ=255)
                                                                                        aaaTCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTgg                                                          >  2:2625224/1‑150 (MQ=255)
                                                                                          aTCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTg                                                        >  2:962151/1‑150 (MQ=255)
                                                                                          aTCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTg                                                        <  1:1746314/150‑1 (MQ=255)
                                                                                          aTCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTg                                                        <  1:624461/150‑1 (MQ=255)
                                                                                              aCGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGa                                                    <  1:1607871/150‑1 (MQ=255)
                                                                                              aCGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGa                                                    <  2:2011086/150‑1 (MQ=255)
                                                                                               cGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAATTTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAg                                                   >  1:1739987/1‑150 (MQ=255)
                                                                                                   gCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGcgac                                               <  1:4020851/150‑1 (MQ=255)
                                                                                                    ccACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGcgacg                                              <  1:2414777/150‑1 (MQ=255)
                                                                                                        ccTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATaa                                          >  2:2907069/1‑150 (MQ=255)
                                                                                                         cTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATaaa                                         >  1:396890/1‑150 (MQ=255)
                                                                                                         cTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATaaa                                         >  1:1630064/1‑150 (MQ=255)
                                                                                                          tGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAg                                        >  2:1608321/1‑150 (MQ=255)
                                                                                                             gTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCg                                     <  1:368645/150‑1 (MQ=255)
                                                                                                                 tATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGccc                                 >  1:3166291/1‑150 (MQ=255)
                                                                                                                  aTCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCg                                <  1:3723946/150‑1 (MQ=255)
                                                                                                                    cGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGtt                              >  1:2288118/1‑150 (MQ=255)
                                                                                                                     ggCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTc                             >  1:177257/1‑150 (MQ=255)
                                                                                                                      gCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCg                            <  2:2224334/150‑1 (MQ=255)
                                                                                                                      gCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCg                            >  2:314169/1‑150 (MQ=255)
                                                                                                                      gCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCg                            >  2:313878/1‑150 (MQ=255)
                                                                                                                      gCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCg                            >  2:2037937/1‑150 (MQ=255)
                                                                                                                        cacaGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCggg                          <  2:2820631/150‑1 (MQ=255)
                                                                                                                        cacaGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCggg                          <  2:791975/150‑1 (MQ=255)
                                                                                                                             tATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGaa                     >  1:1948583/1‑150 (MQ=255)
                                                                                                                              aTGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGaca                    >  2:2386628/1‑148 (MQ=255)
                                                                                                                               tGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAAt                   >  2:447938/1‑150 (MQ=255)
                                                                                                                                gAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATc                  <  2:3801693/150‑1 (MQ=255)
                                                                                                                                 aGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCg                 >  2:3312482/1‑150 (MQ=255)
                                                                                                                                 aGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCg                 >  2:3245076/1‑150 (MQ=255)
                                                                                                                                  gCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGa                <  1:1235867/150‑1 (MQ=255)
                                                                                                                                   cGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGAt               >  1:140589/1‑150 (MQ=255)
                                                                                                                                   cGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGAt               <  1:3291859/150‑1 (MQ=255)
                                                                                                                                    gAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGAtc              >  2:259599/1‑150 (MQ=255)
                                                                                                                                    gAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGAtc              >  1:3090056/1‑150 (MQ=255)
                                                                                                                                    gAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGAtc              <  1:3605185/150‑1 (MQ=255)
                                                                                                                                     aaGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATTGAtct             <  2:765637/150‑1 (MQ=255)
                                                                                                                                        agaTCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCgg          >  2:2255815/1‑150 (MQ=255)
                                                                                                                                             aTTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATc     >  1:2473624/1‑150 (MQ=255)
                                                                                                                                              tttATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCg    <  1:1239369/150‑1 (MQ=255)
                                                                                                                                              tttATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCg    <  1:1239411/150‑1 (MQ=255)
                                                                                                                                              tttATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCg    <  2:1269942/150‑1 (MQ=255)
                                                                                                                                                tATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGgt  >  2:4254030/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AACACGCGTACGGGTGAAATCCGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGT  >  CP009273/3017138‑3017431

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: