Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,118,533 A→C 5.1% intergenic (‑129/+147) 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,5330AC 5.1% 278.0 / 5.3 79intergenic (‑129/+147)wcaK/wzxCcolanic acid biosynthesis protein/putative colanic acid exporter
Reads supporting (aligned to +/- strand):  ref base A (35/40);  new base C (2/2);  total (37/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

CAGAATAAGTAATTTCATTTTTTCCTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACGGATAAGGCGTTTTACGCCGCATCCGGCAACCGTTGTCGGAACCGAAAAC  >  NC_000913/2118387‑2118673
                                                                                                                                                  |                                                                                                                                            
cAGAATAAGTAATTTCATTTTTTCCTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAgg                                                                                                                                            <  1:151628/149‑1 (MQ=255)
   aaTAAGTAATTTCATTTTTTCCTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCg                                                                                                                                          <  4:35802/148‑1 (MQ=255)
   aaTAAGTAATTTCATTTTTTCCTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGt                                                                                                                                         >  1:778982/1‑149 (MQ=255)
   aaTAAGTAATTTCATTTTTTCCTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGt                                                                                                                                         <  2:387572/149‑1 (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                                                                                                                               <  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                                                                                                                               >  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                                                                                                                <  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)
                                                     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:859591/149‑1 (MQ=255)
                                                                    tGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTcc                                                                        <  2:455400/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                                                   <  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                       >  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)
                                                                                                                                                  |                                                                                                                                            
CAGAATAAGTAATTTCATTTTTTCCTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACGGATAAGGCGTTTTACGCCGCATCCGGCAACCGTTGTCGGAACCGAAAAC  >  NC_000913/2118387‑2118673

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: