Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 200,113 T→C 13.3% N1873N (AAT→AAC PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451200,1130TC13.3% 229.3 / 18.0 81N1873N (AAT→AACPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (24/46);  new base C (5/6);  total (29/52)
Fisher's exact test for biased strand distribution p-value = 5.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

CCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCG  >  AE015451/199878‑200342
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                     
ccATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAAt                                                                                                                                                                                                                                       <  1:1449847/236‑1 (MQ=38)
                 ggcgCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCa                                                                                                                                                                                                                       <  1:1429614/235‑1 (MQ=37)
                                aaGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCAt                                                                                                                                                                                                                         >  2:948122/1‑218 (MQ=255)
                                  gTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAgcg                                                                                                                                                                                                     <  1:29658/236‑1 (MQ=255)
                                    tcgAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAgcg                                                                                                                                                                                                     <  1:1232042/232‑1 (MQ=39)
                                       aaaaCCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAAcc                                                                                                                                                                                                         <  1:956901/227‑1 (MQ=40)
                                       aaaaCCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAgcgcgc                                                                                                                                                                                                  >  1:1335875/1‑234 (MQ=40)
                                         aacctgaacctgGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAgcg                                                                                                                                                                                                     <  2:966622/229‑1 (MQ=255)
                                                 cctgGACAAGACCCCGGTCAGCACCGCGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCa                                                                                                                                                                                                                                  <  2:725445/192‑1 (MQ=255)
                                                 cctgGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCa                                                                                                                                                                                                                                  >  1:1089972/1‑192 (MQ=21)
                                                 cctgGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCAGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACAccc                                                                                                                                                                                         <  2:1449068/233‑1 (MQ=255)
                                                 cctgGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGCCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTCACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACAcccg                                                                                                                                                                                        >  1:1431136/1‑233 (MQ=255)
                                                 cctgGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCa                                                                                                                                                                                                                                  <  2:1089969/192‑1 (MQ=21)
                                                 cctgGACAAGACCACGGTCAGCACCAAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTAAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCAc                                                                                                                                                                                       >  1:787460/1‑235 (MQ=255)
                                                         aGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACAcc                                                                                                                                                                                          <  1:1248818/224‑1 (MQ=255)
                                                            ccccGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAgg                                                                                                                                                                            <  1:866414/235‑1 (MQ=255)
                                                              ccGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGa                                                                                                                                                                         <  1:1464334/236‑1 (MQ=255)
                                                               cGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGc                                                                                                                                                                           <  2:906309/233‑1 (MQ=255)
                                                                  tCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAAACCAgcg                                                                                                                                                                                                     >  1:940803/1‑204 (MQ=39)
                                                                  tCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAAACCAgcg                                                                                                                                                                                                     <  2:940801/204‑1 (MQ=38)
                                                                        ccGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCt                                                                                                                                                                 >  2:571198/1‑234 (MQ=255)
                                                                         cGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCa                                                                                                                                                                                                                                  <  2:13262/168‑1 (MQ=40)
                                                                         cGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGc                                                                                                                                                                                                                   >  2:1270098/1‑183 (MQ=39)
                                                                         cGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGAt                                                                                                                                                                                    <  2:1446523/214‑1 (MQ=255)
                                                                         cGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATc                                                                                                                                                              >  1:190402/1‑236 (MQ=255)
                                                                         cGAGGACACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCAt                                                                                                                                                                                                                         >  2:296887/1‑177 (MQ=255)
                                                                            ggTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCGAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCAt                                                                                                                                                                                                                         >  2:1373176/1‑174 (MQ=255)
                                                                             gTCACCGACGAGCCAGGTACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGTTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAgg                                                                                                                                                           <  1:929109/235‑1 (MQ=255)
                                                                              tCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGat                                                                                                                                                         >  1:1447307/1‑235 (MQ=25)
                                                                                  cgacgaGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCa                                                                                                                                                                                                        <  2:747105/185‑1 (MQ=255)
                                                                                  cgacgaGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCt                                                                                                                                                      <  2:1063010/235‑1 (MQ=255)
                                                                                     cgaGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGc                                                                                                                                                   >  1:1320093/1‑235 (MQ=255)
                                                                                        gCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCAc                                                                                                                                                                                       >  1:1018676/1‑196 (MQ=255)
                                                                                           aGGCACTCCAGGCAACGAAGGCGACATTGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAGACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATgc                                                                                                                                                                                  >  1:101747/1‑198 (MQ=255)
                                                                                              cACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAgcgcg                                                                                                                                                                                                   <  2:949447/178‑1 (MQ=255)
                                                                                              cACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAgcgcg                                                                                                                                                                                                   >  1:949449/1‑178 (MQ=255)
                                                                                                   cAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGc                                                                                                                                    <  1:795456/236‑1 (MQ=255)
                                                                                                    aGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGc                                                                                                                                    >  1:1119872/1‑235 (MQ=255)
                                                                                                     ggCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGGCTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGa                                                                                                                                                                                     <  2:50186/185‑1 (MQ=14)
                                                                                                     ggCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGa                                                                                                                                                                                     >  1:50186/1‑185 (MQ=25)
                                                                                                      gCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAg                                                                                                                                                                                                       <  2:1496752/166‑1 (MQ=12)
                                                                                                      gCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGc                                                                                                                                                       <  1:1186104/214‑1 (MQ=255)
                                                                                                       cAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCAt                                                                                                                                  <  2:906851/234‑1 (MQ=255)
                                                                                                        aaCGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCa                                                                                                                                <  1:602376/235‑1 (MQ=255)
                                                                                                              ggCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAg                                                                                                                                                >  2:564562/1‑213 (MQ=255)
                                                                                                                cGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAg                                                                                                                       >  1:1006113/1‑236 (MQ=255)
                                                                                                                     tCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGa                                                                                                                   <  1:976557/235‑1 (MQ=255)
                                                                                                                       gTCAAGGTCACCATCACCGCCGACCAGGCTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCAt                                                                                                                                  <  2:685903/218‑1 (MQ=17)
                                                                                                                             gTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTg                                                                                                           <  1:212622/235‑1 (MQ=255)
                                                                                                                                   aTCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCTTTACCCTGAGCAACAATGCTCAGCTcg                                                                                                                                                                                                                             >  2:632472/1‑114 (MQ=255)
                                                                                                                                               cAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAAcc                                                                                         >  1:375040/1‑235 (MQ=255)
                                                                                                                                                              gagaACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCa                                                                                                                        <  2:502347/189‑1 (MQ=255)
                                                                                                                                                                  aCGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCt                                                                                  <  2:1119869/223‑1 (MQ=255)
                                                                                                                                                                    gTCAACCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGtgc                                                                                                       <  2:375038/200‑1 (MQ=255)
                                                                                                                                                                      ctaaCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGGGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAgg                                                                                                                                                           <  2:1343867/144‑1 (MQ=21)
                                                                                                                                                                        aaCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGc                                                                                                                                    <  2:1260355/167‑1 (MQ=255)
                                                                                                                                                                        aaCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGtt                                                               <  1:571200/236‑1 (MQ=255)
                                                                                                                                                                            gACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGc                                                                                                                <  1:1329998/183‑1 (MQ=255)
                                                                                                                                                                                    ccGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACg                                                                                                                 <  2:975156/174‑1 (MQ=255)
                                                                                                                                                                                           cGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCtt                                                                                       <  1:262906/193‑1 (MQ=255)
                                                                                                                                                                                           cGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGc                                                                                                                                    <  2:753170/148‑1 (MQ=37)
                                                                                                                                                                                            atCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCtt                                                                                       <  2:1041248/191‑1 (MQ=255)
                                                                                                                                                                                             tCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCt                                                                                                                                                      <  2:190401/128‑1 (MQ=39)
                                                                                                                                                                                             tCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTcg                                            <  2:595150/234‑1 (MQ=255)
                                                                                                                                                                                               aaCCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACg                                                                                                                 >  2:640707/1‑163 (MQ=255)
                                                                                                                                                                                                aCCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGaa                                       <  2:1263453/236‑1 (MQ=255)
                                                                                                                                                                                                   aGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAgtg                                    <  1:891538/236‑1 (MQ=255)
                                                                                                                                                                                                      cGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCGGCCGTa                                                                                                                     >  2:565682/1‑152 (MQ=255)
                                                                                                                                                                                                       gCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCa                                                                                                                        <  2:62702/148‑1 (MQ=2)
                                                                                                                                                                                                          tgCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACGCCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAgtgg                                   <  1:149620/229‑1 (MQ=255)
                                                                                                                                                                                                           gCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGtt                                                               <  2:1374375/201‑1 (MQ=255)
                                                                                                                                                                                                                  aCCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACt                     >  1:1115761/1‑236 (MQ=255)
                                                                                                                                                                                                                   ccTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGACGTGGTGGCCAAGCTGACTg                    >  1:351137/1‑236 (MQ=255)
                                                                                                                                                                                                                    cTGGTCGTGACCCTGAGCAAACATGCTCAGGTCACCAGCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCg                                                                  <  2:89993/189‑1 (MQ=255)
                                                                                                                                                                                                                        tCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCGTCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTACACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTg                    <  2:798480/231‑1 (MQ=255)
                                                                                                                                                                                                                         cGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGATGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGa                       >  1:896343/1‑227 (MQ=255)
                                                                                                                                                                                                                          gTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAg                           >  2:1084995/1‑222 (MQ=255)
                                                                                                                                                                                                                              cccTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCaa                            <  2:703590/217‑1 (MQ=255)
                                                                                                                                                                                                                                cTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCa                                                                                                                                                                                                        >  2:1064887/1‑43 (MQ=40)
                                                                                                                                                                                                                                acgAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGTCGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCa                             <  2:813316/212‑1 (MQ=255)
                                                                                                                                                                                                                                      aacaaTGCTCAGGTCACCATCAAGGCTGGCGAAACCAgcgcg                                                                                                                                                                                                   <  1:440068/42‑1 (MQ=38)
                                                                                                                                                                                                                                       acaaTGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGTCGCTGGCCAGATCAGCCTGTGCATTAACTCGGCAGTAGACGCCACTGGTGCTCCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGATCAAGCTGACTGCCACCCCTTCGGTGACCg  <  2:787458/234‑1 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                     
CCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCG  >  AE015451/199878‑200342

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