Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,118,532 A→G 5.9% intergenic (‑128/+148) 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,5320AG 5.9% 263.9 / 2.6 79intergenic (‑128/+148)wcaK/wzxCcolanic acid biosynthesis protein/putative colanic acid exporter
Reads supporting (aligned to +/- strand):  ref base A (34/40);  new base G (3/2);  total (37/42)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCAGAATAAGTAATTTCATTTTTTCCTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACGGATAAGGCGTTTTACGCCGCATCCGGCAACCGTTGTCGGAACCGAAAAC  >  NC_000913/2118385‑2118673
                                                                                                                                                   |                                                                                                                                             
cccAGAATAAGTAATTTCATTTTTTCCTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATa                                                                                                                                               <  3:233929/148‑1 (MQ=255)
  cAGAATAAGTAATTTCATTTTTTCCTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAgg                                                                                                                                            <  1:151628/149‑1 (MQ=255)
     aaTAAGTAATTTCATTTTTTCCTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCg                                                                                                                                          <  4:35802/148‑1 (MQ=255)
     aaTAAGTAATTTCATTTTTTCCTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGt                                                                                                                                         <  2:387572/149‑1 (MQ=255)
     aaTAAGTAATTTCATTTTTTCCTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGt                                                                                                                                         >  1:778982/1‑149 (MQ=255)
            aaTTTCATTTTTTCCTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCAcgc                                                                                                                                   <  3:641493/148‑1 (MQ=255)
              tttCATTTTTTCCTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAg                                                                                                                                             >  2:327853/1‑136 (MQ=255)
              tttCATTTTTTCCTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAg                                                                                                                                             <  1:327853/136‑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                                                                                                                               >  2:426194/1‑134 (MQ=25)
                           tCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATGCGGCGTGAAcgcctta                                                                                                                               <  1:426194/137‑4 (MQ=25)
                           tCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAg                                                                                                                   >  1:967689/1‑149 (MQ=255)
                            cATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGTACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAg                                                                                                                   >  4:303278/1‑148 (MQ=255)
                            cATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATGCGGCGTGAAcgcctta                                                                                                                               >  3:705137/1‑133 (MQ=25)
                            cATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATGCGGCGTGAAcgcctta                                                                                                                               <  4:705137/136‑4 (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                                                                                                                <  3:776332/148‑1 (MQ=255)
                               aaaTTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGc                                                                                                                <  4:106518/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)
                                                       gtttgtttTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAg                                                                                                                                             >  2:1013215/1‑95 (MQ=255)
                                                       gtttgtttTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAg                                                                                                                                             <  1:1013215/95‑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                                                                                >  1:508129/1‑149 (MQ=255)
                                                              ttgttCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCt                                                                                <  2:216703/149‑1 (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                                                                     <  1:634392/149‑1 (MQ=255)
                                                                         cGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGa                                                                     <  2:182796/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                                                                               <  1:500659/128‑1 (MQ=255)
                                                                                    tGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTa                                                                               >  2:500659/1‑128 (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:300965/1‑149 (MQ=255)
                                                                                          gACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGAc                                                    >  2:61895/1‑149 (MQ=255)
                                                                                           aCGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACg                                                   <  2:145828/149‑1 (MQ=255)
                                                                                           aCGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACg                                                   <  1:61895/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                                                        >  1:506646/1‑126 (MQ=255)
                                                                                                             aGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTa                                                        <  2:506646/126‑1 (MQ=255)
                                                                                                                  tACAAGACCCGAGCTCAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACGGATAAGGCGTTTTACGCCGCATc                            >  2:982390/1‑149 (MQ=255)
                                                                                                                        aCCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACGGATAAGGCGTTTTACGCCGCATCCGGCa                       <  4:114184/148‑1 (MQ=255)
                                                                                                                        aCCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACGGATAAGGCGTTTTACGCCGCATCCGGCa                       >  3:486277/1‑148 (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)
                                                                                                                                                   |                                                                                                                                             
CCCAGAATAAGTAATTTCATTTTTTCCTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACGGATAAGGCGTTTTACGCCGCATCCGGCAACCGTTGTCGGAACCGAAAAC  >  NC_000913/2118385‑2118673

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