Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 1,844,166 T→A 87.9% K210* (AAG→TAG)  gacS ← Sensor protein GacS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154511,844,1660TA87.9% 151.6 / 8.8 49K210* (AAG→TAG) gacSSensor protein GacS
Reads supporting (aligned to +/- strand):  ref base T (4/2);  new base A (29/14);  total (33/16)
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.50e-01

ACTTGATGCGGCTGGCCTCAAGTGCCTCTTTGCGCGCCATGTCCAGCTCGATGTTCTGGATTTCGATGGTTTCCAGGTTCTGGCGCACGTCTTCGGTGGCCTGGTCGATGCTGTGTTGCAGTTCTTCGTGGGCATTCTGCAGGGTTTCGGCCATGCGGTTGATGCCGTGGGCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTTGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTTGATGGTGCGGCTCATGCGCAGCGCCAACAGGCCGCTGAGTGCCAGGCAGGCGAGGATCAGCAACAGGCTGGTGAACAGGTTGCGGTAGCCGCGCAGCAGCGTGCCGTCGTGCGACAGCTCGATCTCTACCCAGCCCAGCAGGCGCTCGGCTTCGGCCGGCACGGC  >  AE015451/1843935‑1844373
                                                                                                                                                                                                                                       |                                                                                                                                                                                                               
aCTTGATGCGGCTGGCCTCAAGTGCCTCTTTGCGCGCCATGTCCAGCTCGATGTTCTGGATTTCGATGGTTTCCAGGTTCTGGCGCACGTCTTCGGTGGCCTGGTCGATGCTGTGTTGCAGTTCTTCGTGGGCATTCTGCAGGGTTTCGGCCATGCGGTTGATGCCGTGGGCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTagag                                                                                                                                                                                                              >  1:140249/1‑235 (MQ=255)
        cGGCTGGCCTCAAGTGCCTCTTTGCGCGCCATGTCCAGCTCGATGTTCTGGATTTCGATGGTTTCCAGGTTCTGGCGCACGTCTTCGGTGGCCTGGTCGATGCTGTGTTGCAGTTCTTCGTGGGCATTCTGCAGGGTTTCGGCCATGCGGTTGATGCCGTGGGCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTAGAGCTGATTGa                                                                                                                                                                                                      <  1:571392/235‑1 (MQ=255)
        cGGCTGGCCTCAAGTGCCTCTTTGCGCGCCATGTCCAGCTCGATGTTCTGGATTTCGATGGTTTCCAGGTTCTGGCGCACGTCTTCGGTGGCCTGGTCGATGCTGTGTTGCAGTTCTTCGTGGGCATTCTGCAGGGTTTCGGCCATGCGGTTGATGCCGTGGGCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTAGAGCTGATTGAc                                                                                                                                                                                                     >  1:818697/1‑236 (MQ=255)
         ggCTGGCCTCAAGTGCCTCTTTGCGCGCCATGTCCAGCTCGATGTTCTGGATTTCGATGGTTTCCAGGTTCTGGCGCACGTCTTCGGTGGCCTGGTCGATGCTGTGTTGCAGTTCTTCGTGGGCATTCTGCAGGGTTTCGGCCATGCGGTTGATGCCGTGGGCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTAGAGCTGATTGAc                                                                                                                                                                                                     >  1:1209576/1‑235 (MQ=255)
         ggCTGGCCTCAAGTGCCTCTTTGCGCGCCATGTCCAGCTCGATGTTCTGGATTTCGATGGTTTCCAGGTTCTGGCGCACGTCTTCGGTGGCCTGGTCGATGCTGTGTTGCAGTCCTTCGTGGGCATTCTGCAGGGTTTCGGCCATGCGGTTGATGCCGTGGGCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTAGAGCTGATTGAc                                                                                                                                                                                                     >  1:1121347/1‑235 (MQ=255)
                                  cgcCATGTCCAGCTCGATGTTCTGGATTTCGATGGTCTCCAGGTTCTGGCGCACGTCTTCGGTGGCCTGGTCGATGCTGTGTTGCAGTTCTTCGTGGGCATTCTGCAGGGTTTCGGCCATGCGGTTGATGCCGTGGGCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTAGAGCt                                                                                                                                                                                                            <  2:677042/203‑1 (MQ=255)
                                                cGATGTTCTGGATTTCGATGGTTTCCAGGTTCTGGCGGACGTCTTCGGTGGCCTGGTCGATGCTGTGTTGCAGTTCTTCGTGGGCATTCTGCAGGGTTTCGGCCATGCGGTTGATGCCGTGGGCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTAGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTc                                                                                                                                                                          <  2:582871/223‑1 (MQ=255)
                                                 gATGTTCTGGATTTCGATGGTTTCCTGGTTCTGGCGCACGTCTTCGGTGGCCTGGTCGATGCTGTGTTGCAGTTCTTCGTGGGCATTCTGCAGGGTTTCGGCCATGCGGTTGATGCCGTGGGCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTAGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTAGATAGTGCgg                                                                                                                                                             >  1:681561/1‑235 (MQ=255)
                                                 gATGTTCTGGATTTCGATGGTTTCCAGGTTCTGGCGCACGTCTTCGGTGGCCTGGTCGATGCTGTGTTGCAGTTCTTCGTGGGCATTCTGCAGGGTTTCGGCCATGCGGTTGATGCCGTGGGCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTAGAGCTGATTGACCGCGTGCTTGATGCGTTCGATgg                                                                                                                                                                              >  2:1204315/1‑218 (MQ=255)
                                                  aTGTTCTGGATTTCGATGGTTTCCAGGTTCTGGCGCACGTCTTCGGTGGCCTGGTCGATGCTGTGTTGCAGTTCTTCGTGGGCATTCTGCAGGGTTTCGGCCATGCGGTTGATGCCGTGGGCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTTGAGCTGa                                                                                                                                                                                                          >  2:40656/1‑189 (MQ=255)
                                                                gATGGTTTCCAGGTTCTGGCGCACGTCTTCGGTGGCCTGGTCGATGCTGTGTTGCAGTTCTTCGTGGGCATTCTGCAGGGTTTCGGCCATGCGGTTGATGCCGTGGGCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTAGAGCTGTTTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTTGa                                                                                                                                                                     >  1:182001/1‑212 (MQ=255)
                                                                       tCCAGGTTCTGGCGCACGTCTTCGGTGGCCTGGTCGATGCTGTGTTGCAGTTCTTCGTGGGCATTCTGCAGGGTTTCGGCCATGCGGTTGATGCCGTGGGCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTAGAGCTGATTg                                                                                                                                                                                                       <  2:1012533/171‑1 (MQ=255)
                                                                       tCCAGGTTCTGGCGCACGTCTTCGGTGGCCTGGTCGATGCTGTGTTGCAGTTCTTCGTGGGCATTCTGCAGGGTTTCGGCCATGCGGTTGATGCCGTGGGCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTAGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTTGATGGTGCGGCTCATGCGCAGCGCCAACAGGcc                                                                                                                                      <  1:372505/236‑1 (MQ=255)
                                                                                 ggCGCACGTCTTCGGTGGCCTGGTCGATGCTGTGTTGCAGTTCTTCGTGGGCATTCTGCAGGGTTTCGGCCATGCGGTTGATGCCGTGGGCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTAGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTTGATGGTGCGGCTCATGCGCAGCGCCAACAGGCCGCTGAGTGc                                                                                                                             >  1:988456/1‑235 (MQ=255)
                                                                                                     tGGTCGATGCTGTGTTGCAGTTCTTCGTGGGCATTCTGCAGGGTTTCGGCCATGCGGTTGATGCCGTGGGCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTTGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTTGATGGTGCGGCTCATGCGCAGCGCCAACAGGCCGCTGAGTGCCAGGCAGGCGAg                                                                                                                 >  1:462270/1‑227 (MQ=255)
                                                                                                          gATGCTGTGTTGCAGTTCTTCGTGGGCATTCTGCAGGGTTTCGGCCATGCGGTTGCTGCCGTGGGCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTTGCGCTGATTGACCGCGTGCTTGATGCGTTCGCTGGGGTCGTTGAt                                                                                                                                                                    >  2:69699/1‑171 (MQ=255)
                                                                                                          gATGCTGTGTTGCAGTTCTTCGTGGGCATTCTGCAGGGTTTCGGCCATGCGGTTGATGCCGTGGGCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTAGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTTGATGGTGCGGCTCATGCGCAGCGCCAACAGGCCGCTGAGTGCCAGGCAGGCGAGGATCAGCAACAg                                                                                                     >  1:773365/1‑234 (MQ=255)
                                                                                                          gATGCTGTGTTGCAGTTCTTCGTGGGCATTCTGCAGGGTTTCGGCCATGCGGTTGATGCCGTGGGCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTAGAGCTGATTGACCGCGTGCTTGATGCGTTAGATGGGGTCGTTGATGGTGCGGCTCATgcga                                                                                                                                                    >  2:1085825/1‑186 (MQ=255)
                                                                                                               tgtgTTGCAGTTCTTCGTGGGCATTCTGCAGGGTTTCGGCGATGCGGTTGAAGCCGTGGGCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTAGAGCTGATTGa                                                                                                                                                                                                      <  1:339157/132‑1 (MQ=255)
                                                                                                                    tGCAGTTCTTCGTGGGCATTCTGCAGGGTTTCGGCCATGCGGTTGATGCCGTGGGCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTAGAGCTGCTTGACCGCGTGCTTGATGCGTTCGATGGGGTCGtt                                                                                                                                                                       >  2:256857/1‑158 (MQ=255)
                                                                                                                    tGCAGTTCTTCGTGGGCATTCTGCAGGGTTTCGGCCATGCGGTTGATGCCGTGGGCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGGCCTAGAGCTGATTGCCCGCGt                                                                                                                                                                                                >  2:158103/1‑133 (MQ=255)
                                                                                                                            ttcGTGGGCATTCTGCAGGGTTTCGGCCATGCGGTTGATGCCGTGGGCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTAGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTTGATGGTGCGGCTCATGCGCAGCGCCAACAGGCCGCTGAGGGCCAGGCAGGCGAGGATCAGCAACAGGCTGGTGa                                                                                             >  1:240322/1‑224 (MQ=255)
                                                                                                                              cGTGGGCATTCTGCAGGGTTTCGGCCATGCGGTTGATGCCGTGGGCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTAGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTTGATGGTGCGGCTCATGCGCAGCGCCAACAGGCCGCTGAGTGCCAGGCAGGCGAGGATCAGCAACAGGCTGGTGAACAGGTTGCGGTa                                                                                >  1:127845/1‑235 (MQ=255)
                                                                                                                              cGTGGGCATTCTGCAGGGTTTCGGCCATGCGGTTGATGCCGTGGGCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTAGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTTGATGGTGCGGCTCATGCGCAGCGCCAACAGGCCGCTGAGTGCCAGGCAGGCGAGGATCAGCAACAGGCTGGTGAACAGGTTGCGGTa                                                                                >  1:979748/1‑235 (MQ=255)
                                                                                                                                 gggCATTCTGCAGGGTTTCGGCCATGCGGTTGATGCCGTGGGCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTAGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTTGATGGTGCGGCTCATGCGCAGCGCCAACAGGCCGCTGAGTGCCAGGCAGGCGAGGATCAGCAACAGGCTGGGGAACAGGTTGCGGTAGcc                                                                             >  1:86039/1‑235 (MQ=255)
                                                                                                                                 gggCATTCTGCAGGGTTTCGGCCATGCGGTTGATGCCGTGGGCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTAGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTTGATGGTGCGGCCCATgcgcc                                                                                                                                                   >  2:116285/1‑164 (MQ=255)
                                                                                                                                         tGCAGGGTTTCGGCCATGCGGTTGATGCCGTGGGCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTAGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTTGATGGTGCGGCTCATGCGCAGCGCCAACAGGCCGCTGAGTGCCAGGCAGGCGAGGATCAGCAACAGGCTGGTGAACAg                                                                                         >  1:573176/1‑215 (MQ=255)
                                                                                                                                           cAGGGTTTCGGCCATGCGGTTGATGCCGTGGGCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTAGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTTGATGGTGCGGCTCATGCGCAGCGCCAACAGGCCGCTGAGTGCCAGGCAGGCGAGGATCAGCAACAGGCTGGTGAACAGGTTGCGGTAGCCGCGCAGCAGCg                                                                  <  1:987819/236‑1 (MQ=255)
                                                                                                                                                tttCGGCCATGCGGTTGATGCCGTGGGCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTAGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTTGATGGTGCGGCTCATGCGCAGCGCCAACAGGCCGCTGAGTGCCAGGCAGGCGAGGATCAGCAACAGGCTGGTGAACAGGTTGCGGTAGCCGCGCAGCAGCGTGcc                                                              >  1:493231/1‑235 (MQ=255)
                                                                                                                                                 ttCGGCCATGCGGTTGATGCCGTGGGCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTAGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTTGATGGTGCGGCTCATGCGCAGCGCCAACAGGCCGCTGAGTGCCAGGCAGGCGAGGATCAGCAACAGGCTGGTGAACAGGTTGCGGTAGCCGCGCAGCAGCGTGCcgg                                                            >  2:661171/1‑235 (MQ=255)
                                                                                                                                                      ccATGCGGTTGATGCCGTGGGCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTAGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTTGATGGTGCGGCTCATGCGCAGCGCCAACAGGCCGCTGAGTGCCAGGCAGGCGAGGATCAGCAACAGGCTGGTGAACAGGTTGCGGTAGCCGCGCAGCAGCGTGCcg                                                             >  2:676979/1‑230 (MQ=255)
                                                                                                                                                        aTGCGGTTGATGCCGTGGGCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTTGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTTGATGGTGCGGCTCa                                                                                                                                                         <  2:645268/136‑1 (MQ=255)
                                                                                                                                                           cGGTTGATGCCGTGGGCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTAGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTTGATGGTGCGGCTCATGCGCAGCGCCAACAGGCCGCTGAGTGCCAGGCAGGCGAGGATCAGCAACAGGCTGGTGAACAGGTTGCGGTAGCCGCGCAGCAGCGTGCCGTCGTGCGACAg                                                  <  1:512931/236‑1 (MQ=255)
                                                                                                                                                                   gCCGTGGGCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTAGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTTGATGGTGCGGCTCATGCGCAGCGCCAACAGGCCGCTGAGTGCCAGGCAGGCGAGGATCAGCAACAGGCTGGTGAACAGGTTGCGGTAGCCGCGCAGCAGCGTGCCGTCGTGCGACAg                                                  <  2:1007232/228‑1 (MQ=255)
                                                                                                                                                                        gggCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTAGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTTGATGGTGCGGCTCATGCGCAGCGCCAACAGGCCGCTGAGTGCCAGGCAGGCGAGGATCAGCAACAGGCt                                                                                                  >  2:990030/1‑175 (MQ=255)
                                                                                                                                                                           ccAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTAGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTTGATGGTGCGGCTCATGCGCAGCGCCAACAGGCCGCTGAGTGCCAGGCAGGCGAGGATCAGCAACAGGCTGGTGAACAGGTTGCGGTAGCCGCGCAGCAGCGTGCCGTCGTGCGACAGCTCGATCTCTAcccag                                  >  1:1033422/1‑236 (MQ=255)
                                                                                                                                                                                  gTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTAGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTTGATGGTGCGGCTCATGCGCAGCGCCAACAGGCCGCTGAGTGCCAGGCAGGCGAGGATCAGCAACAGGCTGGTGAACAGGTTGCGGTAGCCGCGCAGCAGCGTGCCGTCGTGCGACAGCTCGATCTCTACccagcccagc                            >  1:691407/1‑235 (MQ=255)
                                                                                                                                                                                   tCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTAGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTTGATGGTGCGGCTCATGCGCAGCGCCAACAGGCCGCTGAGTGCCAGGCAGGCGAGGATCAGCAACAGGCTGGTGAACAGGTTGCgg                                                                                  <  2:238926/180‑1 (MQ=255)
                                                                                                                                                                                      aGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTAGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTTGATGGTGCGGCTCATGCGCAGCGCCAACAGGCCGCTGAGTGCcaggcag                                                                                                                      <  2:290840/141‑1 (MQ=255)
                                                                                                                                                                                           aTGGCTGCCCATGGTCGGCAGGCGTTCTTCGGGGTGGCCGTCCTAGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTTGATGGTGCGGCTCATGCGCAGCGCCAACAGGCCGCTGAGTGCCAGGCAGGCGAGGATCAGCAACAGGCTGGTGAACAGGTTGCGGTAGCCGCGCagcagc                                                                   <  2:543471/187‑1 (MQ=255)
                                                                                                                                                                                                tGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTAGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTTGATGGTGCGGCTCATGCGCAGCGCCAACAGGCCGCTGAGTGCCAGGCAGGCGAGGATCAGCAACAGGCTGGTGAACAGGTTGCGGTAGCCGCGCAGCAGCGTGCCGTCGTGCGACAGCTCGATCTCTACCCAGCCCAGCAGGCGCTc                    >  1:570394/1‑229 (MQ=255)
                                                                                                                                                                                                     aTGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTTGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTTGATGGTGCGGCTCATGCGCAGCGCCAACAGGCCGCTGAGTg                                                                                                                              <  2:940418/118‑1 (MQ=255)
                                                                                                                                                                                                     aTGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTAGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTTGATGGTGCGGCTCATGCGCAGCGCCAACAgg                                                                                                                                        >  2:1208701/1‑108 (MQ=255)
                                                                                                                                                                                                     aTGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTAGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTTGATGGTGCGGCTCATGCGCAGCGCCAACAGGCCGCTg                                                                                                                                  >  2:28368/1‑114 (MQ=255)
                                                                                                                                                                                                     aTGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTAGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTTGATGGTGCGGCTCATGCGCAGCGCCAACAGGACGCTGAGTGCCAGGCAGGCGAGGATCAGCAACAGGCTGGTGAACAGGTTGCGGTAGCCGCGCAGCAGCGTGCCGTCGTGCGACAGCTCGATCTCTACCCAGCCCAGCAGGCGCTCGGCTTcggccgg        >  1:793909/1‑236 (MQ=255)
                                                                                                                                                                                                          cGGCAGGCGTTCTTCGAGGTGGCCGTCCTTGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTTGATGGTGCGGCTCATGCGCAGCGCCAACAGGCCGCTGAGTGCCAGGCAGGCGAGGATCAGCAACAGGCTGGTGAACAGGTTGCGGTGGcc                                                                             >  2:343622/1‑162 (MQ=255)
                                                                                                                                                                                                                      ttcGAGGTGGCCGTCCTAGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTTGATGGTGCGGCTCATGCGCAGCGCCAACAGGCCGCTGAGTGCCAGGCAGGCGAGGATCAGCAACAGGCTGGTGAACAGGTTGCGGTAGCCGCGCAGCAGCGTGCCGTCGTGCGACAGCTCGATCTCTACCCAGCCCAGCAGGCGCTCGGCTTCGGCCGGCACGGc  >  2:594443/1‑225 (MQ=255)
                                                                                                                                                                                                                          agctGGCCGTCCTAGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTTGATGGTGCGGCTCATGCGCAGCGCCAACAGGCCGCTGAGTGCCAGGCAGGCGAGGATCAGCAACAGGCTGGTGAACAGGTTGCGGTAGCCGCGCAGCAGCg                                                                  <  1:976164/154‑1 (MQ=255)
                                                                                                                                                                                                                                 cGTCCTAGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTTGATGGTGCGGCTCATGCGCAGCGCCAACAGGCCGCTGAGTGCcaggc                                                                                                                        >  2:4474/1‑96 (MQ=255)
                                                                                                                                                                                                                                    ccTAGAGCTGATTGACCGCGTGCTTGATGCGTTCg                                                                                                                                                                                  <  1:546136/35‑1 (MQ=255)
                                                                                                                                                                                                                                       |                                                                                                                                                                                                               
ACTTGATGCGGCTGGCCTCAAGTGCCTCTTTGCGCGCCATGTCCAGCTCGATGTTCTGGATTTCGATGGTTTCCAGGTTCTGGCGCACGTCTTCGGTGGCCTGGTCGATGCTGTGTTGCAGTTCTTCGTGGGCATTCTGCAGGGTTTCGGCCATGCGGTTGATGCCGTGGGCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTTGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTTGATGGTGCGGCTCATGCGCAGCGCCAACAGGCCGCTGAGTGCCAGGCAGGCGAGGATCAGCAACAGGCTGGTGAACAGGTTGCGGTAGCCGCGCAGCAGCGTGCCGTCGTGCGACAGCTCGATCTCTACCCAGCCCAGCAGGCGCTCGGCTTCGGCCGGCACGGC  >  AE015451/1843935‑1844373

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