Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,118,540 C→A 5.5% intergenic (‑136/+140) wcaK ← / ← wzxC colanic acid biosynthesis protein/putative colanic acid exporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,118,5400CA 5.5% 247.9 / 3.3 72intergenic (‑136/+140)wcaK/wzxCcolanic acid biosynthesis protein/putative colanic acid exporter
Reads supporting (aligned to +/- strand):  ref base C (33/35);  new base A (2/2);  total (35/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATTTCATTTTTTCCTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACGGATAAGGCGTTTTACGCCGCATCCGGCAACCGTTGTCGGAACCGAAAACAGCAACTCACCCCGC  >  NC_000913/2118397‑2118688
                                                                                                                                               |                                                                                                                                                    
aaTTTCATTTTTTCCTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCAcgc                                                                                                                                                  <  3:641493/148‑1 (MQ=255)
  tttCATTTTTTCCTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCAcgccgc                                                                                                                                               >  1:435705/1‑149 (MQ=255)
       ttttttCCTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCg                                                                                                                                          >  1:23454/1‑149 (MQ=255)
          tttCCTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGc                                                                                                                                        <  3:936395/148‑1 (MQ=255)
          tttCCTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCa                                                                                                                                       <  2:435705/149‑1 (MQ=255)
           ttCCTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATCCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAt                                                                                                                                      >  2:76458/1‑149 (MQ=255)
             ccTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTc                                                                                                                                    <  2:54370/149‑1 (MQ=255)
             ccTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTc                                                                                                                                    <  2:23454/149‑1 (MQ=255)
               tCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATGCGGCGTGAAcgcctta                                                                                                                                              <  1:426194/137‑4 (MQ=25)
               tCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATGCGGCGTGAAcgcctta                                                                                                                                              >  2:426194/1‑134 (MQ=25)
               tCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAg                                                                                                                                  >  1:967689/1‑149 (MQ=255)
                cATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGTACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAg                                                                                                                                  >  4:303278/1‑148 (MQ=255)
                cATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATGCGGCGTGAAcgcctta                                                                                                                                              <  4:705137/136‑4 (MQ=25)
                cATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATGCGGCGTGAAcgcctta                                                                                                                                              >  3:705137/1‑133 (MQ=25)
                cATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAg                                                                                                                                  <  4:214494/148‑1 (MQ=255)
                 aTAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTg                                                                                                                                >  2:83685/1‑149 (MQ=255)
                  tAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTg                                                                                                                                >  3:317532/1‑148 (MQ=255)
                  tAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGc                                                                                                                               >  1:365315/1‑149 (MQ=255)
                   aaaTTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGc                                                                                                                               <  4:106518/148‑1 (MQ=255)
                   aaaTTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGc                                                                                                                               <  3:776332/148‑1 (MQ=255)
                       ttGATGCCAGGTGGGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGa                                                                                                                           >  3:76042/1‑148 (MQ=255)
                       ttGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGa                                                                                                                           <  3:367236/148‑1 (MQ=255)
                              cAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATGAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACg                                                                                                                    >  4:931784/1‑148 (MQ=255)
                                ggTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTg                                                                                                                 >  1:145828/1‑149 (MQ=255)
                                      gCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGAcgcg                                                                                                            <  3:638196/148‑1 (MQ=255)
                                       cTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTc                                                                                                          <  1:441497/149‑1 (MQ=255)
                                        tGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTc                                                                                                          <  4:915543/148‑1 (MQ=255)
                                         gCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGATGACGCGTCt                                                                                                         <  4:597232/148‑1 (MQ=255)
                                           gtttgtttTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTAt                                                                                                      >  2:227824/1‑149 (MQ=255)
                                            tttgttttgtTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTAt                                                                                                      >  4:879572/1‑148 (MQ=255)
                                             ttgttttgttCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGTCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCa                                                                                                    <  1:187691/149‑1 (MQ=255)
                                              tgttttgttCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCa                                                                                                    <  3:189740/148‑1 (MQ=255)
                                                ttttgttCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGc                                                                                                 <  2:850365/149‑1 (MQ=255)
                                                  ttgttCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCt                                                                                               <  2:216703/149‑1 (MQ=255)
                                                  ttgttCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCt                                                                                               >  1:508129/1‑149 (MQ=255)
                                                   tgttCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAGCCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTa                                                                                              <  2:992468/149‑1 (MQ=255)
                                                          tGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTcc                                                                                       <  2:455400/149‑1 (MQ=255)
                                                          tGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTcc                                                                                       <  2:859591/149‑1 (MQ=255)
                                                             cGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCg                                                                                     >  4:1000908/1‑148 (MQ=255)
                                                             cGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGa                                                                                    <  2:182796/149‑1 (MQ=255)
                                                             cGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGa                                                                                    <  1:634392/149‑1 (MQ=255)
                                                              gCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGa                                                                                    >  4:181077/1‑148 (MQ=255)
                                                               cTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGc                                                                                  >  2:829628/1‑149 (MQ=255)
                                                                     gCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGaa                                                                            <  1:83685/149‑1 (MQ=255)
                                                                     gCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGa                                                                             <  3:436605/148‑1 (MQ=255)
                                                                       cTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAAcc                                                                          <  1:202078/149‑1 (MQ=255)
                                                                        tGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTa                                                                                              >  2:500659/1‑128 (MQ=255)
                                                                        tGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTa                                                                                              <  1:500659/128‑1 (MQ=255)
                                                                        tGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAAcc                                                                          >  3:1019572/1‑148 (MQ=255)
                                                                         gATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCg                                                                         >  3:921464/1‑148 (MQ=255)
                                                                          aTGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTa                                                                       <  2:778982/149‑1 (MQ=255)
                                                                              gACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGAc                                                                   >  2:61895/1‑149 (MQ=255)
                                                                              gACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGAc                                                                   >  2:300965/1‑149 (MQ=255)
                                                                               aCGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACg                                                                  <  1:61895/149‑1 (MQ=255)
                                                                               aCGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACg                                                                  <  2:145828/149‑1 (MQ=255)
                                                                                cGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACg                                                                  >  3:538812/1‑148 (MQ=255)
                                                                                 gCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGTATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACgg                                                                 <  4:777569/148‑1 (MQ=255)
                                                                                        cgTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACGGATAAGGc                                                          >  3:114184/1‑148 (MQ=255)
                                                                                               tCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACGGATAAGGCGTTTTAc                                                   >  3:904592/1‑148 (MQ=255)
                                                                                                 aGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTa                                                                       <  2:506646/126‑1 (MQ=255)
                                                                                                 aGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTa                                                                       >  1:506646/1‑126 (MQ=255)
                                                                                                      tACAAGACCCGAGCTCAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACGGATAAGGCGTTTTACGCCGCATc                                           >  2:982390/1‑149 (MQ=255)
                                                                                                            aCCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACGGATAAGGCGTTTTACGCCGCATCCGGCa                                      >  3:486277/1‑148 (MQ=255)
                                                                                                            aCCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACGGATAAGGCGTTTTACGCCGCATCCGGCa                                      <  4:114184/148‑1 (MQ=255)
                                                                                                              ccGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACGGATAAGGCGTTTTACGCCGCATCCGGCAAcc                                   >  2:772106/1‑149 (MQ=255)
                                                                                                              ccGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCATGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACGGATAAGGCGTTTTACGCCGCATCCGGCAAcc                                   >  2:712478/1‑149 (MQ=255)
                                                                                                                 aGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACGGATAAGGCGTTTTACGCCGCATCCGGCAACCGt                                 >  3:1020958/1‑148 (MQ=255)
                                                                                                                    acaGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACGGATAAGGCGTTTTACGCCGCATCCGGCAACCGTTGTc                             <  2:388176/149‑1 (MQ=255)
                                                                                                                        aaCCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACGGATAAGGCGTTTTACGCCGCATCCGGCAACCGTTGTCGGaa                         >  1:558341/1‑149 (MQ=255)
                                                                                                                         aCCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACGGATAAGGCGTTTTACGCCGCATCCGGCAACCGTTGTCGGAAc                        <  2:967689/149‑1 (MQ=255)
                                                                                                                                 aCGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACGGATAAGGCGTTTTACGCCGCATCCGGCAACCGTTGTCGGAACCGAAAAc                 <  3:1000908/148‑1 (MQ=255)
                                                                                                                                            gTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACGGATAAGGCGTTTTACGCCGCATCCGGCAACCGTTGTCGGAACCGAAAACAGCAACTCAccc     >  2:403077/1‑149 (MQ=255)
                                                                                                                                               cACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACGGATAAGGCGTTTTACGCCGCATCCGGCAACCGTTGTCGGAACCGAAAACAGCAACTCACcccgc  <  1:982390/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
AATTTCATTTTTTCCTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACGGATAAGGCGTTTTACGCCGCATCCGGCAACCGTTGTCGGAACCGAAAACAGCAACTCACCCCGC  >  NC_000913/2118397‑2118688

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: