Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,181,764 C→T 5.9% intergenic (‑39/+19) 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,7640CT 5.9% 382.9 / 7.6 120intergenic (‑39/+19)ymfA/potDDUF3592 family inner membrane protein/spermidine/putrescine ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base C (68/43);  new base T (2/5);  total (71/49)
Fisher's exact test for biased strand distribution p-value = 1.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.15e-01

Only 100 of 122 total aligned reads displayed.

GCGTCTTTGCTATAGGTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACTTC  >  NC_000913/1181616‑1181911
                                                                                                                                                    |                                                                                                                                                   
gcgTCTTTGCTATAGGTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATACTTGCTAACTAAAAAACGGGCGGTAATacc                                                                                                                                                    >  1:1610099/1‑150 (MQ=255)
  gTCTTTGCTATAGGTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATaccac                                                                                                                                                  <  1:1513133/150‑1 (MQ=255)
      ttGCTATAGGTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTGGTATTAccgcc                                                                                                                                              <  1:45681/150‑1 (MQ=255)
        gCTATAGGTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCAccgcccc                                                                                                                                            >  2:1510200/1‑149 (MQ=255)
         cTATAGGTTAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTGGTATTAccgcccgt                                                                                                                                           <  1:369542/150‑2 (MQ=255)
         cTATAGGTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGCAATATACTAAAGAAAAATTTACTGTCCTGACTCATTAGCCAAATCATTGCTAACTAAAAAACGGGCGGTAATACCAccgcccgc                                                                                                                                           >  1:37784/1‑150 (MQ=255)
          tataGGTGAACTGGTTAAAAAAGACAACAACATAAATGAGGATCAGTACAAGACCTATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTAATTAGACACATCCTTGCTAACTAAAAAACGGGCGGTGGTATTAccgcccgtt                                                                                                                                          <  2:1610026/150‑3 (MQ=255)
           ataGGTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTGGTATTACCGCCCGttt                                                                                                                                         <  2:1220385/150‑1 (MQ=255)
           ataGGTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCtt                                                                                                                                         <  1:269500/150‑1 (MQ=255)
              ggTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCt                                                                                                                                      <  1:1013005/150‑1 (MQ=255)
              ggTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCt                                                                                                                                      <  2:1091262/150‑1 (MQ=255)
              ggTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCt                                                                                                                                      >  1:114165/1‑150 (MQ=255)
                tGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGa                                                                                                                                    <  1:1261247/150‑1 (MQ=255)
                   aCTAGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAAtt                                                                                                                                 >  2:1676377/1‑150 (MQ=255)
                      ggTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAAc                                                                                                                              >  1:1614296/1‑150 (MQ=255)
                      ggTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAAc                                                                                                                              >  1:1328375/1‑150 (MQ=255)
                      ggTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAAc                                                                                                                              >  1:1796889/1‑150 (MQ=255)
                                  caacaacaGAAATGAGGATCAGTACGAGACCGATACCGAGGAATATACGTAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCATGCTAACTAAAAATCGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCa                                                                                                                  <  2:185493/150‑1 (MQ=255)
                                  caacaacaGAAATGAGGAGCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGGCCGCTTGCTGAATTAACGTCCTGCTTTCa                                                                                                                  <  2:472218/150‑1 (MQ=255)
                                    acaacaGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGc                                                                                                                >  1:430479/1‑150 (MQ=255)
                                    acaacaGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGc                                                                                                                <  1:1018746/150‑1 (MQ=255)
                                    acaacaGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACAACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGc                                                                                                                >  1:500765/1‑150 (MQ=255)
                                     caacaGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCt                                                                                                               >  2:392707/1‑150 (MQ=255)
                                       acaGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTc                                                                                                             >  1:275649/1‑150 (MQ=255)
                                            aaTGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGata                                                                                                        <  2:406332/150‑1 (MQ=255)
                                             aTGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGataa                                                                                                       >  2:1158370/1‑150 (MQ=255)
                                             aTGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGataa                                                                                                       >  1:103727/1‑150 (MQ=255)
                                                 ggATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATAct                                                                                                   >  2:1787925/1‑150 (MQ=255)
                                                  gATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATActc                                                                                                  <  2:1807782/150‑1 (MQ=255)
                                                  gATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATActc                                                                                                  >  2:1092696/1‑150 (MQ=255)
                                                   aTCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATActct                                                                                                 >  2:851568/1‑150 (MQ=255)
                                                    tCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCtt                                                                                                >  2:948829/1‑150 (MQ=255)
                                                      aGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCa                                                                                              <  1:1387356/150‑1 (MQ=255)
                                                       ttACCAGACCGATACCGAGGAATATAAGAAAGAAAACTTTACTGTACTGACTCATTAGCCACATCCTTGCTAACTAAAAAACTGGCGGTAATACCACCGCCCGCTTGCTTAATTAAATTCATGCTTTCAGCTTCTGATAATACTCTTCAt                                                                                             <  2:623435/149‑1 (MQ=255)
                                                           cAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAAt                                                                                         >  2:1366635/1‑150 (MQ=255)
                                                           cAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAAt                                                                                         <  2:912160/150‑1 (MQ=255)
                                                            aGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAAtg                                                                                        >  1:814153/1‑150 (MQ=255)
                                                            aGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAAtg                                                                                        <  1:189420/150‑1 (MQ=255)
                                                               ccGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAAtgctg                                                                                     >  2:514507/1‑150 (MQ=255)
                                                                cGATACCGAGGAATATACGAAAGAAAACTTTCCTGTCCTGACTCATTAGCCACCTCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAAtgctgc                                                                                    >  1:285343/1‑150 (MQ=255)
                                                                     ccGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTgctggct                                                                               >  2:957260/1‑150 (MQ=255)
                                                                     ccGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTgctggct                                                                               >  2:475313/1‑150 (MQ=255)
                                                                          gAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGcc                                                                          <  1:190179/150‑1 (MQ=255)
                                                                              ataCGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACg                                                                      <  1:137828/150‑1 (MQ=255)
                                                                              ataCGAAAGAAAACTTTACTGTCCTGACTCATTACCCACATCCTTGCTAACTAAAAAACGGGCTGTAATACCACCGCCCGCTTGATGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCGGTCTGCGCCAACg                                                                      <  1:170275/150‑1 (MQ=255)
                                                                               taCGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGt                                                                     >  2:1789280/1‑150 (MQ=255)
                                                                               taCGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGt                                                                     >  2:1782177/1‑150 (MQ=255)
                                                                               taCGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGt                                                                     >  2:1259117/1‑150 (MQ=255)
                                                                                aCGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTc                                                                    <  2:551505/150‑1 (MQ=255)
                                                                                 cGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCa                                                                   >  1:1494979/1‑150 (MQ=255)
                                                                                 cGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGACCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCa                                                                   >  1:1235621/1‑150 (MQ=255)
                                                                                 cGAAAGAAAACTTTACTGTCCTCACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACAACAGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCAGCGCCAACGtaa                                                                   >  2:1179361/1‑148 (MQ=255)
                                                                                    aagaaaACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTc                                                                <  1:1148631/150‑1 (MQ=255)
                                                                                       aaaaCTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGc                                                             <  2:354361/150‑1 (MQ=255)
                                                                                         aaCTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCAAACGTCATTCTGCCa                                                           <  2:114165/150‑1 (MQ=255)
                                                                                          aCTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCAt                                                          >  1:639289/1‑150 (MQ=255)
                                                                                              tACTGTCCTGACTCATTAGCCACATCCTTGCTAACTATAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGc                                                      <  1:1154203/150‑1 (MQ=255)
                                                                                              tACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGc                                                      <  1:342700/150‑1 (MQ=255)
                                                                                                 tGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCAt                                                   >  1:149471/1‑150 (MQ=255)
                                                                                                  gTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTTAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCAtt                                                  <  2:761919/150‑1 (MQ=255)
                                                                                                    ccTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCAtttt                                                >  1:608043/1‑150 (MQ=255)
                                                                                                    ccTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGGCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCAtttt                                                <  1:272474/150‑1 (MQ=255)
                                                                                                      tGACTCAATAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTa                                              <  1:1662615/150‑1 (MQ=255)
                                                                                                        aCTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAAt                                            >  1:336232/1‑150 (MQ=255)
                                                                                                         cTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATg                                           <  2:589868/150‑1 (MQ=255)
                                                                                                           cATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCATTTTTAATGGtt                                        <  1:606978/150‑1 (MQ=255)
                                                                                                            aTTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGtt                                        >  2:1443167/1‑150 (MQ=255)
                                                                                                              tAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTc                                      <  2:275649/150‑1 (MQ=255)
                                                                                                               aGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCCGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCa                                     >  1:120083/1‑150 (MQ=255)
                                                                                                               aGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCa                                     <  2:498061/150‑1 (MQ=255)
                                                                                                                gCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGACCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTTCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAg                                    >  1:1092683/1‑150 (MQ=255)
                                                                                                                      tCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCg                              >  1:1386253/1‑150 (MQ=255)
                                                                                                                        cTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCggg                            >  2:672399/1‑150 (MQ=255)
                                                                                                                          tGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTa                          <  1:1787925/150‑1 (MQ=255)
                                                                                                                          tGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTa                          >  1:1088815/1‑150 (MQ=255)
                                                                                                                             taactaaAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTagag                       >  1:446858/1‑150 (MQ=255)
                                                                                                                                ctaaAAAACGGGCGGTAATAACACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAgtgt                    >  1:1292739/1‑150 (MQ=255)
                                                                                                                                 taaAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGtt                   >  1:887053/1‑150 (MQ=255)
                                                                                                                                 taaAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGtt                   >  1:1423137/1‑150 (MQ=255)
                                                                                                                                  aaaaaaCGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGttt                  >  1:48509/1‑150 (MQ=255)
                                                                                                                                  aaaaaaCGGGCGGTAATACCACCGCCCGCTTGATGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGTTAGAGTGttt                  >  2:1661292/1‑150 (MQ=255)
                                                                                                                                  aaaaaaCGGGCGGTAATAACACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGttt                  <  2:37784/150‑1 (MQ=255)
                                                                                                                                     aaaCGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTAt               >  2:526978/1‑150 (MQ=255)
                                                                                                                                     aaaCGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTAt               <  1:269307/150‑1 (MQ=255)
                                                                                                                                     aaaCGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTAt               <  1:841769/150‑1 (MQ=255)
                                                                                                                                       aCGGGCGGTAATACCACCGCCCGCTTTCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCg             <  2:1328375/150‑1 (MQ=255)
                                                                                                                                           gCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCg         >  2:699567/1‑150 (MQ=255)
                                                                                                                                           gCGGTAATACAACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATTATACTCTTCATAAATGCTGCTGGCTGGGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCg         >  1:943809/1‑150 (MQ=255)
                                                                                                                                             ggTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGcc       <  2:103727/150‑1 (MQ=255)
                                                                                                                                             ggTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGcc       <  1:1676377/150‑1 (MQ=255)
                                                                                                                                             ggTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGcc       <  1:1152264/150‑1 (MQ=255)
                                                                                                                                             ggTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGcc       <  2:1494979/150‑1 (MQ=255)
                                                                                                                                             ggTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGcc       <  1:851568/150‑1 (MQ=255)
                                                                                                                                              gTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCa      >  1:823368/1‑150 (MQ=255)
                                                                                                                                              gTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAAAGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCa      >  2:1678383/1‑150 (MQ=255)
                                                                                                                                                aaTACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACt    <  1:1385840/150‑1 (MQ=255)
                                                                                                                                                  tACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCTTCCGGGTAGAGTGTTTTATCGTTCGCCACTTc  >  1:673052/1‑150 (MQ=255)
                                                                                                                                                  tACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACTTc  >  2:549997/1‑150 (MQ=255)
                                                                                                                                                  tACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACTTc  <  1:1070700/150‑1 (MQ=255)
                                                                                                                                                  tACCACCGCCCGCTTGCGGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCATCATAAATGGTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCAGTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCACCTTc  >  2:950907/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GCGTCTTTGCTATAGGTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACTTC  >  NC_000913/1181616‑1181911

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: