Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 1,844,166 T→A 14.3% 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,1660TA14.3% 158.4 / 15.1 55K210* (AAG→TAG) gacSSensor protein GacS
Reads supporting (aligned to +/- strand):  ref base T (29/18);  new base A (6/2);  total (35/20)
Fisher's exact test for biased strand distribution p-value = 6.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.22e-01

CGGCTGGCCTCAAGTGCCTCTTTGCGCGCCATGTCCAGCTCGATGTTCTGGATTTCGATGGTTTCCAGGTTCTGGCGCACGTCTTCGGTGGCCTGGTCGATGCTGTGTTGCAGTTCTTCGTGGGCATTCTGCAGGGTTTCGGCCATGCGGTTGATGCCGTGGGCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTTGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTTGATGGTGCGGCTCATGCGCAGCGCCAACAGGCCGCTGAGTGCCAGGCAGGCGAGGATCAGCAACAGGCTGGTGAACAGGTTGCGGTAGCCGCGCAGCAGCGTGCCGTCGTGCGACAGCTCGATCTCTACCCAGCCCAGCAGGCGCTCGGCTTCGGCCGGCACGGCGTCGGTGGCCA  >  AE015451/1843943‑1844384
                                                                                                                                                                                                                               |                                                                                                                                                                                                                          
cGGCTGGCCTCAAGTGCCTCTTTGCGCGCCATGTCCAGCTCGATGTTCTGGATTTCGATGGTTTCCAGGTTCTGGCGCACGTCTTCGGTGGCCTGGTCGATGCTGTGTTGCAGTTCTTCGTGGGCATTCTGCAGGGTTTCGGCCATGCGGTTGATGCCGTGGGCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTTGAGCTGATTGa                                                                                                                                                                                                                 <  1:515940/235‑1 (MQ=255)
cGGCTGGCCTCAAGTGCCTCTTTGCGCGCCATGTCCAGCTCGATGTTCTGGATTTCGATGGTTTCCAGGTTCTGGCGCACGTCTTCGGTGGCCTGGTCGATGCTGTGTTGCAGTTCTTCGTGGGCATTCTGCAGGGTTTCGGCCATGCGGTTGATGCCGTGGGCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTTGAGCTGATTGa                                                                                                                                                                                                                 <  1:192380/235‑1 (MQ=255)
cGGCTGGCCTCAAGTGCCTCTTTGCGCGCCATGTCCAGCTCGATGTTCTGGATTTCGATGGTTTCCAGGTTCTGGCGCACGTCTTCGGTGGCCTGGTCGATGCTGTGTTGCAGTTCTTCGTGGGCATTCTGCAGGGTTTCGGCCATGCGGTTGATGCCGTGGGCCAGTTCATCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTTGAGCTGATTGa                                                                                                                                                                                                                 <  1:1285300/235‑1 (MQ=255)
         tCAAGTGCCTCTTTGCGCGCCATGTCCAGCTCGATGTTCTGGATTTCGATGGTTTCCAGGTTCTGGCGCACGTCTTCGGTGGCCTGGTCGATGCTGTGTTGCAGTTCTTCGTGGGCATTCTGCAGGGTTTCGGCCATGCGGTTGATGCCGTGGGCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTTGAGCTGATTGACCGCGTGCtt                                                                                                                                                                                                       >  1:1366510/1‑236 (MQ=255)
                  tctTTGCGCGCCATGTCCTGCTCGATGTTCTGGATTTCGATGGTTTCCAGGTTCTGGCGCACGTCTTCGGTGGCCTGGTCGATGCTGTGTTGCAGTTCTTCGTGGGCATTCTGCAGGGTTTCGGCCATGCGGTTGATGCCGTGGGCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTAGAGCTGATTGa                                                                                                                                                                                                                 <  2:1175008/217‑1 (MQ=255)
                   ctTTGCGCGCCATGTCCAGCTCGATGTTCTGGATTTCGATGGTTTCCAGGTTCTGGCGCACGTCTTCGGTGGCCTGGTCGATGCTGTGTTGCAGTTCTTCGTGGGCATTCTGCAGGGTTTCGGCCATGCGGTTGATGCCGTGGGCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTTGAGCTGATTGACCGCGTGCTTGATGCGTTc                                                                                                                                                                                              >  1:1084457/1‑235 (MQ=255)
                             cATGTCCAGCTCGATGTTCTGGATTTCGATGGTTTCCAGGTTCTGGCGCACGTCTTCGGTGGCCTGGTCGATGCTGTGTTGCAGTTCTTCGTGGGCATTCTGCAGGGTTTCGGCCATGCGGTTGATGCCGTGGGCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTTGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGt                                                                                                                                                                                   >  1:606419/1‑236 (MQ=255)
                                          aTGTTCTGGATTTCGATGGTTTCCAGGTTCTGGCGCACGTCTTCGGTGGCCTGGTCGATGCTGTGTTGCAGTTCTTCGTGGGCATTCTGCAGGGTTTCGGCCATGCGGTTGATGCCGTGGGCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTAGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTTGATGGTGCGGCt                                                                                                                                                                      >  1:495065/1‑236 (MQ=255)
                                                        ggtGGTTTCCAGGTTCTGGCGCACGTCTTCGGTGGCCTGGTCGATGCTGTGTTGCAGTTCTTCGTGGGCATTCTGCAGGGTTTCGGCCATGCGGTTGATGCCGTGGGCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTTGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTTGATGGTGCGGCTCa                                                                                                                                                                    <  2:957119/222‑1 (MQ=255)
                                                        gATGGTTTCCAGGTTCTGGCGCACGTCTTCGGTGGCCTGGTCGATGCTGTGTTGCAGTTCTTCGTGGGCATTCTGCAGGGTTTCGGCCATGCGGTTGATGCCGTGGGCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTTGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTTGATGGTGCGGCTCATGCGCAGCGcc                                                                                                                                                         >  1:1257719/1‑235 (MQ=255)
                                                                    gTTCTGGCGCACGTCTTCGGTGGCCTGGTCGATGCTGTGTTGCAGTTCTTCGTGGGCATTCTGCAGGGTTTCGGCCATGCGGTTGATGCCGTGGGCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTTGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTTGATGGTGCGGCTCATGCGCAGCGCCAACAGGcc                                                                                                                                                 <  2:606417/231‑1 (MQ=255)
                                                                          gcgcACGTCTTCGGTGGCCTGGTCGATGCTGTGTTGCAGTTCTTCGTGGGCATTCTGCAGGGTTTCGGCCATGCGGTTGATGCCGTGGGCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTTGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTTGATGGTGCGGCTCATGCGCAGCGCCAACAGGCCGCTGAGTGCca                                                                                                                                      >  1:1120196/1‑236 (MQ=255)
                                                                            gcACGTCTTCGGTGGCCTGGTCGATGCTGTGTTGCAGTTCTTCGTGGGCATTCTGCAGGGTTTCGGCCATGCGGTTGATGCCGTGGGCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTTGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTTGATGGTGCGGCTCATGCGCAGCGCCAACAGGCCGCTGAGTGCcag                                                                                                                                     >  1:1058827/1‑235 (MQ=255)
                                                                            gcACGTCTTCGGTGGCCTGGTCGATGCTGTGTTGCAGTTCTTCGTGGGCATTCTGCAGGGTTTCGGCCATGCGGTTGATGCCGTGGGCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTTGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTTGATGGTGCGGCTCATGCGCAGCGCCAACAGGCCGCTGAGTGCcag                                                                                                                                     >  1:522377/1‑235 (MQ=255)
                                                                                   ttCGGTGGCCTGGTCGATGCTGTGTGGCAGTTCTTCGTGGGCATTCTGCAGGGTTTCGGCCATGCGGTTGATGCCGTGGGCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTTGAGCTGATTGACCGCGTGCTTGATGCGTTCg                                                                                                                                                                                             <  2:106611/172‑1 (MQ=255)
                                                                                               gTCGATGCTGTGTTGCAGTTCTTCGTGGGCATTCTGCAGGGTTTCGGCCATGCGGTTGATGCCGTGGGCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTTGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGGTGATGGTGCGGCt                                                                                                                                                                      >  2:245900/1‑183 (MQ=255)
                                                                                                     gCTGTGTTGCAGTTCTTCGTGGGCATTCTGCAGGGTTTCGGCCATGCGGTTGATGCCGTGGGCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTTGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTTGATGGTGCGGc                                                                                                                                                                       >  1:1398787/1‑176 (MQ=255)
                                                                                                       tgtgTTGCAGTTCTTCGTGGGCATTCTGCAGGGTTTCGGCCATGCGGTTGATGCCGTGGGCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTTGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTTGATgg                                                                                                                                                                             <  2:1220207/168‑1 (MQ=255)
                                                                                                        gtgtTGCAGTTCTTCGTGGGCATTCTGCAGGGTTTCGGCCATGCGGTTGATGCCGTGGGCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTTGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTTGATGGTGCGGCTCATGCGCAGCGCCAACAGGcc                                                                                                                                                 <  2:981153/195‑1 (MQ=255)
                                                                                                        gggtTGCAGTTCTTCGTGGGCATTCTGCAGGGTTTCGGCCATGCGGTTGATGCCGTGGGCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTTGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTTGATGGTGCGGCTCATGCGCAGCGCCAACAGGCCGCTGAGTGCCAGGCAGGCGAGGATCAGCAACAGGCTgg                                                                                                           <  2:793522/231‑1 (MQ=255)
                                                                                                            tGCAGTTCTTCGTGGGCATTCTGCAGGGTTTCGGCCATGCGGTTGATGCCGTGGGCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTTGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTTGATGGTGCGGCTCATGCGCAGCGCCAACAGGCCGCTGAGTGCCAGGCAGGCGAGGATCAGCAACAGGCTGGTGAACAg                                                                                                    >  1:647046/1‑236 (MQ=255)
                                                                                                            tGCAGTTCTTCGTGGGCATTCTGCAGGGTTTCGGCCATGCGGTTGATGCCGTGGGCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTGCGAGGTGGCCGTCCTTGAGCTGATTGACCGCGTGCTTGATGCGTTCGATggggcc                                                                                                                                                                                     >  2:238347/1‑153 (MQ=255)
                                                                                                                 ttcttcGTGGGCATTCTGCAGGGTTTCGGCCATGCGGTTGATGCCGTGGGCCAGTTCGTCGAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTTGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTTGATGGTGCGGCTCATGCGCAGCGCCAACAGGCCGCt                                                                                                                                              >  1:730401/1‑189 (MQ=255)
                                                                                                                 ttcttcGTGGGCATTCTGCAGGGTTTCGGCCATGCGGTTGATGCCGTGGGCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTTGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTTGATGGTGCGGCTCATGCGCAGCGCCAACAGGCCg                                                                                                                                                >  1:883000/1‑187 (MQ=255)
                                                                                                                  tcttcGTTGGCATTCTGCAGGGTTTCGGCCATGCGGTTGATGCCGTGGGCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTTGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTTGATGGTGCGGCTCATGCGCAGCGCCAACAGGCCGCTGAGTGCcaggcag                                                                                                                                 <  1:29132/201‑1 (MQ=255)
                                                                                                                      cGTGGGCATTCTGCAGGGTTTCGGCCATGCGGTTGATGCCGTGGGCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTTGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTTGATGGTGCGGCTCATGCGCAGCGCCAACAGGCCGCTGAGTGCCAGGCAGGCGAGGATCAGCAACAGGCTGGTGAACAGGTTGCGGTa                                                                                           >  1:859772/1‑235 (MQ=255)
                                                                                                                      cGTGGGCATTCTGCAGGGTTTCGGCCATGCGGTTGATGCCGTGGGCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTAGAGCTGATTGACCGCGTGCTTGATGCGTTCGAGGGGGTCGTTGATGGTGCGGCTc                                                                                                                                                                     >  2:793801/1‑161 (MQ=255)
                                                                                                                         gggCATTCTGCAGGGTTTCGGCCATGCGGTTGATGCCGTGGGCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTTGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTTGATGGTGCGGCTCATGCGCAGCGCCAACAGGCCGCTGAGTGCCAGGCAGGCGAGGATCAGCAACa                                                                                                                 >  1:598827/1‑210 (MQ=255)
                                                                                                                                   cAGGGTTTCGGCCATGCGGTTGATGCCGTGGGCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTTGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTTGATGGTGCGGCTCATGCGCAGCGCCAACAGGcc                                                                                                                                                 <  2:745099/168‑1 (MQ=255)
                                                                                                                                              ccATGCGGTTGATGCCGTGGGCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTAGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTTGATGGTGCGGCTCATGCGCAGCGCCAACAGGCCGCTGAGTGCCAGGCAGGCGAGGATCAGCAACAGGCTGGTGAACAGGTTGCGGTAGCCGCGCAGCAGCGTGCcgtcg                                                                     >  2:861427/1‑233 (MQ=255)
                                                                                                                                                   cGGTTGATGCCGTGGGCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTTGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTTGATGGTGCGGCTCATGCGCAGCGCCAACAGGCCGCTGAGTGCCAGGCAGGCGAGGATCAGCAACAGGCTGGTGAACAGGTTGCGGTAGCCGCGCAGCAGCGTGCCGTCGTGCGACa                                                              >  1:920579/1‑235 (MQ=255)
                                                                                                                                                   cGGTTGATGCCGTGGGCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTTGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTTGATGGTGCGGCTCATGCGCAGCGCCAACAGGCCGCTGAGTGCCAGGCAGGCGAGGATCAGCAACAGGCTGGTGAACAGGTTGCGGTAGCCGCGCAGCAGCGTGCCGTCGTGCGACa                                                              >  2:1036564/1‑235 (MQ=255)
                                                                                                                                                         aTGCCGTGGGCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCACGCGTTCTTCGAGGTGGCCGTCCTAGAGCTg                                                                                                                                                                                                                      <  2:525148/77‑1 (MQ=255)
                                                                                                                                                                gggCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTTGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTTGATGGTGCGGCTCATGCGCAGCGCCAACAGGCCGCTGAGTGCCAGGCAGGCGAGGATCAg                                                                                                                      <  2:130752/166‑1 (MQ=255)
                                                                                                                                                                gggCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTTGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTTGATGGTGCGGCTCATGCGCAGCGCCAACAGGCCGCTGAGTGCCAGGCAGGCGAGGATCAGCAACAGGCTGGTGAACAGGTTGCGGTAGCCGCGCAGCAGCGTGCCGTCGTGCGACAg                                                             <  2:785401/223‑1 (MQ=255)
                                                                                                                                                                 ggCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTTGAGc                                                                                                                                                                                                                        >  2:160761/1‑67 (MQ=255)
                                                                                                                                                                 ggCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTTGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTTGATGGTGCGGCTCATGCGCAGCGCCAACAGGCCGCTGAGTGCCAGGCAGGCGAGGATCAGCAACAGGCTGGTGAACAGGTTGCGGTAGCCGCGCAGCAGCGTGCCGTCGTGCGACAGCTCGATCTCTAcc                                                >  1:678334/1‑235 (MQ=255)
                                                                                                                                                                  gCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTTGAGCTGATTGACCGCGTGCTTGATGCGTTCg                                                                                                                                                                                             <  2:925527/93‑1 (MQ=255)
                                                                                                                                                                  gCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTTGAGCTGATTGACCGCGTGCTTGATGCGTTCg                                                                                                                                                                                             >  1:925529/1‑93 (MQ=255)
                                                                                                                                                                   ccAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTTGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTTGATGGTGCGGCTCATGCGCAGCGCCAACAGGCCGCTGAGTGCCAGGCAGGCGAGGATCAGCAACAGGCTGGTGAACAGGTTGCGGTAGCCGCGCAGCAGCGTGCCGTCGTGCGACAGCTCGATCTCTAcccag                                             >  1:115368/1‑236 (MQ=255)
                                                                                                                                                                     aGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTAGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTTGATGGTGc                                                                                                                                                                          >  1:1144766/1‑109 (MQ=255)
                                                                                                                                                                      gTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTTGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTTGATGGTGCGGCTCATGCGCAGCGCCAACAGGCCGCTGAGTGCCAGGCAGGCGAGGATCAGCAACAGGCTGGt                                                                                                          <  2:203981/172‑1 (MQ=255)
                                                                                                                                                                              aGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTTGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTTGATGGGGCGGCCCATGCGCAGCGCCAACAGGCCGCTGAGTGCCAGGCAGGCGAGGATCAGCAACAGGCTGGt                                                                                                          >  1:464351/1‑164 (MQ=255)
                                                                                                                                                                                cTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTTGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTTGATGGTGCGGCTCATGCGCAGCGCCAACAGGCCGCTGAGTGCCAGGCAGGCGAGGATCAGCAACAGGCTGGTGAACAGGTTGCGGTAGCCGCGCAGCAGCGTGCCGTCGTGCGACAg                                                             <  2:888276/207‑1 (MQ=255)
                                                                                                                                                                                        tGCCCATGTTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTAGAGCTGATTGCCCGCGTGCTCGATGCGTTCGCTGGGGTCGt                                                                                                                                                                                   >  2:83638/1‑81 (MQ=255)
                                                                                                                                                                                            cATTGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTTGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTTGATGGTGCGGCTCATGCGCAGCGCCAACAGGCCGCTGAGTGCCAGGCAGGCGAGGATCAGCAACAGGCTGGTGAACa                                                                                                     <  2:178333/155‑1 (MQ=255)
                                                                                                                                                                                             aTGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTTGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTTGATGGTGCGGCTCATGCGCAGCGCCAACAGGCCGCTGAGTGCCAGGCAGGCGAGGATCAGCAACAGGCTGGTGAACAGGTTGCggc                                                                                            >  2:917457/1‑162 (MQ=255)
                                                                                                                                                                                             aTGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTTGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTTGATGGTGCGGCTCATGCGCAGCGCCAACAGGCCGCTGAGTGCCAGGCAGGCGAGGATCAGCAACAGGCTGGTGAACAGGTTGCGGTAGCCGCGCAGCAGCGTGCCGTCGTGCGACAGCTCGATCTCTACCCAGCCCAGCAGGCGCTCGGCTTcggccg                    >  2:1262368/1‑235 (MQ=255)
                                                                                                                                                                                             aTGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTTGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTTGATGGTGCGGCTCATGCGCAGCGCCAACAGGCCGCTGAGTGCCAGGCAGGCGAGGATCAGCAACAGGCTGGTGAACAGGTTGCGGTAGCCGCGCAGCAGCGTGCCGTCGTGCGACAGCTCGATCTCTACCCAGCCCAGCAGGCGCTCGGCTTcggccg                    >  1:1402957/1‑235 (MQ=255)
                                                                                                                                                                                             aTGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTTGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTTGATGGTGCGGCTCATGCGCAGCGCCAACAGGCCGCTGAGTGCCAGGCAGGCGAGGATCAGCAACAGGCTGGTGAACAGGTTGCGGTAGCCGCGCAGCAGCGTGCCGTCGTGCGACAGCTCGATCGCTAcccagccca                                         >  1:1087145/1‑214 (MQ=255)
                                                                                                                                                                                             aTGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTTGAGCTGAGTGACCGCGGGCTTGATGCGTTCGATGGGGGCGTTGATGGTGCGGCTCAGGCGCAGCGCCAACAGGCCGCTGAGTGCCAGGCAGGCGGGGATCAGCAACCGGCt                                                                                                             >  2:177190/1‑146 (MQ=255)
                                                                                                                                                                                             aTGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTAGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTTGATGGTGCGGCTCATGCGCAGCGCCAACAGGCCGCTGAGTGCCAGGCAGGCGAGGATCAGCAACAGGCTGGTGAACAGGTTGCGGTAGCCGCGCAGCAGCGTGCCGTCGTGCGACAGCTCGATCTCTACCCAGCCCAGCAGGCGCGCGGc                            >  1:165679/1‑227 (MQ=255)
                                                                                                                                                                                                  cGGCAGGCGTTCTTCGAGGTGGCCGTCCTTGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTTGATGGTGCGGCTCATGCGCAGCGCCAACAGGCCGCTGAGTGCCAGGCAGGCGAGGATCAGCAACAGGCTGGTGAACAGGTTGCGGTAGCCGCGCAGCAGCGTGCCGTCGTGCGACAGCTCGATCTCTACCCAGCCCAGCAGGCGCTCGGCTTCGGCCGGCACg               >  1:49590/1‑235 (MQ=255)
                                                                                                                                                                                                             cttcGAGGTGGCCGTCCTTGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTTGATGGTGCGGCTCATGCGCAGCGCCAACAGGCCGCTGAGTGCCAGGCAGGCGAGGATCAGCAACAGGCTGGTGAACAGGTTGCGGTAGCCGCGCAGCAGCGTGCCGTCGTGCGACAGCTCGATCTCTACCCAGCCCAGCAGGCGCTCg                              <  2:703441/209‑1 (MQ=255)
                                                                                                                                                                                                              ttcGAGGTGGCCGTCCTTGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTTGATGGTGCGGCTCATGCGCAGCGCCAACAGGCCGCTGAGTGCCAGGCAGGCGAGGATCAGCAACAGGCTGGTGAACAGGTTGCGGTAGCCGCGCAGCAGCGTGCCGTCGTGCGACAGCTCGATCTCTACCCAGCCCAGCAGGCGCTCGGCTTCGGCCGGCACGGCGTCGGTGGCCa  >  1:1250405/1‑236 (MQ=255)
                                                                                                                                                                                                                               |                                                                                                                                                                                                                          
CGGCTGGCCTCAAGTGCCTCTTTGCGCGCCATGTCCAGCTCGATGTTCTGGATTTCGATGGTTTCCAGGTTCTGGCGCACGTCTTCGGTGGCCTGGTCGATGCTGTGTTGCAGTTCTTCGTGGGCATTCTGCAGGGTTTCGGCCATGCGGTTGATGCCGTGGGCCAGTTCGTCCAGCTCATGGCTGCCCATGGTCGGCAGGCGTTCTTCGAGGTGGCCGTCCTTGAGCTGATTGACCGCGTGCTTGATGCGTTCGATGGGGTCGTTGATGGTGCGGCTCATGCGCAGCGCCAACAGGCCGCTGAGTGCCAGGCAGGCGAGGATCAGCAACAGGCTGGTGAACAGGTTGCGGTAGCCGCGCAGCAGCGTGCCGTCGTGCGACAGCTCGATCTCTACCCAGCCCAGCAGGCGCTCGGCTTCGGCCGGCACGGCGTCGGTGGCCA  >  AE015451/1843943‑1844384

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