Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,035,421 G→T 94.0% intergenic (‑62/‑8) mobA ← / → yihD molybdopterin‑guanine dinucleotide synthase/DUF1040 protein YihD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,035,4210GT94.0% 813.8 / 42.1 268intergenic (‑62/‑8)mobA/yihDmolybdopterin‑guanine dinucleotide synthase/DUF1040 protein YihD
Reads supporting (aligned to +/- strand):  ref base G (5/11);  new base T (125/127);  total (130/138)
Fisher's exact test for biased strand distribution p-value = 2.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.72e-01

Only 100 of 276 total aligned reads displayed.

TCAAGCAATCCTTTATCTACGCCGCCCATTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGAGCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGTTTTGACGGCGAACTGGCGGATTTGACGGATG  >  CP009273/4035274‑4035570
                                                                                                                                                   |                                                                                                                                                     
 cAAGCAATCCTTTATCTACGCCGCCCATTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCAt                                                                                                                                                    <  1:3158784/150‑1 (MQ=255)
 cAAGCAATCCTTTATCTACGCCGCCCATTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCAt                                                                                                                                                    >  1:4213458/1‑150 (MQ=255)
 cAAGCAATCCTTTATCTACGCCGCCCATTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCAt                                                                                                                                                    >  1:1380030/1‑150 (MQ=255)
     cAATCCTTTATCTACGCCGCCCATTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACt                                                                                                                                                >  2:3552222/1‑150 (MQ=255)
     cAATCCTTTATCTACGCCGCCCATTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACt                                                                                                                                                >  1:2792158/1‑150 (MQ=255)
       aTCCTTTATCTACGCCGCCCATTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTAt                                                                                                                                              >  1:1957229/1‑150 (MQ=255)
          cTTTATCTACGCCGCCCATTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGaa                                                                                                                                           <  2:3569689/150‑1 (MQ=255)
             tATCTACGCCGCCCATTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATg                                                                                                                                        >  2:3904631/1‑150 (MQ=255)
             tATCTACGCCGCCCATTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATg                                                                                                                                        <  1:2958995/150‑1 (MQ=255)
                 tACGCCGCCCATTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTaaa                                                                                                                                    <  2:4535252/150‑1 (MQ=255)
                  aCGCCGCCCATTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGAGCATCACTATGAAATGTAAAc                                                                                                                                   >  1:889844/1‑150 (MQ=255)
                   cgccgcCCATTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACg                                                                                                                                  <  2:1237824/150‑1 (MQ=255)
                     ccgccCATTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTc                                                                                                                                >  2:921034/1‑150 (MQ=255)
                       gccCATTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCtg                                                                                                                              >  2:2913753/1‑150 (MQ=255)
                       gccCATTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCtg                                                                                                                              >  2:3910923/1‑150 (MQ=255)
                         ccATTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCtgaa                                                                                                                            >  2:1946464/1‑150 (MQ=255)
                           aTTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCtgaatg                                                                                                                          >  2:3659814/1‑150 (MQ=255)
                           aTTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGAGCATCACTATGAAATGTAAACGTCtgaatg                                                                                                                          >  2:303928/1‑150 (MQ=255)
                             tCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCtgaatgaa                                                                                                                        <  1:2565791/150‑1 (MQ=255)
                              cGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAg                                                                                                                       >  2:4136388/1‑150 (MQ=255)
                               gTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGt                                                                                                                      <  2:870689/150‑1 (MQ=255)
                                 cTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTa                                                                                                                    <  2:1862352/150‑1 (MQ=255)
                                 cTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTa                                                                                                                    <  2:1497980/150‑1 (MQ=255)
                                 cTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTa                                                                                                                    <  1:1364332/150‑1 (MQ=255)
                                  tGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTAt                                                                                                                   >  2:4559760/1‑150 (MQ=255)
                                   ggCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTAtt                                                                                                                  >  2:3385460/1‑150 (MQ=255)
                                    gCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTg                                                                                                                 <  1:2129947/150‑1 (MQ=255)
                                      tttACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGaa                                                                                                               >  2:1629319/1‑150 (MQ=255)
                                      tttACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGaa                                                                                                               >  2:1051941/1‑150 (MQ=255)
                                      tttACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGaa                                                                                                               >  2:480974/1‑150 (MQ=255)
                                      tttACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGaa                                                                                                               >  2:2792221/1‑150 (MQ=255)
                                       ttACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAAc                                                                                                              <  2:512196/150‑1 (MQ=255)
                                         aCCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGAGCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAActc                                                                                                            <  2:2046164/150‑1 (MQ=255)
                                          ccGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAActcc                                                                                                           <  1:1798347/150‑1 (MQ=255)
                                                  cAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTccagcca                                                                                                   <  1:4094475/150‑1 (MQ=255)
                                                     cacaACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGTGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCagccagcc                                                                                                <  1:1252248/150‑1 (MQ=255)
                                                     cacaACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCagccagcc                                                                                                <  1:3095137/150‑1 (MQ=255)
                                                       caACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTg                                                                                              >  2:3163473/1‑150 (MQ=255)
                                                        aaCGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTgg                                                                                             >  2:1594851/1‑150 (MQ=255)
                                                        aaCGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTgg                                                                                             >  2:1008077/1‑150 (MQ=255)
                                                        aaCGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGAGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCAATATGAAATGTAAACGTCTGAATGAAGTAATTGAACTCCTCCAGCCAGCCTgg                                                                                             >  2:2101137/1‑150 (MQ=255)
                                                             cTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCaaaa                                                                                        >  1:1005031/1‑150 (MQ=255)
                                                               gTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAag                                                                                      <  2:2137086/150‑1 (MQ=255)
                                                                 tATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAagag                                                                                    <  1:3515829/150‑1 (MQ=255)
                                                                  aTCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGc                                                                                   >  2:3598079/1‑150 (MQ=255)
                                                                     gtcgtcATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAg                                                                                <  1:602710/150‑1 (MQ=255)
                                                                     gtcgtcATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGAGCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAg                                                                                <  2:3798723/150‑1 (MQ=255)
                                                                         tcatcaGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCt                                                                            >  2:3228793/1‑150 (MQ=255)
                                                                         tcatcaGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCt                                                                            >  2:4196978/1‑150 (MQ=255)
                                                                          catcaGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCtt                                                                           >  1:2288457/1‑150 (MQ=255)
                                                                           atcaGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTa                                                                          <  2:2575561/150‑1 (MQ=255)
                                                                            tcaGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTaa                                                                         >  2:17266/1‑150 (MQ=255)
                                                                             caGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAAc                                                                        <  2:62350/150‑1 (MQ=255)
                                                                             caGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAAc                                                                        >  2:1505209/1‑150 (MQ=255)
                                                                             caGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAAc                                                                        <  2:4562440/150‑1 (MQ=255)
                                                                              aGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAAcc                                                                       >  1:1915127/1‑150 (MQ=255)
                                                                              aGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAAcc                                                                       >  1:2472336/1‑150 (MQ=255)
                                                                              aGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAAcc                                                                       >  1:3342223/1‑150 (MQ=255)
                                                                               gATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCAACACTATGAAATGTAAAAGTCTGAATGAAGTTATTGAACTCATCCAGCCAGACTGGCAAAAAGAGCCAGACCTTAACct                                                                      >  2:3716885/1‑150 (MQ=255)
                                                                                aTTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACctg                                                                     >  1:172207/1‑150 (MQ=255)
                                                                                   cACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACctgctg                                                                  <  1:3435303/150‑1 (MQ=255)
                                                                                    aCCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCtgctgc                                                                 <  1:1658864/150‑1 (MQ=255)
                                                                                            cGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTg                                                         >  2:606252/1‑150 (MQ=255)
                                                                                            cGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTg                                                         >  2:1398961/1‑150 (MQ=255)
                                                                                              ccTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCa                                                       >  1:3607404/1‑150 (MQ=255)
                                                                                              ccTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCa                                                       >  2:29406/1‑150 (MQ=255)
                                                                                                 ctTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGaa                                                    >  1:3950466/1‑150 (MQ=255)
                                                                                                 ctTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGaa                                                    >  1:2215814/1‑150 (MQ=255)
                                                                                                  ttttATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGaaa                                                   >  2:2487612/1‑150 (MQ=255)
                                                                                                   tttATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAAc                                                  >  1:3682954/1‑150 (MQ=255)
                                                                                                      atggtgGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTgg                                               <  2:2450291/147‑1 (MQ=255)
                                                                                                      aTTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTgg                                               <  2:1862420/150‑1 (MQ=255)
                                                                                                      aTTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTgg                                               <  2:3087047/150‑1 (MQ=255)
                                                                                                       ttGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGAGCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGc                                              <  1:375352/150‑1 (MQ=255)
                                                                                                          tgGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGaa                                           <  1:2871975/150‑1 (MQ=255)
                                                                                                           gggATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGaaa                                          >  1:4263498/1‑150 (MQ=255)
                                                                                                           gggATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGaaa                                          >  2:1236517/1‑150 (MQ=255)
                                                                                                           gggATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGaaa                                          >  2:3550397/1‑150 (MQ=255)
                                                                                                            ggATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAag                                         <  1:812082/150‑1 (MQ=255)
                                                                                                              aTTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAagag                                       <  1:3365019/150‑1 (MQ=255)
                                                                                                               ttGACCCTGCTAACGTGTCTGTCTCAAGAGTAANGATCATCACNATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGt                                      <  1:1014145/150‑1 (MQ=255)
                                                                                                                    ccTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGt                                 <  2:1870097/150‑1 (MQ=255)
                                                                                                                     cTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGtt                                <  2:3446948/150‑1 (MQ=255)
                                                                                                                       gCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGtttt                              <  2:282126/150‑1 (MQ=255)
                                                                                                                         tAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGTTTTGa                            >  1:752848/1‑150 (MQ=255)
                                                                                                                         tAACGTGTCTGTCTCAAGAGTAAGGAGCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGTTTTGa                            >  1:4181180/1‑150 (MQ=255)
                                                                                                                               gtctgtctCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGTTTTGACGGCGa                      >  1:1911011/1‑150 (MQ=255)
                                                                                                                               gtctgtctCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGTTTTGACGGCGa                      >  1:1205075/1‑150 (MQ=255)
                                                                                                                                 ctgtctCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGTTTTGACGGCGAAc                    >  1:4062582/1‑150 (MQ=255)
                                                                                                                                  tgtctCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGTTTTGACGGCGAACt                   >  1:494603/1‑150 (MQ=255)
                                                                                                                                      tcAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTACTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAATAGTCAGGTTTTGACGGCGAACTGGCg               >  2:2239382/1‑150 (MQ=255)
                                                                                                                                       cAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGTTTTGACGGCGAACTGGCgg              <  2:2194271/150‑1 (MQ=255)
                                                                                                                                          gagTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGTTTTGACGGCGAACTGGCGGAtt           >  1:4395978/1‑150 (MQ=255)
                                                                                                                                          gagTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGTTTTGACGGCGAACTGGCGGAtt           >  1:1469091/1‑150 (MQ=255)
                                                                                                                                           agTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGTTTTGACGGCGAACTGGCGGAttt          <  2:3289207/150‑1 (MQ=255)
                                                                                                                                           agTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGTTTTGACGGCGAACTGGCGGAttt          <  1:4520457/150‑1 (MQ=255)
                                                                                                                                           agTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGTTTTGACGGCGAACTGGCGGAttt          <  2:850592/150‑1 (MQ=255)
                                                                                                                                              aaGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCGTGGCAAAAAGAGCCAGACCTTAACCTTCTGCTATTTTTTCAGAATCTGGCGAATGAGTCAGGTTTTGACGGCGAACTGGCGTATTAGAc       <  2:2556124/150‑1 (MQ=255)
                                                                                                                                               aGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGTTTTGACGGCGAACTGGCGGATTTGACg      >  2:198988/1‑150 (MQ=255)
                                                                                                                                                   gCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGTTTTGACGGCGAACTGGCGGATTTGACGGATg  <  2:1456483/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TCAAGCAATCCTTTATCTACGCCGCCCATTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGAGCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGTTTTGACGGCGAACTGGCGGATTTGACGGATG  >  CP009273/4035274‑4035570

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: