Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,533,940 C→T 13.4% A970T (GCG→ACG)  narZ ← nitrate reductase 2 (NRZ), alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,533,9400CT13.4% 1031.8 / 123.3 335A970T (GCG→ACG) narZnitrate reductase 2 (NRZ), alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (147/142);  new base T (26/19);  total (174/161)
Fisher's exact test for biased strand distribution p-value = 4.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 339 total aligned reads displayed.

TGATCGCTTTCCAGACCGGACCAGGTGGGGCTGGAGATAATTTTACGCGGCTGCGCCTGAATATCGCGAAAGCGAATCTTCTCGTCCTCTTTGTGCAGCGCCAGATGGGTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCTGTGCCGGGCCATCACGCTTGGTGTAATTGAGTTTGCCGAGGAAATCGATTTCATCCTGGG  >  CP009273/1533792‑1534089
                                                                                                                                                    |                                                                                                                                                     
    cGCTTTCCAGACCGGACCAGGTGGGGCTGGAGATAATTTTACGCGGCTGCGCCTGAATATCGCGAAAGCGAATCTTCTCGTCCTCTTTGTGCAGCGCCAGATGGGTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGCTTTGa                                                                                                                                                  <  1:880190/150‑1 (MQ=255)
     gCTTTCCAGACCGGACCAGGTGGGGCTGGAGATAATTTTACGCGGCTGCGCCTGAATATCGCGAAAGCGAATCTTCTCGTCCTCTTTGTGCAGCGCCAGATGGGTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGCTTTGAc                                                                                                                                                 >  2:4173404/1‑150 (MQ=255)
         tCCAGACCGGACCAGGTGGGGCTGGAGATAATTTTACGCGGCTGCGCCTGAATATCGCGAAAGCGAATCTTCTCGTCCTCTTTGTGCAGCGCCAGATGGGTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCa                                                                                                                                             <  1:540439/150‑1 (MQ=255)
            aGACCGGACCAGGTGGGGCTGGAGATAATTTTACGCGGCTGCGCCTGAATATCGCGAAAGCGAATCTTCTCGTCCTCTTTGTGCAGCGCCAGATGGGTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACa                                                                                                                                          <  1:3091169/150‑1 (MQ=255)
              aCCGGACCAGGTGGGGCTGGAGATAATTTTACGCGGCTGCGCCTGAATATCGCGAAAGCGAATCTTCTCGTCCTCTTTGTGCAGCGCCAGATGGGTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACATg                                                                                                                                        <  1:2209313/150‑1 (MQ=255)
                  gACCAGGTGGGGCTGGAGATAATTTTACGCGGCTGCGCCTGAATATCGCGAAAGCGAATCTTCTCGTCCTCTTTGTGCAGCGCCAGATGGGTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCg                                                                                                                                    <  1:4802974/150‑1 (MQ=255)
                  gACCAGGTGGGGCTGGAGATAATTTTACGCGGCTGCGCCTGAATATCGCGAAAGCGAATCTTCTCGTCCTCTTTGTGCAGCGCCAGATGGGTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCg                                                                                                                                    <  1:3012684/150‑1 (MQ=255)
                   aCCAGGTGGGGCTGGAGATAATTTTACGCGGCTGCGCCTGAATATCGCGAAAGCGAATCTTCTCGTCCTCTTTGTGCAGCGCCAGATGGGTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCGt                                                                                                                                   >  1:3203094/1‑150 (MQ=255)
                     cAGGTGGGGCTGGAGATAATTTTACGCGGCTGCGCCTGAATATCGCGAAAGCGAATCTTCTCGTCCTCTTTGTGCAGCGCCAGATGGGTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTg                                                                                                                                 <  2:2495199/150‑1 (MQ=255)
                      aGGTGGGGCTGGAGATAATTTTACGCGGCTGCGCCTGAATATCGCGAAAGCGAATCTTCTCGTCCTCTTTGTGCAGCGCCAGATGGGTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTgg                                                                                                                                >  1:822665/1‑150 (MQ=255)
                         tGGGGCTGGAGATAATTTTACGCGGCTGCGCCTGAATATCGCGAAAGCGAATCTTCTCGTCCTCTTTGTGCAGCGCCAGATGGGTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGttt                                                                                                                             <  1:2377282/150‑1 (MQ=255)
                          ggggCTGGAGATAATTTTACGCGGCTGCGCCTGAATATCGCGAAAGCGAATCTTCTCGTCCTCTTTGTGCAGCGCCAGATGGGTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTc                                                                                                                            <  2:1481645/150‑1 (MQ=255)
                           gggCTGGAGATAATTTTACGCGGCTGCGCCTGAATATCGCGAAAGCGAATCTTCTCGTCCTCTTTGTGCAGCGCCAGATGGGTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGTTTTGACTGCAACATGACCGTTGGTTTcc                                                                                                                           >  1:4406514/1‑150 (MQ=255)
                           gggCTGGAGATAATTTTACGCGGCTGCGCCTGAATATCGCGAAAGCGAATCTTCTCGTCCTCTTTGTGCAGCGCCAGATGGGTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGTTTTGACTGCAACATGACCGTTGGTTTcc                                                                                                                           >  2:3821102/1‑150 (MQ=255)
                           gggCTGGAGATAATTTTACGCGGCTGCGCCTGAATATCGCGAAAGCGAATCTTCTCGTCCTCTTTGTGCAGCGCCAGATGGGTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTcc                                                                                                                           >  1:1530567/1‑150 (MQ=255)
                            ggCTGGAGATAATTTTACGCGGCTGCGCCTGAATATCGCGAAAGCGAATCTTCTCGTCCTCTTTGTGCAGCGCCAGATGGGTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCg                                                                                                                          >  1:167557/1‑150 (MQ=255)
                            ggCTGGAGATAATTTTACGCGGCTGCGCCTGAATATCGCGAAAGCGAAACTTCTCGTCATCTTTGTGCAGCTCCAGATTGGAATGATCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGCTATCACTGCAACATGAACGATGGTTTCCg                                                                                                                          >  2:275226/1‑150 (MQ=255)
                               tGGAGATAATTTTACGCGGCTGCGCCTGAATATCGCGAAAGCGAATCTTCTCGTCCTCTTTGTGCAGCGCCAGATGGGTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGgtg                                                                                                                       <  2:3769935/150‑1 (MQ=255)
                                    aTAATTTTACGCGGCTGCGCCTGAATATCGCGAAAGCGAATCTTCTCGTCCTCTTTGTGCAGCGCCAGATGGGTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGTTTTGACTGCAACATGACCGTTGGTTTCCGgtgccagt                                                                                                                  <  2:410350/150‑1 (MQ=255)
                                    aTAATTTTACGCGGCTGCGCCTGAATATCGCGAAAGCGAATCTTCTCGTCCTCTTTGTGCAGCGCCAGATGGGTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGTTTTGACTGCAACATGACCGTTGGTTTCCGgtgccagt                                                                                                                  <  2:38711/150‑1 (MQ=255)
                                     tAATTTTACGCGGCTGCGCCTGAATATCGCGAAAGCGAATCTTCTCGTCCTCTTTGTGCAGCGCCATATGGGTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGgtgccagtg                                                                                                                 <  1:1401493/150‑1 (MQ=255)
                                     tAATTTTACGCGGCTGCGCCTGAATATCGCGAAAGCGAATCTTCTCGTCCTCTTTGTGCAGCGCCAGATGGGTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGgtgccagtg                                                                                                                 <  2:27175/150‑1 (MQ=255)
                                      aaTTTTACGCGGCTGCGCCTGAATATCGCGAAAGCGAATCTTCTCGTCCTCTTTGTGCAGCGCCAGATGGGTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGgtgccagtgc                                                                                                                <  2:1868419/150‑1 (MQ=255)
                                             cgcgGCTGCGCCTGAATATCGCGAAAGCGAATCTTCTCGTCCTCTTTGTGCAGCGCCAGATGGGTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATc                                                                                                         >  1:3948857/1‑150 (MQ=255)
                                               cgGCTGCGCCTGAATATCGCGAAAGCGAATCTTCTCGTCCTCTTTGTGCAGCGCCAGATGGGTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCAc                                                                                                       >  1:202308/1‑150 (MQ=255)
                                                ggCTGCGCCTGAATATCGCGAAAGCGAATCTTCTCGTCCTCTTTGTGCAGCGCCAGATGGGTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACt                                                                                                      >  2:1377466/1‑150 (MQ=255)
                                                      gcCTGAATATCGCGAAAGCGAATCTTCTCGTCCTCTTTGTGCAGCGCCAGATGGGTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGaa                                                                                                >  1:3986546/1‑150 (MQ=255)
                                                        cTGAATATCGCGAAAGCGAATCTTCTCGTCCTCTTTGTGCAGCGCCAGATGGGTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGc                                                                                              <  2:4723468/150‑1 (MQ=255)
                                                        cTGAATATCGCGAAAGCGAATCTTCTCGTCCTCTTTGTGCAGCGCCAGATGGGTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGc                                                                                              <  2:2365154/150‑1 (MQ=255)
                                                         tGAATATCGCGAAAGCGAATCTTCTCGTCCTCTTTGTGCAGCGCCAGATGGGTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCg                                                                                             <  2:1052008/150‑1 (MQ=255)
                                                           aaTATCGCGAAAGCGAATCTTCTCGTCCTCTTTGTGCAGCGCCAGATGGGTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTc                                                                                           >  2:241214/1‑150 (MQ=255)
                                                            atatCGCGAAAGCGAATCTTCTCGTCCTCTTTGTGCAGCGCCAGATGGGTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCa                                                                                          >  1:2497620/1‑150 (MQ=255)
                                                            atatCGCGAAAGCGAATCTTCTCGTCCTCTTTGTGCAGCGCCAGATGGGTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCa                                                                                          <  1:770238/150‑1 (MQ=255)
                                                                cgcgAAAGCGAATCTTCTCGTCCTCTTTGTGCAGCGCCAGATGGGTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGTTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATgg                                                                                      >  1:3264590/1‑150 (MQ=255)
                                                                  cgAAAGCGAATCTTCTCGTCCTCTTTGTGCAGCGCCAGATGGGTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCg                                                                                    <  1:3735873/150‑1 (MQ=255)
                                                                   gAAAGCGAATCTTCTCGTCCTCTTTGTGCAGCGCCAGATGGGTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCgg                                                                                   >  2:1760948/1‑150 (MQ=255)
                                                                    aaaGCGAATCTTCTCGTCCTCTTTGTGCAGCGCCAGATGGGTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGgt                                                                                  <  1:4319001/150‑1 (MQ=255)
                                                                    aaaGCGAATCTTCTCGTCCTCTTTGTGCAGCGCCAGATGGGTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGgt                                                                                  <  1:1852745/150‑1 (MQ=255)
                                                                    aaaGCGAATCTTCTCGTCCTCTTTGTGCAGCGCCAGATGGGTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGgt                                                                                  >  1:1367680/1‑150 (MQ=255)
                                                                      aGCGAATCTTCTCGTCCTCTTTGTGCAGCGCCAGATGGGTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGgtgt                                                                                <  2:152779/150‑1 (MQ=255)
                                                                        cGAATCTTCTCGTCCTCTTTGTGCAGCGCCAGATGGGTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCACCATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCa                                                                              <  2:3979274/150‑1 (MQ=255)
                                                                         gAATCTTCTCGTCCTCTTTGTGCAGCGCCAGATGGGTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCaa                                                                             <  1:2634747/150‑1 (MQ=255)
                                                                          aaTCTTCTCGTCCTCTTTGTGCAGCGCCAGATGGGTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAAt                                                                            >  2:2351703/1‑150 (MQ=255)
                                                                          aaTCTTCTCGTCCTCTTTGTGCAGCGCCAGATGGGTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAAt                                                                            >  1:1882263/1‑150 (MQ=255)
                                                                           aTCTTCTCGTCCTCTTTGTGCAGCGCCAGATGGGTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATc                                                                           >  2:4022649/1‑150 (MQ=255)
                                                                            tcttctCGTCCTCTTTGTGCAGCGCCAGATGGGTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCa                                                                          <  1:3058437/150‑1 (MQ=255)
                                                                             cttctCGTCCTCTTTGTGCAGCGCCAGATGGGTATGGTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGTTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCAg                                                                         <  1:1255681/150‑1 (MQ=255)
                                                                              ttctCGTCCTCTTTGTGCAGCGCCAGATGGGTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCAGc                                                                        <  1:440803/150‑1 (MQ=255)
                                                                               tctcGTCCTCTTTGTGCAGCGCCAGATGGGTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCAGCg                                                                       >  2:821090/1‑150 (MQ=255)
                                                                               tctcGTCCTCTTTGTGCAGCGCCAGATGGGTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCAGCg                                                                       >  2:2033333/1‑150 (MQ=255)
                                                                                       tctTTGTGCAGCGCCAGATGGGTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCAGCGGACGcccc                                                               >  2:3098086/1‑150 (MQ=255)
                                                                                        ctTTGTGCAGCGCCAGATGGGTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCt                                                              >  2:2813468/1‑150 (MQ=255)
                                                                                        ctTTGTGCAGCGCCAGATGGGTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCt                                                              >  1:3537814/1‑150 (MQ=255)
                                                                                         tttGTGCAGCGCCAGATGGGTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGTTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCtg                                                             >  2:3199668/1‑150 (MQ=255)
                                                                                          ttGTGCAGCGCCAGATGGGTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCtgt                                                            <  1:2100792/150‑1 (MQ=255)
                                                                                           tgtgCAGCGCCAGATGGGTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCtgtg                                                           >  1:5074035/1‑150 (MQ=255)
                                                                                              gCAGCGCCAGATGGGTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCTGTGCCg                                                        <  1:4105656/150‑1 (MQ=255)
                                                                                              gCAGCGCCAGATGGGTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCTGTGCCg                                                        >  2:391969/1‑150 (MQ=255)
                                                                                               cAGCGCCAGATGGGTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCTGTGCCgg                                                       <  2:4913652/150‑1 (MQ=255)
                                                                                                    ccAGATGGGTATGTTCGCGCCCGGTTATCTCGCCCAGCGCCTGCCACGTTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCTGTGCCGGGCCAt                                                  >  1:5097207/1‑150 (MQ=255)
                                                                                                    ccAGATGGGTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCTGTGCCGGGCCAt                                                  >  1:214708/1‑150 (MQ=255)
                                                                                                    ccAGATGGGTATGTTCGCGCCCGGTGAGCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGGGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCTGTGCCGGGCCAt                                                  <  2:4604586/150‑1 (MQ=255)
                                                                                                     cAGATGGGTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCTGTGCCGGGCCATc                                                 >  2:5023695/1‑150 (MQ=255)
                                                                                                      aGATGGGTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCTGTGCCGGGCCATCa                                                <  1:1942579/150‑1 (MQ=255)
                                                                                                       gATGGGTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCTGTGCCGGGCCATCAc                                               <  2:4124054/150‑1 (MQ=255)
                                                                                                       gATGGGTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCTGTGCCGGGCCATCAc                                               >  2:348702/1‑150 (MQ=255)
                                                                                                        aTGGGTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCTGTGCCGGGCCATCACg                                              <  1:33851/150‑1 (MQ=255)
                                                                                                        aTGGGTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATAAGCGGACGCCCCTGTGCCGGGCCATCACg                                              <  2:751340/150‑1 (MQ=255)
                                                                                                           ggTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGTTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCTGTGCCGGGCCATCACGCtt                                           <  2:862800/150‑1 (MQ=255)
                                                                                                            gTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCTGTGCCGGGCCATCACGCTTg                                          <  1:4545090/150‑1 (MQ=255)
                                                                                                            gTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCTGTGCCGGGCCATCACGCTTg                                          <  1:1990234/150‑1 (MQ=255)
                                                                                                              aTGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCTGTGCCGGGCCATCACGCTTGgt                                        >  1:1184348/1‑150 (MQ=255)
                                                                                                              aTGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCTGTGCCGGGCCATCACGCTTGgt                                        >  2:532211/1‑150 (MQ=255)
                                                                                                               tGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGTTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCTGTGCCGGGCCATCACGCTTGgtg                                       >  1:2173676/1‑150 (MQ=255)
                                                                                                               tGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCTGTGCCGGGCCATCACGCTTGgtg                                       <  1:2134771/150‑1 (MQ=255)
                                                                                                                    gcgcCCGGTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCTGTGCCGGGCCATCACGCTTGGTGTAAtt                                  <  2:2180285/150‑1 (MQ=255)
                                                                                                                      gcCCGGTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGACTCACTTCTGAAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCTGTGCCGGGCCATCACGCTTGGTGTACTTGa                                >  2:4411486/1‑150 (MQ=255)
                                                                                                                       cccGGTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCTGTGCCGGGCCATCACGCTTGGTGTAATTGAg                               <  2:3104791/150‑1 (MQ=255)
                                                                                                                           gTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCTGTGCCGGGCCATCACGCTTGGTGTAATTGAGTTTg                           <  2:3441425/150‑1 (MQ=255)
                                                                                                                              aTCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCTGTGCCGGGCCATCACGCTTGGTGTAATTGAGTTTGCCg                        >  1:4684080/1‑150 (MQ=255)
                                                                                                                               tctcGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCTGTGCCGGGCCATCACGCTTGGTGTAATTGAGTTTGCCGa                       >  2:2988135/1‑150 (MQ=255)
                                                                                                                                 tcGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCTGTGCCGGGCCATCACGCTTGGTGTAATTGAGTTTGCCGAgg                     >  1:987573/1‑150 (MQ=255)
                                                                                                                                      cAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCTGTGCCGGGCCATCACGCTTGGTGTAATTGAGTTTGCCGAGGAAATc                <  1:1461700/150‑1 (MQ=255)
                                                                                                                                      cAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCTGTGCCGGGCCATCACGCTTGGTGTAATTGAGTTTGCCGAGGAAATc                <  1:3602129/150‑1 (MQ=255)
                                                                                                                                      cAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCTGTGCCGGGCCATCACGCTTGGTGTAATTGAGTTTGCCGAGGAAATc                <  2:4303426/150‑1 (MQ=255)
                                                                                                                                       aGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCTGTGCCGGGCCATCACGCTTGGTGTAATTGAGTTTGCCGAGGAAATCg               >  1:1072316/1‑150 (MQ=255)
                                                                                                                                       aGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCTGTGCCGGGCCATCACGCTTGGTGTAATTGAGTTTGCCGAGGAAATCg               >  2:2218838/1‑150 (MQ=255)
                                                                                                                                         cgcCTGCCACGTTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCTGTGCCGGGCCATCACGCTTGGTGTAATTGAGTTTGCCGAGGAAATCGAt             >  1:2800985/1‑150 (MQ=255)
                                                                                                                                         cgcCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCTGTGCCGGGCCATCACGCTTGGTGTAATTGAGTTTGCCGAGGAAATCGAt             <  2:1272037/150‑1 (MQ=255)
                                                                                                                                         cgcCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCTGTGCCGGGCCATCACGCTTGGTGTAATTGAGTTTGCCGAGGAAATCGAt             <  2:1620604/150‑1 (MQ=255)
                                                                                                                                          gcCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCTGTGCCGGGCCATCACGCTTGGTGTAATTGAGTTTGCCGAGGAAATCGAtt            <  1:2908804/150‑1 (MQ=255)
                                                                                                                                           ccTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCTGTGCCGGGCCATCACGCTTGGTGTAATTGAGTTTGCCGAGGAAATCGAttt           <  2:5096012/150‑1 (MQ=255)
                                                                                                                                            cTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCTGTGCCGGGCCATCACGCTTGGTGTAATTGAGTTTGCCGAGGAAATCGATTTc          >  2:3344825/1‑150 (MQ=255)
                                                                                                                                             tGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCTGTGCCGGGCCATCACGCTTGGTGTAATTGAGTTTGCCGAGGAAATCGATTTCa         >  1:4978162/1‑150 (MQ=255)
                                                                                                                                              gCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCTGTGCCGGGCCATCACGCTTGGTGTAATTGAGTTTGCCGAGGAAATCGATTTCAt        <  2:111295/150‑1 (MQ=255)
                                                                                                                                               ccACGTTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCTGTGCCGGGCCATCACGCTTGGTGTAATTGAGTTTGCCGAGGAAATCGATTTCATc       <  2:2641097/150‑1 (MQ=255)
                                                                                                                                                  cGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCTGTGCCGGGCCATCACGCTTGGTGTAATTGAGTTTGCCGAGGAAATCGATTTCATCCTg    >  2:4813365/1‑150 (MQ=255)
                                                                                                                                                  cGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCTGTGCCGGGCCATCACGCTTGGTGTAATTGAGTTTGCCGAGGAAATCGATTTCATCCTg    <  2:3726026/150‑1 (MQ=255)
                                                                                                                                                    cTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCTGTGCCGGGCCATCACGCTTGGTGTAATTGAGTTTGCCGAGGAAATCGATTTCATCCTggg  >  1:484648/1‑150 (MQ=255)
                                                                                                                                                    cTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCTGTGCCGGGCCATCACGCTTGGTGTAATTGAGTTTGCCGAGGAAATCGATTTCATCCTggg  <  1:3781137/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TGATCGCTTTCCAGACCGGACCAGGTGGGGCTGGAGATAATTTTACGCGGCTGCGCCTGAATATCGCGAAAGCGAATCTTCTCGTCCTCTTTGTGCAGCGCCAGATGGGTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCTGTGCCGGGCCATCACGCTTGGTGTAATTGAGTTTGCCGAGGAAATCGATTTCATCCTGGG  >  CP009273/1533792‑1534089

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: