Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,181,761 A→G 5.8% intergenic (‑36/+22) ymfA ← / ← potD DUF3592 family inner membrane protein/spermidine/putrescine 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 5.8% 204.2 / 3.3 69intergenic (‑36/+22)ymfA/potDDUF3592 family inner membrane protein/spermidine/putrescine ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base A (29/36);  new base G (2/2);  total (31/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.25e-01

ATCGCGTCTTTGCTATAGGTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACT  >  NC_000913/1181613‑1181909
                                                                                                                                                    |                                                                                                                                                    
aTCGCGTCTTTGCTATAGGTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAAt                                                                                                                                                     >  1:455330/1‑150 (MQ=255)
 tCGCGTCTTTGCTATAGGTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATa                                                                                                                                                    >  2:172534/1‑150 (MQ=255)
     gTCTTTGCTATAGGTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATaccac                                                                                                                                                <  1:1059112/150‑1 (MQ=255)
     gTCTTTGCTATAGGTGAACTGGGTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATaccac                                                                                                                                                <  1:1064468/150‑1 (MQ=255)
      tCTTTGCTATAGGTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAACCGGGCGGTAAAaccacc                                                                                                                                               >  1:975287/1‑150 (MQ=255)
      tCTTTGCTATAGGTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATaccacc                                                                                                                                               <  1:111272/150‑1 (MQ=255)
      tCTATGCTATAGGTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATaccacc                                                                                                                                               <  2:425380/150‑1 (MQ=255)
         ttGCTATAGGTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCAccgcc                                                                                                                                            <  2:603668/150‑1 (MQ=255)
         ttGCTATAGGTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCAccgcc                                                                                                                                            >  1:1184160/1‑150 (MQ=255)
         ttGCTATAGGTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACAAccgcc                                                                                                                                            >  1:308836/1‑150 (MQ=255)
          tGCTATAGGTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACAGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCAccgccc                                                                                                                                           <  2:198317/150‑1 (MQ=255)
            cTATAGGTGAATTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCAccgcccgc                                                                                                                                         >  2:1128555/1‑150 (MQ=255)
            cTATAGGTGAACGGGTAATAAAAGACAACAACATAAATGAGGATCAGTACCAGACCGATACCGAGGAATAAACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTGGTATTAccgcccgt                                                                                                                                         <  1:958275/150‑2 (MQ=255)
               taGGGGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTg                                                                                                                                      <  1:136241/150‑1 (MQ=255)
                aGGTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGc                                                                                                                                     <  2:961552/150‑1 (MQ=255)
                aGGTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGACGAATATACGAAAGAAAACTTTACTGGCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGc                                                                                                                                     <  1:56359/150‑1 (MQ=255)
                 ggTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAGAAACGGGCGGTGGTATTACCGCCCGtttttt                                                                                                                                    <  1:339921/150‑4 (MQ=255)
                  gTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTg                                                                                                                                   <  1:552398/150‑1 (MQ=255)
                    gAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGaa                                                                                                                                 >  2:5854/1‑150 (MQ=255)
                          gTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACg                                                                                                                           <  2:1133924/150‑1 (MQ=255)
                           ttATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGt                                                                                                                          <  1:295435/150‑1 (MQ=255)
                            tataAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTc                                                                                                                         >  2:323888/1‑150 (MQ=255)
                             ataAAAGACAACAACAGAAATGAGGATCAGTACCAGACCTATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTcc                                                                                                                        <  1:956627/150‑1 (MQ=255)
                             ataAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTACCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGAcc                                                                                                                        <  1:956630/150‑1 (MQ=255)
                              taAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCt                                                                                                                       >  2:70678/1‑150 (MQ=255)
                                    acaacaacaGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCTCCGCCCGCTTGCTGAATTAACGTCCTGCTTTc                                                                                                                 <  1:1195050/150‑1 (MQ=255)
                                     caacaacaGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCa                                                                                                                <  1:570818/150‑1 (MQ=255)
                                      aacaacaGAAATTAGTATCAGTACCAGACCTATACCGAGGAATATACGAAAGTAAACTTTACTGTCCTTAATCCTTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAg                                                                                                               <  2:362251/150‑1 (MQ=255)
                                       acaacaGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGc                                                                                                              <  1:1061712/150‑1 (MQ=255)
                                              aaaTGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGat                                                                                                       <  2:1231064/150‑1 (MQ=255)
                                              aaaTGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGat                                                                                                       <  2:1162932/150‑1 (MQ=255)
                                              aaaTGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGat                                                                                                       <  2:1157737/150‑1 (MQ=255)
                                                  taGGATCAGTACCAGACCTATTCCGAGGAATATACGAAAGTAAACTTAACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACTTGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGataata                                                                                                   <  1:78696/149‑1 (MQ=255)
                                                     gATCAGTACCAGACCGATACAGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCCTGCTAACTAAAAAACGGGCGTTAATAACAACGCCCGCTTGCTCAATAAACGTCATACTTACAGCTTCGGAAAATActg                                                                                                >  2:819184/1‑149 (MQ=255)
                                                       tCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCtt                                                                                              >  1:240146/1‑150 (MQ=255)
                                                       tCAGTACCAGACCGATACCGAGGAATAAACTAAAAAAAACTTTAAAGTAATCACTCATTAACCACATCCTAGATAACTAAAAAACGGGCGGAAAAACAACCGCCCGCTTGCTGAATTAACGTACAGATATCACCATCTGATAATAcacat                                                                                              >  2:285760/1‑146 (MQ=255)
                                                        cAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTc                                                                                             >  2:944652/1‑150 (MQ=255)
                                                         aGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCa                                                                                            <  1:377820/150‑1 (MQ=255)
                                                               aGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAAtg                                                                                      <  1:167460/150‑1 (MQ=255)
                                                                          gAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTgc                                                                           >  1:544078/1‑150 (MQ=255)
                                                                              aaTATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCTCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCGTCTGATAATTCTCTTCATAAATGCTGCAGGCTGCGccc                                                                       >  2:857539/1‑149 (MQ=255)
                                                                                 ataCGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACg                                                                    <  2:929846/150‑1 (MQ=255)
                                                                                  taCGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGt                                                                   >  1:1057936/1‑150 (MQ=255)
                                                                                     gaaagaaaACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCAAACGTCAt                                                                <  2:1132971/150‑1 (MQ=255)
                                                                                         gaaaACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAAATAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTg                                                            >  1:941085/1‑150 (MQ=255)
                                                                                                  aCTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGcc                                                   <  2:134712/150‑1 (MQ=255)
                                                                                                  aCTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGcc                                                   <  2:934257/150‑1 (MQ=255)
                                                                                                  aCTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGcc                                                   <  2:31732/150‑1 (MQ=255)
                                                                                                       ccTGACTCATTAGACACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGATTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCTTTCTGCAATTCGAAAtttt                                              >  2:818207/1‑150 (MQ=255)
                                                                                                         tGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTa                                            >  1:609177/1‑150 (MQ=255)
                                                                                                            cTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATg                                         <  2:478618/150‑1 (MQ=255)
                                                                                                               aTTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGtt                                      >  2:718995/1‑150 (MQ=255)
                                                                                                               aTTAGCCACATCCTAGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTTCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGtt                                      <  1:962616/150‑1 (MQ=255)
                                                                                                                   gCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCAGTTTTAATGGTTTCAg                                  >  2:675089/1‑150 (MQ=255)
                                                                                                                    ccACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGc                                 <  1:121176/150‑1 (MQ=255)
                                                                                                                     cacaTCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCa                                >  1:924551/1‑150 (MQ=255)
                                                                                                                       caTCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATc                              >  1:582445/1‑150 (MQ=255)
                                                                                                                           cTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCggg                          >  1:353079/1‑150 (MQ=255)
                                                                                                                           cTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCggg                          >  1:358091/1‑150 (MQ=255)
                                                                                                                                 aactaaAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAgt                    <  1:323888/150‑1 (MQ=255)
                                                                                                                                   ctaaAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAAGTAACGTCCTGCTTTCAGCTTCTGCTAATACTCTTCATAAGTGCTGCTGGCTGCGCCAACGTCATTCTGCCAATCGCCATTTTTTATGTTTTCAGCTTCCGGGTAGAgagt                  >  2:792349/1‑150 (MQ=255)
                                                                                                                                        aaaCGGGAGGTAATACGACCGCCCGCTTGCTGAATTACCGTCCTGCATTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTAAt             <  1:500292/150‑1 (MQ=255)
                                                                                                                                          gcGGGCGGTGGTATTACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATGCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCg           >  2:95065/2‑150 (MQ=255)
                                                                                                                                           cGGGCGGTGGTATTACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGt          >  2:727260/1‑150 (MQ=255)
                                                                                                                                           cGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGTGTGTTTTATCGt          <  1:383830/150‑1 (MQ=255)
                                                                                                                                           cGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGt          <  1:423138/150‑1 (MQ=255)
                                                                                                                                            gggCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATAGTCTTCATAAATGCTTCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTATTGGTTTCAGCATCCGGGTAGAGTGTTTTATCGtt         >  2:650873/1‑150 (MQ=255)
                                                                                                                                            gggCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGtt         >  1:1696/1‑150 (MQ=255)
                                                                                                                                                ggTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGcc     >  1:135325/1‑150 (MQ=255)
                                                                                                                                                   aaTACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACt  <  2:66287/150‑1 (MQ=255)
                                                                                                                                                   aaTACCACCGCCCGCTTGCTGAATTAACGTCCTGCTATCAGCTTCTGATAATACTCTACATAAATGCTGCTGGCTGCGCCAACGACATACTGCCAATCCCCAGTTTTAATGGATTACGCATCCGGGTAGAGTGTTTTATCGTTCGCCACt  >  1:1203899/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ATCGCGTCTTTGCTATAGGTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACT  >  NC_000913/1181613‑1181909

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