Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,563,266 A→T 5.4% intergenic (+149/‑326) yffP → / → yffQ CPZ‑55 prophage; uncharacterized protein/CPZ‑55 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,563,2660AT 5.4% 237.2 / 2.4 85intergenic (+149/‑326)yffP/yffQCPZ‑55 prophage; uncharacterized protein/CPZ‑55 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (47/33);  new base T (3/2);  total (50/35)
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.48e-01

TTTTCGGACAAACGCGCGCGTAGCCGACGTCAGCGGCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCCAACAAGATGGGAG  >  NC_000913/2563119‑2563415
                                                                                                                                                   |                                                                                                                                                     
ttttCGGACAAACGCGCGCGTAGCCGACGTCAGCGGCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAAc                                                                                                                                                     <  2:432291/150‑1 (MQ=255)
    cGGACAAACGCGCGCGTAGCCGACGTCAGCGGCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACACTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTAc                                                                                                                                                 >  2:1412297/1‑150 (MQ=255)
     ggACAAACGCGCGCGTAGCCGACGTCAGCGGCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTAcc                                                                                                                                                <  2:1773700/150‑1 (MQ=255)
                cgcgTAGCCGACGTCAGCGGCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTtct                                                                                                                                     >  2:855131/1‑150 (MQ=255)
                 gcgTAGCCGACGTCAGCGGCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTtctc                                                                                                                                    >  2:574837/1‑150 (MQ=255)
                   gTAGCCGACGTCAGCGGCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTcac                                                                                                                                  <  1:523838/150‑1 (MQ=255)
                    tAGCCGACGTCAGCGGCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCGGGTAGGTTGCTAACTGCCTGGGCCtttattttt                                                                                                                                 <  1:1689426/150‑7 (MQ=37)
                    tAGCCGACGTCAGCGGCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTcaca                                                                                                                                 <  2:1314008/150‑1 (MQ=255)
                      gCCGACGTCAGCGGCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAAt                                                                                                                               >  1:1594132/1‑150 (MQ=255)
                       ccGACGTCAGCGGCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATc                                                                                                                              >  2:448593/1‑150 (MQ=255)
                        cGACGTCAGCGGCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATcc                                                                                                                             <  1:476136/150‑1 (MQ=255)
                        cGACGTCAGCGGCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATcc                                                                                                                             <  2:1830205/150‑1 (MQ=255)
                            gTCAGCGGCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGc                                                                                                                         >  1:834772/1‑150 (MQ=255)
                              cAGCGGCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTg                                                                                                                       <  2:1412882/150‑1 (MQ=255)
                                 cGGCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCAAAATCCCCGCTGACa                                                                                                                    >  1:394393/1‑150 (MQ=255)
                                   gCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACAtt                                                                                                                  >  1:277598/1‑150 (MQ=255)
                                   gCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCGGCAAATTTTCCAGATTTCACGTACCCGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCACCGTTGACagt                                                                                                                  >  2:971674/1‑148 (MQ=255)
                                    cAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTCCCCGGGCCTTTCTCACAATCCCCGCTGCCATta                                                                                                                 >  1:1262758/1‑150 (MQ=255)
                                    cAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATta                                                                                                                 <  1:558986/150‑1 (MQ=255)
                                       gACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTGTCACAATCCCCGCTGACATTatat                                                                                                              >  1:647602/1‑150 (MQ=255)
                                          tGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCAc                                                                                                           <  1:583521/150‑1 (MQ=255)
                                            gCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCAccc                                                                                                         <  1:1099798/150‑1 (MQ=255)
                                              ggTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTCCACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGCCATTATATCACCCAt                                                                                                       >  1:951111/1‑150 (MQ=255)
                                               gTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATc                                                                                                      >  2:615586/1‑150 (MQ=255)
                                               gTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATc                                                                                                      >  2:567413/1‑150 (MQ=255)
                                                 cATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTACCCGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGAAGTTAGCAAACTACCCGGTCCTTTCTCACAATCGCCACTGACATTATATCACCAATCtt                                                                                                    >  2:1394877/1‑150 (MQ=255)
                                                 cATTTCCCCGCTCCTACGCGAAAAAACCCTCACAACTACCCTGCAATTTTTCCAGATTTCACGTACCCGCGCCAAAAATAAAGGCCCAGGCAGTTTGCATCCTACGCGGGCCTTTCTCACGATCCCCGCTGACATTATATCTCCCATCtt                                                                                                    >  2:1209353/1‑150 (MQ=255)
                                                    ttCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTcc                                                                                                 >  2:539089/1‑150 (MQ=255)
                                                         cGCTCCTACGCGAAAAATCCCTCACAACTACAATGCAAATTTTGCAGATTTCACGTACAAGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCTGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGtttt                                                                                            <  2:1459380/150‑1 (MQ=255)
                                                           cTCCTACGAGAAAAATCACTCAAAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATACAGGCCAAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCa                                                                                          <  2:1340601/150‑1 (MQ=255)
                                                            tCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAg                                                                                         <  1:899400/150‑1 (MQ=255)
                                                            tCCTAAGCGAAAAATCCCTCACAACTATTCTGCAAATTTTCCAGATTTCACGTAAAAGAGCCAAAAATAAAGGCACAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAg                                                                                         <  2:899710/150‑1 (MQ=255)
                                                              cTACGCGAAAAATCCCTCACATCTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCAAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGt                                                                                       >  2:1392845/1‑150 (MQ=255)
                                                              cTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACATACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGt                                                                                       <  1:494563/150‑1 (MQ=255)
                                                                aCGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTgg                                                                                     <  1:352177/150‑1 (MQ=255)
                                                                aCGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTgg                                                                                     >  1:1251621/1‑150 (MQ=255)
                                                                   cgAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCCCGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGgag                                                                                  <  2:1000348/150‑1 (MQ=255)
                                                                   cgAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGgag                                                                                  >  1:528334/1‑150 (MQ=255)
                                                                     aaaaaTCCCTCACAACTACCATGCAACTTTTCCAGATTTCACGTACACGCGCCAAAAATTCAGGCCCACGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGAc                                                                                <  1:277496/150‑1 (MQ=255)
                                                                       aaaTCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTg                                                                              <  2:1448067/150‑1 (MQ=255)
                                                                          tCCCTCACATCTACCTTTCAAATTTTCGAGATTTCACGTACACTCGCCAAACATAAAGCCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATg                                                                           <  2:298369/150‑1 (MQ=255)
                                                                          tCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATg                                                                           >  2:1569466/1‑150 (MQ=255)
                                                                               cacaACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGACAAAAATAAAGGCCCAGTCAGTTAGCAACCTACCCGGGCCTTTCTCACAATACCCGCTGACATTATATCACCCATCTATCCTTTTTCATGTGGTAGACTGATGgccgc                                                                      >  1:942438/1‑150 (MQ=255)
                                                                                acaACTACCCTGCAAATTTTCCAGCATTCACGTACACGCGCCAAAAATCCAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTCTATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGgccgcc                                                                     <  1:992770/150‑1 (MQ=255)
                                                                                 caACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGccgccg                                                                    >  1:1247240/1‑150 (MQ=255)
                                                                                        ccTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGCCTGATGGCCGCCGACGTTGc                                                             >  1:436373/1‑150 (MQ=255)
                                                                                         cTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGca                                                            >  2:426086/1‑150 (MQ=255)
                                                                                         cTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGca                                                            >  2:1510792/1‑150 (MQ=255)
                                                                                         cTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGca                                                            >  2:689705/1‑150 (MQ=255)
                                                                                         cTGCAAATTTTCCAGATTTCACGAACACGCGCCAAAAATAAAGGCCCAGGCAGTTTGCAACCTACCCGGGCCTTTCTCACAATCCCCGCAGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCTCCGACCTTGca                                                            >  2:708945/1‑150 (MQ=255)
                                                                                          tGCAAATTTTCCAGATTTCACGTACACGCGCAAAAAATAAAGGCCCAGGCAGTTAGCAAACTACCCGGGCCTTTATCACAATCCCCGCTGACATTATATCACCAAGCGTTCCGTTTTCAGGTGGTAGACTGATGTCCGCCGCAGTTGcac                                                           >  2:165863/1‑150 (MQ=255)
                                                                                               aatttCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTg                                                      <  1:10429/148‑1 (MQ=255)
                                                                                               aTTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTg                                                      <  1:642520/150‑1 (MQ=255)
                                                                                                  ttCAAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCATTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGAGGGCCGCCGACGTTGCACAATTGAtt                                                   <  2:262452/150‑1 (MQ=255)
                                                                                                   tCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTa                                                  >  1:1328611/1‑150 (MQ=255)
                                                                                                     cAGATTTCACGTACACGCGCCAAAAATAAAGGCCCCGGCAGTTAGCAACCAACCCGGGACTTTCTCAGAATCCCCGCTGACATTATATCACCCATCAGTCCGGTTTCAGGTGGGAGACGGATGGCCGCCGACGTTGCACAATTGaataat                                                >  2:735164/1‑150 (MQ=255)
                                                                                                         tttCTCGTACACGCGCCAATAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCt                                            <  1:1404489/150‑1 (MQ=255)
                                                                                                           gcACGTACACTCGCCAAAAATAATTGCAAATGCATTTATCAACCTACCCGGGCATTTCTCAGAATGCCCGCTGACATTATATTACCCATCTTTCCGTTTTCAGGTGGGAGACTTATGGCCGGCGACGTTGCACAATTGATTAATTTCTGa                                          <  2:281063/149‑1 (MQ=255)
                                                                                                            cACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGaa                                         <  2:1097892/150‑1 (MQ=255)
                                                                                                                tACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGCTGGCCGCCGACGTTGCACAATTGATTAATGTCTGAACAtt                                     >  1:1682599/1‑150 (MQ=255)
                                                                                                                tACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACAtt                                     >  2:645870/1‑150 (MQ=255)
                                                                                                                 tctcgctCCAAATATAAAGGCCCAGGCAGTTAGCTACCTACCCGGTCCTTTCTCACAATCCCCGCTGATATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACATTTGTTTAATTTCTGAACATTc                                    <  2:1277235/147‑1 (MQ=255)
                                                                                                                 acacGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACAtta                                    >  2:1467116/1‑149 (MQ=255)
                                                                                                                  cacGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCa                                   >  2:1292441/1‑150 (MQ=255)
                                                                                                                  cacGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCa                                   <  1:14457/150‑1 (MQ=255)
                                                                                                                    cgcgCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCaaa                                 >  1:589887/1‑150 (MQ=255)
                                                                                                                         cAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCGTTTCTCACAATCCACGCTGACATTATATCACCCATCTTTACGGTTTCAGGTGGGTGACTGATGGCCGCCGACGTTGCACATTTGATTAATTTCTGAACATTCAAATCAAc                            >  1:1659357/1‑150 (MQ=255)
                                                                                                                         cAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAAc                            <  1:539089/150‑1 (MQ=255)
                                                                                                                         cAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTAGATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCaaa                            >  1:1691653/1‑149 (MQ=255)
                                                                                                                           aaaaTAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTc                          <  1:1431173/150‑1 (MQ=255)
                                                                                                                            aaaTAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACACATCTTTACGTTTTCAGGTGGTAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTcc                         >  2:1084452/1‑150 (MQ=255)
                                                                                                                             gagaaaGGCCCGGGTAGGTTGCTAACTGCCTGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTccc                        >  1:546911/4‑150 (MQ=38)
                                                                                                                             gagaaaGGCCCGGGTAGGTTGCTAACTGCCTGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTccc                        >  2:1273934/4‑150 (MQ=38)
                                                                                                                             aaTAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGGTGCACAATTGAGTAATTTCTGAACATTCAAATCAACTccc                        >  1:1898657/1‑150 (MQ=255)
                                                                                                                             aaTAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCCTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTccc                        >  1:612757/1‑150 (MQ=255)
                                                                                                                                   ggCCCGGGTAGGTTGCTAACTGCCTGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCaaa                  >  1:711337/1‑150 (MQ=255)
                                                                                                                                   ggCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCaaa                  >  2:1561996/1‑150 (MQ=255)
                                                                                                                                   ggCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCaaa                  >  2:1029037/1‑150 (MQ=255)
                                                                                                                                       cAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCca              <  1:611882/150‑1 (MQ=255)
                                                                                                                                       cAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCca              >  2:1397827/1‑150 (MQ=255)
                                                                                                                                             gTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCCAACAAGa        >  1:910271/1‑150 (MQ=255)
                                                                                                                                               tAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCCAACAAGATg      >  2:1424495/1‑150 (MQ=255)
                                                                                                                                               tAGCAACCTACCCGGGCCTTTCTCACAAACCCCGCTGACATAATATCACACATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCAAAATTGATTAATTTCTGAACATTCAAATAAACTCCCACCAAAGCCAACAAGAtt      >  1:1681374/1‑149 (MQ=255)
                                                                                                                                                  cAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGTTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCCAACAAGATGGGa   <  1:183678/150‑1 (MQ=255)
                                                                                                                                                  cAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCCAACAAGATGGGa   <  1:1555380/150‑1 (MQ=255)
                                                                                                                                                   aaCCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCCAACAAGATGGGAg  <  1:615586/150‑1 (MQ=255)
                                                                                                                                                   aaCCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCCAACAAGATGGGAg  >  1:175227/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TTTTCGGACAAACGCGCGCGTAGCCGACGTCAGCGGCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCCAACAAGATGGGAG  >  NC_000913/2563119‑2563415

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑