Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 195,424 G→A 10.7% V310V (GTG→GTA PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451195,4240GA10.7% 138.1 / 2.1 54V310V (GTG→GTAPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base G (28/20);  new base A (4/2);  total (32/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCCGGCAAACGACACACCTTACATCGACGGTGGCCAGATTTCGGCCACCGTGACAGGTAGTACCGGTGGTGGCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCG  >  AE015451/195195‑195652
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                    
tcccGGCAAACGACACACCTTACATCGACGGTGGCCAGATTTCGGCCACCGTGACAGGTAGTACCGGTGGTGGCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGgtgaccgtg                                                                                                                                                                                                                                <  2:1136655/234‑1 (MQ=255)
       aaaCGACACACCTTACATCGACGGTGGCCAGATTTCGGCCACCGTGACAGGTAGTACCGGTGGTGGCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGt                                                                                                                                                                                                                         <  1:8235/236‑1 (MQ=255)
                          gACGGTGGCCAGATTTCGGCCACCGTGACAGGTAGTACCGGTGGTGGCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCa                                                                                                                                                                                                               <  2:767198/227‑1 (MQ=255)
                            cGGTGGCCAGATTTCGGCCACCGTGACAGGTAGTACCGGTGGTGTCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGgtgaccgtg                                                                                                                                                                                                                                <  1:496452/208‑1 (MQ=255)
                               tGGCCAGATTTCGGCCACCGTGACAGGTAGTACCGGTGGTGGCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGgtaaccg                                                                                                                                                                                                                                  >  1:1124824/1‑203 (MQ=255)
                               tGGCCAGATTTCGGCCACCGTGACAGGTAGTACCGGTGGTGGCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGACACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGtt                                                                                                                                                                                                  >  1:884438/1‑235 (MQ=255)
                                      aTTTCGGCCACCGTGACAGGTAGTACCGGTGGTGGCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCg                                                                                                                                                                                           <  1:1165164/235‑1 (MQ=255)
                                        ttCGGCCACCGTGACAGGTAGTACCGGTGGTGGCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGa                                                                                                                                                                                                <  2:838634/228‑1 (MQ=255)
                                                   tGACAGGTAGTACCGGTGGTGGCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGcc                                                                                                                                                                                                                >  2:881544/1‑201 (MQ=255)
                                                              aCCGGTGGTGGCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAAc                                                                                                                                                                                                                    <  1:786443/186‑1 (MQ=255)
                                                              aCCGGTGGTGGCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAAc                                                                                                                                                                                                                    >  2:786440/1‑186 (MQ=255)
                                                                cGGTGGTGGCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGACCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTc                                                                                                                                                                 >  1:449875/1‑235 (MQ=255)
                                                                    ggtggCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGgtcaccg                                                                                                                                                                                                                                  >  1:491553/1‑166 (MQ=255)
                                                                      tggCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAg                                                                                                                                                                                >  1:610666/1‑214 (MQ=255)
                                                                      tggCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAg                                                                                                                                                                                <  2:610663/214‑1 (MQ=255)
                                                                      tggCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGAccc                                                                                                                                                          >  1:624703/1‑236 (MQ=255)
                                                                          gggCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCCGACGCCGgtgaccgt                                                                                                                                                                                                                                 >  1:136961/1‑161 (MQ=255)
                                                                          gggCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAAc                                                                                                                                                                                                                    <  2:815042/174‑1 (MQ=255)
                                                                              tGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCa                                                                                                                                                    <  2:507647/234‑1 (MQ=255)
                                                                                           gCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGcc                                                                                                                                                     <  2:884435/220‑1 (MQ=255)
                                                                                             cGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCa                                                                                                                                                    <  2:1108043/219‑1 (MQ=255)
                                                                                                   cGCCAGCGGTTACCCAAGTGCCCGATACCATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACg                                                                                                                                                                                                                   <  1:238098/150‑1 (MQ=255)
                                                                                                             tccccAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGt                                                                                                                                                                                                            <  1:212856/145‑1 (MQ=255)
                                                                                                                 cAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGcc                                                                                                                                                                                                                >  2:832628/1‑139 (MQ=255)
                                                                                                                      gACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAAt                                                                                                                                                  <  1:881547/196‑1 (MQ=11)
                                                                                                                       aCCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCg                                                                                                          >  1:691634/1‑235 (MQ=255)
                                                                                                                             aCAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGCCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGc                                                                                                   >  1:431879/1‑236 (MQ=255)
                                                                                                                             aCAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGACCAgg                                                                                                                                                                                                             >  1:1217989/1‑130 (MQ=255)
                                                                                                                                aTCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGtt                                                                                                                               <  2:806792/205‑1 (MQ=11)
                                                                                                                                    acacTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGATCATCACCGTTGAAGCCGGCAAAAcc                                                                                                                                                                                   >  1:1260888/1‑149 (MQ=255)
                                                                                                                                                 gACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAg                                                                                                                     <  2:866193/198‑1 (MQ=255)
                                                                                                                                                   cccTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCa                                                                                                                                                                                                               >  1:347596/1‑106 (MQ=255)
                                                                                                                                                   cccTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCATGCCCAGACGCCGGTGACCGCGACGCTGCCCAACGGc                                                                                                                                                                                                                 >  1:1192181/1‑104 (MQ=14)
                                                                                                                                                     cTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGAccc                                                                                                                                                          >  1:804357/1‑157 (MQ=255)
                                                                                                                                                         ccGCCTCGCCAAGCGTGGCCGAAGGTGGCGTGATCACCTTCACCGTGACCCTGAGCAACCCGGCCCAGACGCCGGTGACCGCGACCCTGTCCAAc                                                                                                                                                                                                                    <  2:301156/95‑1 (MQ=255)
                                                                                                                                                               cGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGaa                                                                  <  2:610834/235‑1 (MQ=255)
                                                                                                                                                                  cAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCa                                                                                                                                                    >  1:856523/1‑150 (MQ=255)
                                                                                                                                                                  cAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCa                                                                                                                                                    <  2:856520/150‑1 (MQ=255)
                                                                                                                                                                          aCCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAAc                                                                                                                                                                                                                    >  1:1237930/1‑78 (MQ=255)
                                                                                                                                                                          aCCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAAc                                                                                                                                                                                                                    <  2:1237927/78‑1 (MQ=255)
                                                                                                                                                                               aGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGa                                                 <  1:463744/236‑1 (MQ=255)
                                                                                                                                                                                       tGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAAc                                                                                                                                                                                                                    <  2:232676/65‑1 (MQ=25)
                                                                                                                                                                                        gATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGCCgt                                                                                                                                                                     >  1:286580/1‑110 (MQ=37)
                                                                                                                                                                                            aCCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGCCCGTCAGCgc                                                                                                                  >  1:56476/1‑157 (MQ=18)
                                                                                                                                                                                                       gACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGa                                    >  2:894643/1‑225 (MQ=255)
                                                                                                                                                                                                         cccTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCg                                                                    >  1:515506/1‑191 (MQ=16)
                                                                                                                                                                                                           cTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGACCGCCACcggcggca                                                                                        >  1:126274/1‑167 (MQ=17)
                                                                                                                                                                                                           cTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACcggcggc                                                                                         >  2:472240/1‑167 (MQ=18)
                                                                                                                                                                                                           cTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCAACACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCAcc                                                                                               >  1:342096/1‑162 (MQ=18)
                                                                                                                                                                                                           cTGAGCAACCCTGCCCAGACGCCGGTGACCGTGAACCTGTCCAACGGCCAAGCCATCACCGTTGAAGCTGGCAACTGCACCGACAGCCTCGATCTCCAGACTCCCGCGAATGACGtatg                                                                                                                                          >  1:1202705/1‑116 (MQ=2)
                                                                                                                                                                                                                 aaCCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACgg                                                                                                                                                                                                                  >  2:706188/1‑41 (MQ=255)
                                                                                                                                                                                                                 aaCCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCgaga                                                                                                       >  1:946531/1‑148 (MQ=34)
                                                                                                                                                                                                                 aaCCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAGCGCCACcggcggca                                                                                        >  1:459338/1‑161 (MQ=17)
                                                                                                                                                                                                                 aaCCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAACTTCg                                                                                  >  1:401937/1‑169 (MQ=18)
                                                                                                                                                                                                                 aaCCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGCTCATCACCGTTGAACCCGGCAACACCCAg                                                                                                                                                                                >  1:462517/1‑75 (MQ=255)
                                                                                                                                                                                                                      tGCCCAGACGCCGGTGACCGTGACCCTGTCCTACg                                                                                                                                                                                                                   >  1:447468/1‑35 (MQ=255)
                                                                                                                                                                                                                           aGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCGGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGGCTACAACAACGGTTCGACCGCCAGCGTCACCATCGAGAAAGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGcc                                                     >  1:25499/1‑187 (MQ=255)
                                                                                                                                                                                                                           aGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGt                                                                                                                        <  2:1135567/121‑1 (MQ=33)
                                                                                                                                                                                                                           aGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGt                                                                                                                        >  1:1135570/1‑121 (MQ=33)
                                                                                                                                                                                                                                     gaccgtgaccCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCg  >  1:578142/1‑229 (MQ=255)
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                    
TTCCGGCAAACGACACACCTTACATCGACGGTGGCCAGATTTCGGCCACCGTGACAGGTAGTACCGGTGGTGGCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCG  >  AE015451/195195‑195652

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: