Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,181,761 A→G 6.1% intergenic (‑36/+22) ymfA ← / ← potD putative inner membrane protein/spermidine preferential ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,181,7610AG 6.1% 446.4 / 13.1 147intergenic (‑36/+22)ymfA/potDputative inner membrane protein/spermidine preferential ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base A (72/65);  new base G (4/5);  total (76/71)
Fisher's exact test for biased strand distribution p-value = 7.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

Only 100 of 148 total aligned reads displayed.

AATCGCGTCTTTGCTATAGGTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCA  >  NC_000913/1181612‑1181907
                                                                                                                                                     |                                                                                                                                                  
 aTCGCGTCTTTGCTATAGGTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAAt                                                                                                                                                   >  1:1693605/1‑150 (MQ=255)
  tCGCGTCTTTGCTATAGGTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAAtn                                                                                                                                                  <  1:709290/150‑2 (MQ=255)
    gcgTCTTTGCTATAGGTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATacc                                                                                                                                                >  2:2463882/1‑150 (MQ=255)
    gcgTCTTTGCTATAGGTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATacc                                                                                                                                                >  2:394998/1‑150 (MQ=255)
       tCTTTGCTATAGGTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATaccacc                                                                                                                                             <  1:56050/150‑1 (MQ=255)
         tttGCTATAGGTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCAccgc                                                                                                                                           >  2:1441002/1‑150 (MQ=255)
           tGCTATAGGTGAACTGGTTATGAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCAccgccc                                                                                                                                         <  2:199211/150‑1 (MQ=255)
            gCTATAGGTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTGGTATTAccgcccg                                                                                                                                        <  1:822001/150‑1 (MQ=255)
            gCTATAGGTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTGGTATTAccgcccgt                                                                                                                                       <  2:679400/150‑2 (MQ=255)
             cTATAGGTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTGGTATTAccgcccgt                                                                                                                                       <  1:582327/150‑2 (MQ=255)
              tataGGTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCt                                                                                                                                      >  1:485970/1‑150 (MQ=255)
              tataGGTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCt                                                                                                                                      >  1:255772/1‑150 (MQ=255)
              tataGGTGAACTGGTTATAAAACTCAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCt                                                                                                                                      <  2:549642/150‑1 (MQ=255)
                  ggTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTGGTATTACCGCCCGtttttt                                                                                                                                  <  2:1033574/150‑4 (MQ=255)
                  ggTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCt                                                                                                                                  >  2:868481/1‑150 (MQ=255)
                  ggTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCt                                                                                                                                  >  2:2368039/1‑150 (MQ=255)
                  ggTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCt                                                                                                                                  >  1:1354866/1‑150 (MQ=255)
                   gTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCAGCGCCCGCTTGCTg                                                                                                                                 <  2:113689/150‑1 (MQ=255)
                   gTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTg                                                                                                                                 <  1:1950665/150‑1 (MQ=255)
                      aaCTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGTAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAAt                                                                                                                              <  1:1650363/150‑1 (MQ=255)
                      aaCTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAAt                                                                                                                              <  2:2155514/150‑1 (MQ=255)
                       aCTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAAtt                                                                                                                             <  1:204974/150‑1 (MQ=255)
                       aCTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCACCCGCTTGCTGAAtt                                                                                                                             >  1:118084/1‑150 (MQ=255)
                           gTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACg                                                                                                                         <  1:477014/150‑1 (MQ=255)
                              ataAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTcc                                                                                                                      >  1:2053705/1‑150 (MQ=255)
                                aaaaGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTg                                                                                                                    >  2:295328/1‑150 (MQ=255)
                                       aacaacaGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAg                                                                                                             <  2:1384801/150‑1 (MQ=255)
                                           acaGAAATGAGGATCAGTACCAGACCGATACCTAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCTGTAATACCACCGCCCGCTTGCTGAATTAACGTCATGCTTTCAGCTTc                                                                                                         <  1:2270978/150‑1 (MQ=255)
                                           acaGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTc                                                                                                         <  2:2327828/150‑1 (MQ=255)
                                            caGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCt                                                                                                        <  1:1693294/150‑1 (MQ=255)
                                                aaTGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGata                                                                                                    >  2:2180487/1‑150 (MQ=255)
                                                   gAGGATCAGTACCAGACCGATACCGAGTAATATACGAAAGAAAACTTTACTGTCCTGAATCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGGTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGataata                                                                                                 <  2:627933/150‑1 (MQ=255)
                                                    aGGATCAGTACCAGACCGATACCGAGCAAAATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCAGCTTTCAGCTTCTGATAATAc                                                                                                <  1:2127757/150‑1 (MQ=255)
                                                     ggATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATAct                                                                                               >  1:523007/1‑150 (MQ=255)
                                                     ggATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATAct                                                                                               >  1:1581438/1‑150 (MQ=255)
                                                       aTCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATActct                                                                                             >  1:2327081/1‑150 (MQ=255)
                                                       aTCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAAATTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAGAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATActct                                                                                             >  2:1562295/1‑150 (MQ=255)
                                                           gTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCATCTGATAATACTCTTCAt                                                                                         >  1:1094776/1‑150 (MQ=255)
                                                            tACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATa                                                                                        >  2:1274305/1‑150 (MQ=255)
                                                            tACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATa                                                                                        <  2:1418604/150‑1 (MQ=255)
                                                               cAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATAACACCGCCAGCTTGCTGAATTAACGTAATGCTTTCAGCTTCTTATAATACTCTTAATAAAt                                                                                     >  1:1326594/1‑150 (MQ=255)
                                                                aGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAAtg                                                                                    >  2:707265/1‑150 (MQ=255)
                                                                   ccGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAAtgctg                                                                                 >  1:1313089/1‑150 (MQ=255)
                                                                     gATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATgctgct                                                                               >  2:1303655/1‑150 (MQ=255)
                                                                     gATACCGAGGAATATACGAAAGAAAACTTTACTGTCATGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGATGAATTAAAGTAATGCTTTCAGCTTCTGATAATACTCTTCATAAATtctgct                                                                               >  1:950310/1‑150 (MQ=255)
                                                                           gAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTgc                                                                         >  1:897054/1‑150 (MQ=255)
                                                                               aaTATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCa                                                                     >  2:1126139/1‑150 (MQ=255)
                                                                                  ataCGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTACTGATAATACTCTTCATAAATGCTACTGGCTGCGCCAACg                                                                  >  2:2287749/1‑150 (MQ=255)
                                                                                   taCGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGt                                                                 <  2:2328906/150‑1 (MQ=255)
                                                                                    aCGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTc                                                                <  1:2368039/150‑1 (MQ=255)
                                                                                     cGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCa                                                               >  2:7452/1‑150 (MQ=255)
                                                                                      gaaagaaaACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCAt                                                              >  1:2153864/1‑150 (MQ=255)
                                                                                      gaaagaaaACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCAt                                                              <  2:2053705/150‑1 (MQ=255)
                                                                                        aagaaaACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTc                                                            >  1:1244780/1‑150 (MQ=255)
                                                                                           aaaaCTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGc                                                         <  1:336545/150‑1 (MQ=255)
                                                                                            aaaCTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGcc                                                        <  2:786750/150‑1 (MQ=255)
                                                                                             aaCTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCa                                                       <  1:2451211/150‑1 (MQ=255)
                                                                                                    cTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCa                                                >  1:2422342/1‑150 (MQ=255)
                                                                                                      gTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCAtt                                              >  1:1399124/1‑150 (MQ=255)
                                                                                                       tCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCAttt                                             <  2:1092078/150‑1 (MQ=255)
                                                                                                        ccTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCAtttt                                            <  1:707265/150‑1 (MQ=255)
                                                                                                        ccTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCAtttt                                            <  1:1919625/150‑1 (MQ=255)
                                                                                                          tGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTa                                          <  1:1441002/150‑1 (MQ=255)
                                                                                                           gACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTaa                                         >  1:813533/1‑150 (MQ=255)
                                                                                                             cTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATg                                       <  1:1145596/150‑1 (MQ=255)
                                                                                                              tCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCAGTCGCCATTTTTAATgg                                      <  1:1025200/150‑1 (MQ=255)
                                                                                                               cATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGt                                     >  2:1630255/1‑150 (MQ=255)
                                                                                                                   aGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCa                                 <  2:2100824/150‑1 (MQ=255)
                                                                                                                   aGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCa                                 <  1:1918444/150‑1 (MQ=255)
                                                                                                                    gCCACATCCTTGCTAACTAAAAACCGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAg                                >  1:2350190/1‑150 (MQ=255)
                                                                                                                    gCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAg                                >  1:518561/1‑150 (MQ=255)
                                                                                                                    gCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAg                                <  2:1693605/150‑1 (MQ=255)
                                                                                                                    gCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAg                                <  2:513316/150‑1 (MQ=255)
                                                                                                                    gCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAg                                <  1:1113744/150‑1 (MQ=255)
                                                                                                                      cacaTCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTGCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCa                              <  2:1208103/150‑1 (MQ=255)
                                                                                                                      cacaTCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCa                              >  2:1360750/1‑150 (MQ=255)
                                                                                                                      cacaTCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCa                              >  2:1747626/1‑150 (MQ=255)
                                                                                                                          tCCTTGCTAACTAAAGAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCg                          <  1:1562295/150‑1 (MQ=255)
                                                                                                                            cTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCggg                        <  1:1360750/150‑1 (MQ=255)
                                                                                                                              tGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTa                      <  1:1274305/150‑1 (MQ=255)
                                                                                                                               gCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTag                     <  1:1630255/150‑1 (MQ=255)
                                                                                                                                  aactaaAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAgt                  >  1:1035880/1‑150 (MQ=255)
                                                                                                                                    ctaaAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAgtgt                >  1:2019977/1‑150 (MQ=255)
                                                                                                                                    ctaaAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAgtgt                >  1:2305340/1‑150 (MQ=255)
                                                                                                                                     taaAAAACGGGCGGTAATACCACCTCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGtt               <  2:118084/150‑1 (MQ=255)
                                                                                                                                     taaAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGtt               >  1:177984/1‑150 (MQ=255)
                                                                                                                                     taaAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGtt               >  2:1784903/1‑150 (MQ=255)
                                                                                                                                     taaAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGtt               >  1:887285/1‑150 (MQ=255)
                                                                                                                                      aaaaaaCGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATTTTGCTGTCTGCGCTAACGTGATTCTGGAATTCGACATTTTTTATGTTTTCTGCATCAGTGTAGTGTGttt              <  2:391168/150‑1 (MQ=255)
                                                                                                                                      aaaaaaCGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGttt              >  2:1083031/1‑150 (MQ=255)
                                                                                                                                      aaaaaaCGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGttt              >  1:1068968/1‑150 (MQ=255)
                                                                                                                                        caagcGGGCGGTGGTATTACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTa            >  1:72761/5‑150 (MQ=255)
                                                                                                                                        caagcGGGCGGTGGTATTACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTa            >  2:240033/5‑150 (MQ=255)
                                                                                                                                        aaaaCGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTa            <  1:1522820/150‑1 (MQ=255)
                                                                                                                                        aaaaCGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTa            <  1:1531841/150‑1 (MQ=255)
                                                                                                                                         aaaCGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTAt           <  2:1313089/150‑1 (MQ=255)
                                                                                                                                          agcGGGCGGTGGTATTACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATc          >  1:1178103/3‑150 (MQ=255)
                                                                                                                                              ggCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTc      >  1:2199276/1‑150 (MQ=255)
                                                                                                                                                 ggTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGcc   <  2:1277049/150‑1 (MQ=255)
                                                                                                                                                  gTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCa  >  1:1147813/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
AATCGCGTCTTTGCTATAGGTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCA  >  NC_000913/1181612‑1181907

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