Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 215,395 A→G 9.3% P6967P (CCA→CCG PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451215,3950AG 9.3% 157.3 / 3.1 61P6967P (CCA→CCGPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base A (20/35);  new base G (2/4);  total (22/39)
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.44e-01

GCCGGTCAGCACGTCGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGCCAAGCTGACCGCGACCCCTTCGGTCACCGAAGGTGGCGAGATCACTTACACCATC  >  AE015451/215182‑215611
                                                                                                                                                                                                                     |                                                                                                                                                                                                                        
gCCGGTCAGCACGTCGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGc                                                                                                                                                                                                     <  1:803570/235‑1 (MQ=35)
gCCGGGCAGCACGCCGGTGCCCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTgc                                                                                                                                                                                                            <  1:203728/228‑1 (MQ=17)
  cGGTCAGCACGTCGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGAt                                                                                                                                                                                                  <  1:68502/236‑1 (MQ=35)
          aCGTCGGTGACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCCAAACCAGCGCGCCATac                                                                                                                                                                                                                       >  1:393095/1‑207 (MQ=11)
                             cAGGTACGCCGGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATa                                                                                                                                                                                                                        >  1:157728/1‑187 (MQ=33)
                                          gaaCGAAGGCGACTTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAacg                                                                                                                                                                                     <  2:349319/208‑1 (MQ=38)
                                          cAACGAAGGCGACTTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAacg                                                                                                                                                                                     <  2:361616/209‑1 (MQ=38)
                                                aGGCGACTTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCg                                                                                                                                                    <  1:558082/236‑1 (MQ=38)
                                                      cTTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAGGGCGATGACGttta                                                                                                                                                                                          <  2:547254/192‑4 (MQ=12)
                                                      cTTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAGGGCGATGACGTTTa                                                                                                                                                                                          >  1:547256/1‑192 (MQ=17)
                                                                 tCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGt                                                                                                                                   <  1:272454/236‑1 (MQ=38)
                                                                   aCCATCCCCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCATACAcc                                                                                                                                                                                                                    <  2:258603/153‑1 (MQ=9)
                                                                     cATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGcc                                                                                                                                                                            <  1:624900/191‑1 (MQ=38)
                                                                     cATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAACCCGACCTTCACCGTTCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCg                                                                                                                                                                                                    <  2:143661/167‑1 (MQ=255)
                                                                      aTCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTggg                                                                                                                                                               >  2:685165/1‑203 (MQ=37)
                                                                           cgccgACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCg                                                                                                                                                    <  2:557058/209‑1 (MQ=38)
                                                                             ccgACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGaa                                                                                                                       >  1:718035/1‑236 (MQ=37)
                                                                              cgACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGaaa                                                                                                                      >  1:752841/1‑236 (MQ=37)
                                                                                aCCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGaa                                                                                                                       >  2:270786/1‑233 (MQ=35)
                                                                                    gACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAg                                                                                                                                                                                                       <  1:281750/149‑1 (MQ=35)
                                                                                     aCCTCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAacg                                                                                                                                                                                     <  2:393093/166‑1 (MQ=255)
                                                                                          ggTGGCCGAGAACGTCAAGCCGACCTTCACCGTTCACGTCAACCCGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCt                                                                                                                                                                                                              <  2:111150/136‑1 (MQ=21)
                                                                                               ccGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCg                                                                                                       <  1:594632/234‑1 (MQ=17)
                                                                                                 gagaACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGc                                                                                                   <  1:695399/236‑1 (MQ=34)
                                                                                                  agaACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGc                                                                                                        >  2:342531/1‑230 (MQ=34)
                                                                                                                aCCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTcc                                                                                                                                                     >  1:361144/1‑171 (MQ=34)
                                                                                                                aCCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTcc                                                                                                                                                     <  2:361142/171‑1 (MQ=34)
                                                                                                                     ccccGGGCACGGCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCCGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGAt                                                                                                                                                                                                  <  2:243585/119‑1 (MQ=255)
                                                                                                                         gTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACACCGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGAt                                                                                                                                                                                                  <  1:443632/117‑1 (MQ=25)
                                                                                                                             aCGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCt                                                                                                                                                  >  2:152802/1‑161 (MQ=33)
                                                                                                                               gTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGACGTAGCCgg                                                                                                                                    >  1:793131/1‑173 (MQ=21)
                                                                                                                                  aaCCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACcgacga                                                                  <  1:673784/236‑1 (MQ=17)
                                                                                                                                    ccAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCg                                                                                                                                                                                                    <  2:92207/104‑1 (MQ=35)
                                                                                                                                    ccAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTcc                                                                                                                                                     <  2:684710/151‑1 (MQ=33)
                                                                                                                                     cAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACcgacga                                                                  <  2:427369/233‑1 (MQ=33)
                                                                                                                                       gCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACAACCGACGAagta                                                              >  1:4671/1‑235 (MQ=17)
                                                                                                                                       gCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGcc                                                                                                                                                                   >  1:687426/1‑134 (MQ=34)
                                                                                                                                           cTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAg                                                                                                                                                                                                       <  1:653845/94‑1 (MQ=37)
                                                                                                                                           cTGGCCCACGACCGGGGCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCa                                                                                                                                                                           <  2:159608/122‑1 (MQ=25)
                                                                                                                                                     aCCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATa                                                                                                                                                                                                                        >  1:175920/1‑67 (MQ=35)
                                                                                                                                                                  catGAGCAACAGCGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCaa                                                                                      <  1:262601/182‑1 (MQ=18)
                                                                                                                                                                   cTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGc                                                                                                   <  2:483893/170‑1 (MQ=255)
                                                                                                                                                                   cTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGc                                                                                                   >  1:483895/1‑170 (MQ=255)
                                                                                                                                                                     gAGCAACAACGCCCAGGTCACCAGCAAGGCCGGCGAAACCAGCGCGCCATACAACCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCa                                                                                                                                                                           <  2:143701/96‑1 (MQ=16)
                                                                                                                                                                       gCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTg                                                                                                                                                                 >  2:352/1‑104 (MQ=32)
                                                                                                                                                                         aacaacGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAacg                                                                                                                                                                                     >  1:677445/1‑82 (MQ=33)
                                                                                                                                                                          acaacGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGCCAAGCTGACCGCGACCCCTTCGGTCACCGaa                          >  1:138203/1‑236 (MQ=32)
                                                                                                                                                                          acaacGACCAGGTCACCAACAAGGGCGGCGAAACCCGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGCATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATac                                                                           <  2:642011/187‑1 (MQ=255)
                                                                                                                                                                             acGCCCAGGTCACCACCAAGGCCGGCGAAACCAGCGCGCCATACACCCCCGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTg                                                                                                                 <  2:456273/145‑1 (MQ=255)
                                                                                                                                                                               gCCCCGGTCGCCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGGCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCg                                                                                                                                 <  1:146037/128‑1 (MQ=12)
                                                                                                                                                                                cccAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGTGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTg                                                                                                           <  2:763620/149‑1 (MQ=17)
                                                                                                                                                                                cccAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGCCAAGCTGACCGCGACCCCTTCGGTCACCGAAGGTGGc                    <  1:364692/236‑1 (MQ=31)
                                                                                                                                                                                  cAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGCCAAGCTGACCGCGACCCCTTCGGTCACCGAAGGTGGCga                  <  1:793339/236‑1 (MQ=31)
                                                                                                                                                                                      tCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCAAGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGCCAAGCTGACCGCGACCCCTTCGGTCACCGAAGGt                       >  1:553146/1‑227 (MQ=14)
                                                                                                                                                                                        aCCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGCCAAGCTGACCGCGACCCCTTCGGTCACCGAAGGTGGCGAGAt               >  2:682832/1‑233 (MQ=31)
                                                                                                                                                                                           aTCAAGGCTGGCGAAACCAGCGCGCCGTGCACCCACGCTGCGCAGGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCGCGTCCGCTGTGGATg                                                                                                                                           <  2:152411/106‑1 (MQ=11)
                                                                                                                                                                                             cAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCGGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGCCCAAGCTGACcg                                              >  2:612123/1‑197 (MQ=14)
                                                                                                                                                                                                  ccGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGt                                                                                                                                          >  2:662141/1‑100 (MQ=255)
                                                                                                                                                                                                  ccGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGt                                                                                                                                          <  1:662143/100‑1 (MQ=16)
                                                                                                                                                                                                  ccGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGCCAAGCTGACCGCGACCCCTTCGGTCACCGAAGGTGGCGAGATCACTTACACCATc  <  1:268673/236‑1 (MQ=31)
                                                                                                                                                                                                  ccGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGCCAAGCTGACCGCGACCCCTTCGGTCACCGAAGGTGGCGAGATCACTTACACCATc  <  1:691855/236‑1 (MQ=31)
                                                                                                                                                                                                                     |                                                                                                                                                                                                                        
GCCGGTCAGCACGTCGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGCCAAGCTGACCGCGACCCCTTCGGTCACCGAAGGTGGCGAGATCACTTACACCATC  >  AE015451/215182‑215611

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