Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 213,820 T→C 8.6% G6442G (GGT→GGC)  PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451213,8200TC 8.6% 175.8 / 2.9 56G6442G (GGT→GGC) PP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (27/23);  new base C (3/2);  total (31/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

GACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCC  >  AE015451/213589‑214039
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                           
gACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTg                                                                                                                                                                                                                            <  2:513952/233‑1 (MQ=255)
               cGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTc                                                                                                                                                                                                           <  2:452191/235‑1 (MQ=40)
                                        aTCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAg                                                                                                                                                                                      >  2:590604/1‑231 (MQ=255)
                                            cGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTc                                                                                                                                                                             >  1:298398/1‑236 (MQ=255)
                                                  ccAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGcgcg                                                                                                                                                                                                                          <  2:224700/185‑1 (MQ=18)
                                                        aaGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTca                                                                                                                                                                    >  1:26280/1‑232 (MQ=255)
                                                        aaGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGt                                                                                                                                                                  <  1:409510/235‑1 (MQ=255)
                                                               ccTGCCGATCAACAACCCCAGTGCAGTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACcgac                                                                                                                                                           <  1:244674/235‑1 (MQ=255)
                                                                    cGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGcc                                                                                                                                                      <  1:85888/235‑1 (MQ=255)
                                                                               cacaGTGCAGTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGACGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAAc                                                                                                                                                                                                    <  2:171983/178‑1 (MQ=255)
                                                                                    tGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCg                                                                                                                                                                                           <  1:690407/182‑1 (MQ=255)
                                                                                    tGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCg                                                                                                                                                                                           >  2:690404/1‑182 (MQ=255)
                                                                                       aTTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGa                                                                                                                                                                                          >  2:765166/1‑180 (MQ=255)
                                                                                       aTTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTAGGCAAGTTCGAgac                                                                                                                                                                                                      >  1:35772/1‑166 (MQ=255)
                                                                                        ttGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTc                                                                                                                                                                                                           <  2:244484/162‑1 (MQ=255)
                                                                                        ttGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTc                                                                                                                                                                                                           >  1:244486/1‑162 (MQ=255)
                                                                                        ttGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGAcc                                                                                                                                                                                   >  2:733792/1‑186 (MQ=255)
                                                                                        ttGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTCATGTCGATCGCAACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGAcc                                                                                                                                                                                   <  1:733795/186‑1 (MQ=255)
                                                                                           aCCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCac                                                                                                                                           >  2:644215/1‑222 (MQ=255)
                                                                                              ttCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACcga                                                                                                                                                            >  1:641010/1‑203 (MQ=255)
                                                                                                 aCGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGac                                                                                                                                                                                                <  2:750359/164‑1 (MQ=255)
                                                                                                 aCGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTCGGTAAGTTCGAGCAACTGac                                                                                                                                                                                                >  1:750362/1‑164 (MQ=255)
                                                                                                   gCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCa                                                                                                                                        >  1:77534/1‑218 (MQ=255)
                                                                                                        gCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACGTCGGTCACCGACGAGCCGGGTACCCCAGGCAACg                                                                                                                                     >  1:201731/1‑216 (MQ=255)
                                                                                                          gACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCa                                                                                                                <  1:550452/235‑1 (MQ=255)
                                                                                                                    ccGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGGGCGGACGTCGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACtt                                                                                                                                                                      >  1:192758/1‑171 (MQ=255)
                                                                                                                            aCCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGa                                                                                              <  1:410151/235‑1 (MQ=255)
                                                                                                                                   cGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGCCCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTgc                                                                                                                                                                                                 >  2:114699/1‑128 (MQ=255)
                                                                                                                                     gCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTgg                                                                                     >  1:786357/1‑235 (MQ=255)
                                                                                                                                                  tGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGc                                                                                    >  1:418742/1‑223 (MQ=255)
                                                                                                                                                         tGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTc                                                                                                                       <  2:298395/181‑1 (MQ=255)
                                                                                                                                                               cGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTc                                                                          <  1:31726/220‑1 (MQ=255)
                                                                                                                                                                        tCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACGTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGACGGCGACCGGGTCAAgg                                                                                                                   >  2:771804/1‑170 (MQ=255)
                                                                                                                                                                          cgcaCAAGGTCTACGTCGGGAGCAACGACCCTGTCGCGATGTCGGTCGGTACCGTGGACGGTGCGGAGGTTGGCAAGTTCGAGCAACTGACACt                                                                                                                                                                                             <  2:412504/91‑1 (MQ=11)
                                                                                                                                                                          cGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAAc                                               >  1:289873/1‑236 (MQ=255)
                                                                                                                                                                           ggACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGt                                                                                                                  >  2:732993/1‑168 (MQ=255)
                                                                                                                                                                           ggACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTc                                                                          <  1:472118/208‑1 (MQ=255)
                                                                                                                                                                                               ccAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGc                                                                                                                                <  1:617286/134‑1 (MQ=255)
                                                                                                                                                                                               ccAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGc                                                                                                                                >  2:617283/1‑134 (MQ=255)
                                                                                                                                                                                                cAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGAcc                                                                                                                                                                                   <  2:544734/82‑1 (MQ=255)
                                                                                                                                                                                                cAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGAcccc                                                                                                                                                                                 <  1:395159/84‑1 (MQ=255)
                                                                                                                                                                                                cAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGAcc                                                                                                                                                                                   >  1:544737/1‑82 (MQ=9)
                                                                                                                                                                                                     aCCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTg                    <  1:407563/236‑1 (MQ=255)
                                                                                                                                                                                                          gTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACTCCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCGCACGACCTGGTCGTGcccc                >  1:680153/1‑235 (MQ=255)
                                                                                                                                                                                                            cGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTa                                                                                                                                                 <  2:798720/104‑1 (MQ=255)
                                                                                                                                                                                                            cGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCt                                       >  2:514998/1‑210 (MQ=255)
                                                                                                                                                                                                            cGTGATGTCGATCGCTACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTa                                                                                                                                                 >  1:798723/1‑104 (MQ=255)
                                                                                                                                                                                                             gTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGAc                                                                  >  1:6446/1‑182 (MQ=255)
                                                                                                                                                                                                              tGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCAcc                                                                                                        <  2:331671/143‑1 (MQ=255)
                                                                                                                                                                                                              tGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCAcc                                                                                                        >  1:331674/1‑143 (MQ=255)
                                                                                                                                                                                                                aTGTCGATCGCTACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGca          <  1:913105/235‑1 (MQ=255)
                                                                                                                                                                                                                   tcgatcgaTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGcaaca       <  1:251259/235‑1 (MQ=255)
                                                                                                                                                                                                                   tCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGcaaca       >  1:629249/1‑235 (MQ=255)
                                                                                                                                                                                                                     gATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCATCTGACACTCGACAAGACCCCGGTCAGCACGTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTAGCCGAGAACGTCaa                                                                        >  2:878597/1‑168 (MQ=255)
                                                                                                                                                                                                                     gATCGCTACCGATGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCt                                                                <  1:788097/176‑1 (MQ=255)
                                                                                                                                                                                                                        cGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGcc  >  1:365052/1‑235 (MQ=255)
                                                                                                                                                                                                                        cGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTGACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCt               >  1:217476/1‑222 (MQ=255)
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                           
GACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCC  >  AE015451/213589‑214039

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