Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 195,732 A→C 14.2% N413T (AAC→ACC)  PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451195,7320AC14.2% 264.6 / 20.1 102N413T (AAC→ACC) PP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base A (42/44);  new base C (10/6);  total (52/50)
Fisher's exact test for biased strand distribution p-value = 4.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 103 total aligned reads displayed.

CCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCAC  >  AE015451/195504‑195936
                                                                                                                                                                                                                                    |                                                                                                                                                                                                            
cgaaTGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTc                                                                                                                                                                                                        >  1:978854/3‑235 (MQ=2)
           tACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAAc                                                                                                                                                                                            <  1:1030308/236‑1 (MQ=37)
                      ttCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACCATTACCGtt                                                                                                                                                                                  >  1:1022811/1‑235 (MQ=16)
                           ccGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCGTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGtt                                                                                                                                                                                  <  1:699344/230‑1 (MQ=255)
                           ccGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGtt                                                                                                                                                                                  >  2:699339/1‑230 (MQ=255)
                                                    cacaGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGccg                                                                                                                                                                           <  1:1253905/212‑3 (MQ=1)
                                                    cACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGAtt                                                                                                                                                   >  1:1217132/1‑236 (MQ=255)
                                                     acaggcggTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGc                                                                                                                                                          >  1:939015/4‑228 (MQ=255)
                                                         gcggTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAg                                                                                                                                              <  1:1213388/236‑1 (MQ=255)
                                                                 ttCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCgg                                                                                                                                                                           <  1:117630/199‑1 (MQ=255)
                                                                      gCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATg                                                                                                                                 <  2:665193/236‑1 (MQ=255)
                                                                            gACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAggg                                                                                                                                                              >  1:314550/1‑201 (MQ=40)
                                                                                       cGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCgg                                                                                                                                                                           <  2:923170/177‑1 (MQ=255)
                                                                                       cGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCgg                                                                                                                                                                           >  1:923175/1‑177 (MQ=9)
                                                                                            ccGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAg                                                                                                                                                           >  1:391655/1‑188 (MQ=37)
                                                                                                      ccACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCgg                                                                                                                                                                           <  1:760282/162‑1 (MQ=38)
                                                                                                         cGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAg                                                                                                                                              <  1:1002841/188‑1 (MQ=37)
                                                                                                             cAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGAc                                                                                                                               <  2:600456/199‑1 (MQ=38)
                                                                                                             cAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGAc                                                                                                                               >  1:600461/1‑199 (MQ=37)
                                                                                                             cAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCATCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCgg                                                                                                                                                                           <  1:911188/155‑1 (MQ=255)
                                                                                                              aaCGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATc                                                                                           <  2:1217127/234‑1 (MQ=17)
                                                                                                                  aCTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAAc                                                                                                                                                                                            <  1:905776/133‑1 (MQ=255)
                                                                                                                  aCTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGt                                                                                                                             <  1:814392/196‑1 (MQ=39)
                                                                                                                  aCTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGt                                                                                                                             >  2:814387/1‑196 (MQ=38)
                                                                                                                  aCTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCgaga                                                                                       >  1:458042/1‑234 (MQ=38)
                                                                                                                  aCTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAAc                                                                                                                                                                                            >  2:905771/1‑133 (MQ=19)
                                                                                                                   cTCGGTCGACACCGCCACCGCGACCCTGACGGCGAGCCCGTCGGTCTCCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCgg                                                                                                                                                                           <  2:391650/149‑1 (MQ=255)
                                                                                                                      agTCGACACTACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAAc                                                                                                                                                                                            <  2:1236299/128‑1 (MQ=255)
                                                                                                                      agTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCa                                                                                                                                    <  1:10009/184‑1 (MQ=37)
                                                                                                                        tCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAAc                                                                                                                                                                                            <  1:29719/127‑1 (MQ=40)
                                                                                                                         cGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGCCGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGCTTACCGTTGAAGCGGGTAAAACCCAg                                                                                                                                                                >  2:916835/1‑154 (MQ=255)
                                                                                                                           acacCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGGAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTacaa                                                                                                                       >  1:401238/1‑193 (MQ=37)
                                                                                                                                 ccaccGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCgg                                                                                                                                                                           <  1:106937/135‑1 (MQ=40)
                                                                                                                                  taccGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGcga                                                                                                                                    <  2:939010/172‑3 (MQ=18)
                                                                                                                                            cTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCgg                                                                                                                                                                           >  1:887903/1‑124 (MQ=25)
                                                                                                                                            cTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAAc                                                                                                                                                                                            >  1:1277922/1‑107 (MQ=25)
                                                                                                                                               aCGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCgg                                                                                                                                                                           >  2:304388/1‑121 (MQ=40)
                                                                                                                                               aCGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCgg                                                                                                                                                                           <  1:304392/121‑1 (MQ=40)
                                                                                                                                                       cccGTCGGTCACCGAAGGCGGCGTGCTCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACGATCACCGTTGAAGCGGGAAAAACCCAGAGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACACCGGCTCGa                                                                                                            >  1:875756/1‑176 (MQ=11)
                                                                                                                                                            cGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAActgtct                                            >  1:1263133/1‑230 (MQ=255)
                                                                                                                                                             cgTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAAc                                                                                                                                                                                                    <  2:482700/81‑1 (MQ=255)
                                                                                                                                                              gTCACCGAAGGCGACGTGATCACTTACACCGTACCCTTAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACTGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACa                                                              >  1:679877/1‑214 (MQ=255)
                                                                                                                                                               tCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACgg                                                                                                                 <  1:458036/163‑1 (MQ=39)
                                                                                                                                                                      aGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCgg                                                                                                                                                                           <  1:1035117/98‑1 (MQ=255)
                                                                                                                                                                           gcgTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAAc                                                                                                                                                                                            <  1:552710/76‑1 (MQ=255)
                                                                                                                                                                               gATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGACCCTGt                                                                                                                                                                                                         <  2:1039181/59‑1 (MQ=255)
                                                                                                                                                                               gATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGACCCTGt                                                                                                                                                                                                         >  1:1039186/1‑59 (MQ=255)
                                                                                                                                                                                aTCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATg                                                                                                                                 <  2:1040538/130‑1 (MQ=25)
                                                                                                                                                                                aTCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATg                                                                                                                                 >  1:1040543/1‑130 (MQ=40)
                                                                                                                                                                                   acatACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACg                     <  1:1031076/232‑1 (MQ=255)
                                                                                                                                                                                    cctACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACa                                                              <  2:883489/191‑1 (MQ=37)
                                                                                                                                                                                        cacCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCAGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGCTCGAcc                                                                                                          >  2:506150/1‑145 (MQ=11)
                                                                                                                                                                                        cacCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCg                                                                                          >  2:89657/1‑161 (MQ=38)
                                                                                                                                                                                         acCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGTAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACTGCTCGACCGTCAGCGTCACCATc                                                                                           <  2:936523/159‑1 (MQ=255)
                                                                                                                                                                                         acCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACAGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCg                <  1:701942/234‑1 (MQ=255)
                                                                                                                                                                                              gACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGaa                                                                <  2:1263128/181‑1 (MQ=40)
                                                                                                                                                                                                cccTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGa                                                                                                                                                                                <  1:1033494/67‑1 (MQ=11)
                                                                                                                                                                                                cccTGAGCAACCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGa                                                                                                                                                                                >  1:664557/1‑67 (MQ=25)
                                                                                                                                                                                                cccTGAGCAACCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACGGTCAGCGACAcctc                                                                                            >  2:383868/1‑149 (MQ=9)
                                                                                                                                                                                                cccTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGa                                                                                                                                                                                <  2:664552/67‑1 (MQ=255)
                                                                                                                                                                                                cccGGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCa                                                                                                                                    <  1:302453/111‑1 (MQ=255)
                                                                                                                                                                                                cacTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAg                                                                                                                                              <  2:395339/99‑1 (MQ=255)
                                                                                                                                                                                                  cTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCa                                                                                 >  2:856812/1‑160 (MQ=255)
                                                                                                                                                                                                  cTGAGCAACCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCCCAGACCACGATCAACGACCCGGTc             >  1:404143/1‑228 (MQ=255)
                                                                                                                                                                                                  cTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTacaa                                                                                                                       >  1:1073364/1‑122 (MQ=38)
                                                                                                                                                                                                  cTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTacaa                                                                                                                       <  2:1073359/122‑1 (MQ=38)
                                                                                                                                                                                                  cTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGt                                                                                                        >  1:373183/1‑137 (MQ=40)
                                                                                                                                                                                                  cTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGt                                                                                                        <  2:373178/137‑1 (MQ=40)
                                                                                                                                                                                                    gAGCAACCCTGCCCAGACGCCGGTGACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGTAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTgg                                                                             <  2:897882/162‑1 (MQ=255)
                                                                                                                                                                                                     aGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACCGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGAcaccaccac  <  1:1199664/236‑1 (MQ=255)
                                                                                                                                                                                                        aaCCCTGCCCAGACGCCGGTGACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGccg                                                                                                                                                                           >  2:875988/1‑62 (MQ=255)
                                                                                                                                                                                                        aaCCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACg                                                                                                                                                                                                   >  2:513708/1‑40 (MQ=255)
                                                                                                                                                                                                        aaCCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAGACGATTACCGTTGAAGCGGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGCAACTTCGaa                                                                >  2:743042/1‑171 (MQ=255)
                                                                                                                                                                                                        aaCCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAg                                                                                                                                                                              >  1:326456/1‑61 (MQ=255)
                                                                                                                                                                                                        aaCCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCgg                                                                                                                                                                           <  1:875993/64‑1 (MQ=255)
                                                                                                                                                                                                        aaCCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGt                                                                                                                                                        >  1:1152902/1‑83 (MQ=255)
                                                                                                                                                                                                        aaCCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGcc                                                                                                                                     >  2:333499/1‑102 (MQ=40)
                                                                                                                                                                                                        aaCCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACaaccac                                                                                                                   >  1:232509/1‑120 (MQ=255)
                                                                                                                                                                                                             tGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAACCGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCa                              >  1:1126549/1‑200 (MQ=255)
                                                                                                                                                                                                             tGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTg                                                                                                                                                                                 >  2:420523/1‑53 (MQ=255)
                                                                                                                                                                                                             tGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAAt                                                                                                                                  >  2:141623/1‑100 (MQ=37)
                                                                                                                                                                                                             tGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAAt                                                                                                                                  <  1:141626/100‑1 (MQ=37)
                                                                                                                                                                                                             tGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACGGTCa                                                                                                      >  1:1125694/1‑128 (MQ=25)
                                                                                                                                                                                                             tGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCg                                        <  2:458037/190‑1 (MQ=255)
                                                                                                                                                                                                             tGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAACGCCACGGGCGGCAACTTCGAACAACTGACCCCGAACCCGAc                                            >  2:1242246/1‑186 (MQ=255)
                                                                                                                                                                                                               cccAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCg                                              <  1:26627/182‑1 (MQ=255)
                                                                                                                                                                                                                  aGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAAccc                                                                                                                                                                  >  1:1262572/1‑63 (MQ=255)
                                                                                                                                                                                                                  aGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCCCGACCGt                                                                                                        >  2:280500/1‑121 (MQ=38)
                                                                                                                                                                                                                     cGCCAGTAACCGGGAACCTCTCCAGCGTCCAAGCGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAg                                                                                                                                              <  2:759222/80‑1 (MQ=14)
                                                                                                                                                                                                                       ccGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGccg                                                                                                                                                                           >  2:825821/1‑47 (MQ=255)
                                                                                                                                                                                                                       ccGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAg                                                                                                                                              >  1:966240/1‑78 (MQ=39)
                                                                                                                                                                                                                       ccGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCCATGACCTTAAAAACAACTCCGCGACCACCAGc                                                                                                    >  1:992665/1‑120 (MQ=11)
                                                                                                                                                                                                                       ccGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAACTGACCCCGAACCCGa                                             >  2:108322/1‑175 (MQ=255)
                                                                                                                                                                                                                        cGGTAACCGTGAACCTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTgg                                                                             <  2:975958/142‑1 (MQ=21)
                                                                                                                                                                                                                            aaCCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCgca                                                                                                                                                                          >  1:702881/1‑43 (MQ=255)
                                                                                                                                                                                                                            aaCCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACAGGCGGCa                                                                        >  2:117040/1‑143 (MQ=37)
                                                                                                                                                                                                                              ccGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCgg                                                                                                                                                                           <  2:1151558/42‑1 (MQ=255)
                                                                                                                                                                                                                                 tGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGcga                                                                                                                                    >  1:518809/1‑76 (MQ=25)
                                                                                                                                                                                                                                 tGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGcga                                                                                                                                    <  2:518804/78‑3 (MQ=255)
                                                                                                                                                                                                                                   aaCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGCCCTTTTCCAGACCCCGGCCACTACCGTCCAGAACAACAGCTCGACTGTCAGCGTCACCATCgaga                                                                                       >  1:861872/1‑121 (MQ=11)
                                                                                                                                                                                                                                    |                                                                                                                                                                                                            
CCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCAC  >  AE015451/195504‑195936

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: