Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 195,424 G→A 12.5% 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,4240GA12.5% 119.7 / 10.0 48V310V (GTG→GTAPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base G (24/18);  new base A (3/3);  total (27/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.17e-01

GTTCCGGCAAACGACACACCTTACATCGACGGTGGCCAGATTTCGGCCACCGTGACAGGTAGTACCGGTGGTGGCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGAC  >  AE015451/195194‑195648
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                
gTTCCGGCAAACGACACACCTTACATCGACGGTGGCCAGATTTCGGCCACCGTGACAGGTAGTACCGGTGGTGGCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGgtgaccg                                                                                                                                                                                                                              <  1:657840/235‑1 (MQ=255)
       cAAACGACACACCTTACATCGACGGTGGCCAGATTTCGGCCACCGTGACAGGTAGTACCGGTGGTACCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGgtaaccgtgac                                                                                                                                                                                                                          >  2:670672/1‑232 (MQ=255)
             acacacCTTACATCGACGGTGGCCAGATTTCGGCCACCGTGACAGGTAGTACCGGTGGTGGCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAAc                                                                                                                                                                                                                <  1:728597/236‑1 (MQ=255)
                 acCTTACATCGACGGTGGCCAGATTTCGGCCACCGTGACAGGTAGTACCGGTGGTGGCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAAc                                                                                                                                                                                                                <  2:986841/232‑1 (MQ=255)
                     tACATCGACGGTGGCCAGATTTCGGCCACCGTGACAGGTAGTACCGGTGGTGGCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGt                                                                                                                                                                                                        <  1:604638/236‑1 (MQ=255)
                                tGGCCAGATTTCGGCCACCGTGACAGGTAGTACCGGTGGTGGCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTg                                                                                                                                                                                             >  1:554684/1‑236 (MQ=255)
                                     aGATTTCGGCCACCGTGACAGGTAGTACCGGTGGTGGCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGAcccc                                                                                                                                                                                                                       >  1:588326/1‑204 (MQ=255)
                                     aGATTTCGGCCACCGTGACAGGTAGTACCGGTGGTGGCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGa                                                                                                                                                                                            <  2:1156689/232‑1 (MQ=255)
                                       aTTTCGGCCACCGTGACAGGTAGTACCGGTGGTGGCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCg                                                                                                                                                                                       <  1:28775/235‑1 (MQ=255)
                                       aTTCCGGCCACCGTGACAGGTAGTACCGGTGGTGGCGGGCTGACGGTGACATTGCCGAAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAAc                                                                                                                                                                                                                <  1:368791/210‑1 (MQ=255)
                                        tttCGGCCACCGTGACAGGTAGTACCGGTGGTGGCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCg                                                                                                                                                                                       <  2:1181852/234‑1 (MQ=255)
                                                                     ggtggCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGAc                                                                                                                                                        >  1:87229/1‑236 (MQ=255)
                                                                       tggCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCa                                                                                                                                                                                                                  >  1:733013/1‑176 (MQ=255)
                                                                             ccTGACGCTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCATGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCaaac                                                                                                                                                                                                        <  2:654089/179‑4 (MQ=18)
                                                                                           tGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCgg                                                                                                                                                   <  2:1213127/219‑1 (MQ=255)
                                                                                            gCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCCACGGCCAgg                                                                                                                                                                                                         >  1:292422/1‑164 (MQ=255)
                                                                                             ccGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACCCTACAACCGCGCCCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGc                                                                                                                                                                                         <  1:405148/179‑1 (MQ=255)
                                                                                                     gCCAGCGGTTCCCCAAGTGCCCGATACAATCGCCACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCt                                                                                                                                                                                                                       <  1:283088/141‑1 (MQ=17)
                                                                                                           ggTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGCGCCCGTGACCCTGTCCCACGGCCAg                                                                                                                                                                                                          >  2:399766/1‑148 (MQ=255)
                                                                                                                  cAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGGAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGgtgaccgtg                                                                                                                                                                                                                            <  2:690081/123‑1 (MQ=255)
                                                                                                                  cAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGt                                                                                                                                                                                                                     >  1:1021142/1‑130 (MQ=255)
                                                                                                                        aCCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCa                                                                                                                  >  1:94603/1‑223 (MQ=11)
                                                                                                                              aCAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGc                                                                                               >  1:261451/1‑236 (MQ=255)
                                                                                                                                  tCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGAc                                                                                                                                           <  1:300291/188‑1 (MQ=255)
                                                                                                                                          aCAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGt                                                                                                                                         >  2:10009/1‑182 (MQ=255)
                                                                                                                                          aCAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGCCATCACCGTTGAAGcc                                                                                                                                                                                        >  2:1031071/1‑135 (MQ=255)
                                                                                                                                          aCAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGCGACCGTGACCCTGCCCAACGGCCAGGTCATCACTGTTGCAGCCGGCAAAACCCAGGGCAGCGTGGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTGCGACCGTCAGc                                                                                                                >  1:423528/1‑207 (MQ=255)
                                                                                                                                            aaCCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGCCATCACCGTTGAAGCGGGCAAAACCCCGGGCac                                                                                                                                                                       >  1:1101423/1‑149 (MQ=255)
                                                                                                                                               cgcgACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGt                                                                                                                    <  2:485633/198‑1 (MQ=255)
                                                                                                                                                  gACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGAc                                                                                                                                           <  1:258814/172‑1 (MQ=255)
                                                                                                                                                      cTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACCATCAcc                                                                                                                                                                                                 >  1:482705/1‑113 (MQ=255)
                                                                                                                                                            gCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCg                                                                                                                                                                                       <  2:758186/118‑1 (MQ=255)
                                                                                                                                                            gCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCg                                                                                                                                                                                       >  1:758191/1‑118 (MQ=255)
                                                                                                                                                                 gCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACg                                                                                                                                                                                                               >  1:1112927/1‑89 (MQ=255)
                                                                                                                                                                                  gTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGAcc                                           <  1:1129752/236‑1 (MQ=255)
                                                                                                                                                                                   tGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCg                                                          <  1:344160/220‑1 (MQ=255)
                                                                                                                                                                                           tCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCGGCTCAGACCACGATCAACg                                 <  1:1078025/236‑1 (MQ=255)
                                                                                                                                                                                             aCCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGc                                                                                               >  1:528451/1‑173 (MQ=33)
                                                                                                                                                                                                       tGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGCCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGc                                                                           >  1:149671/1‑183 (MQ=38)
                                                                                                                                                                                                        gACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGAACCTGTCCAACGGCCAAATGATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCg                        >  1:108612/1‑233 (MQ=255)
                                                                                                                                                                                                          cccTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAAc                                                                                                                                                                                                                >  1:1056095/1‑47 (MQ=255)
                                                                                                                                                                                                          cccTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCATGTCTTAACCATTGATGCGGGCCCACCCCAGGGCAGCGTCGATTGCCAGACCCCGTCCAATTACGt                                                                                                                                         >  1:1051131/1‑118 (MQ=14)
                                                                                                                                                                                                          cccTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGCCCAGAtcctca                                                                                                                                                                                                   >  1:342999/1‑60 (MQ=14)
                                                                                                                                                                                                           ccTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGattcca                  >  1:894033/1‑231 (MQ=255)
                                                                                                                                                                                                            cTGAGCAACCCTGCCCAGACGCCGgtga                                                                                                                                                                                                                                 >  2:545896/1‑28 (MQ=255)
                                                                                                                                                                                                            cTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTg                                                                                                                                                                                                                      >  1:1182145/1‑39 (MQ=255)
                                                                                                                                                                                                              gAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCg                            <  2:894028/223‑1 (MQ=255)
                                                                                                                                                                                                                gCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACgggcg                                                                                       >  1:903096/1‑162 (MQ=21)
                                                                                                                                                                                                                  aaCCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTc                         >  1:1158345/1‑222 (MQ=255)
                                                                                                                                                                                                                   aCCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACtggcggca                                                                                    <  1:564760/162‑3 (MQ=9)
                                                                                                                                                                                                                            aGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCg                                                                                                                                                                                       >  2:834018/1‑54 (MQ=35)
                                                                                                                                                                                                                                 ccGGTGACCGTGACCCTGGCCAACGGCGAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAAt                                                                                                                                              <  2:93525/90‑1 (MQ=255)
                                                                                                                                                                                                                                      gaccgtgaccCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGAc  >  1:1202491/1‑225 (MQ=39)
                                                                                                                                                                                                                                      gaccgtgaccCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGCCTACAACAACGGTGCGACCGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAATTTCGAGCAACTGACCCCGAAccc                                                           >  1:568070/1‑168 (MQ=255)
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                
GTTCCGGCAAACGACACACCTTACATCGACGGTGGCCAGATTTCGGCCACCGTGACAGGTAGTACCGGTGGTGGCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGAC  >  AE015451/195194‑195648

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