Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 213,808 T→A 9.8% A6438A (GCT→GCA)  PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451213,8080TA 9.8% 194.3 / 6.9 60A6438A (GCT→GCA) PP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (31/23);  new base A (3/3);  total (34/26)
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

CAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCC  >  AE015451/213577‑214039
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                       
cAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCg                                                                                                                                                                                                                                     >  1:452194/1‑236 (MQ=255)
            gACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTg                                                                                                                                                                                                                            <  2:513952/233‑1 (MQ=255)
                           cGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTc                                                                                                                                                                                                           <  2:452191/235‑1 (MQ=40)
                                                    aTCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAg                                                                                                                                                                                                                               >  1:62150/1‑190 (MQ=39)
                                                    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)
                                                                                                               gCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCTGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCg                                                                                                                                                                                                                                     >  2:342331/1‑125 (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                                                                                                                                >  2:617283/1‑134 (MQ=255)
                                                                                                                                                                                                           ccAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGc                                                                                                                                <  1:617286/134‑1 (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                                                                                                        >  1:331674/1‑143 (MQ=255)
                                                                                                                                                                                                                          tGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCAcc                                                                                                        <  2:331671/143‑1 (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)
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                       
CAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCC  >  AE015451/213577‑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: ‑