Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 201,202 C→G 6.1% A2236A (GCC→GCG PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451201,2020CG 6.1% 223.8 / 3.6 66A2236A (GCC→GCGPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base C (32/29);  new base G (2/2);  total (35/31)
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.92e-01

CTGACCGCGACCCCTTCAGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGC  >  AE015451/200975‑201429
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                   
cTGACCGCGACCCCTTCAGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAg                                                                                                                                                                                                                                >  1:379210/1‑233 (MQ=255)
 tGACCGCGACCCCTTCAGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAgcg                                                                                                                                                                                                                              <  2:666962/234‑1 (MQ=255)
                 aGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGt                                                                                                                                                                                                            <  1:723816/236‑1 (MQ=255)
                  gTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGt                                                                                                                                                                                                            <  1:283857/235‑1 (MQ=255)
                   tCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGtt                                                                                                                                                                                                           <  1:702228/235‑1 (MQ=255)
                                tGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAaacctgaac                                                                                                                                                                                              >  1:775505/1‑235 (MQ=255)
                                                aTCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGTGCTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGCACCTGGACAAGACCCCGGt                                                                                                                                                                             >  1:102170/1‑236 (MQ=255)
                                                 tCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACTAATGCGCCGGTGGTCAATGCCATCGACGCCGTc                                                                                                                                                                                                                                  >  2:737566/1‑182 (MQ=255)
                                                    cccTGACCAACAAAGATGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAACACCTGAACCTGGACAAGACCCCGGACAg                                                                                                                                                                          >  1:119993/1‑235 (MQ=255)
                                                    cccTGACCAACAAAGATGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAg                                                                                                                                                                          >  1:799463/1‑235 (MQ=40)
                                                    cccTGACCAACAAAGATGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGc                                                                                                                                                                         <  1:293608/236‑1 (MQ=40)
                                                                      gTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAg                                                                                                                                                       >  1:523097/1‑236 (MQ=255)
                                                                      gTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACTAATGCGCCGGTGGGCAATGCCATCGACGCAGt                                                                                                                                                                                                                                   >  1:510403/1‑160 (MQ=255)
                                                                        cTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGcc                                                                                                                                                     >  1:284670/1‑236 (MQ=255)
                                                                           ccGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAg                                                                                                                                                   <  1:422148/235‑1 (MQ=255)
                                                                              aTCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCCACTGCAACCGCCCCAGACAACGTCTATGTCGGCACCAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAaacct                                                                                                                                                                                                  >  2:782689/1‑185 (MQ=25)
                                                                                 aacaacCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACt                                                                                                                                              >  1:513785/1‑234 (MQ=255)
                                                                                            ggAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGa                                                                                                                                   <  1:193839/234‑1 (MQ=255)
                                                                                                    aCTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACa                                                                                                                           >  2:591202/1‑234 (MQ=255)
                                                                                                       tCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCg                                                                                                                        <  2:755596/234‑1 (MQ=255)
                                                                                                             ggACCGACTGCACCACTGTCGTCGTGGCAGCCAACAGCCCCACGGGCTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGc                                                                                                                                                      <  1:226080/197‑1 (MQ=255)
                                                                                                                    cGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGc                                                                                                                                        <  2:102856/205‑1 (MQ=255)
                                                                                                                           acTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACTAATGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGaa                                                                                                                                  >  1:178616/1‑204 (MQ=255)
                                                                                                                                 gtcgtGGCAGCCAACAGCACCACGGGCTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCa                                                                                                            <  1:591580/220‑1 (MQ=255)
                                                                                                                                       gCAGCCAACAGCACCACGGGCTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACTAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAaacctgaacct                                                                                                                                                                                            <  1:553145/134‑1 (MQ=17)
                                                                                                                                       gCAGCCAACAGCACCACGGGCTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAaacctgaacct                                                                                                                                                                                            >  2:553143/1‑134 (MQ=37)
                                                                                                                                         aGCCAACAGCACCACGGGCTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGAGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTACAGGCAACGAAgg                                                                                                                                >  2:746057/1‑192 (MQ=255)
                                                                                                                                            cAACAGCACCACGGGCTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGc                                                                                                                                                 >  1:578989/1‑172 (MQ=255)
                                                                                                                                            cAACAGCACCACGGGCTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGc                                                                                                                                                 <  2:578987/172‑1 (MQ=255)
                                                                                                                                            cAACAGCACCACGGGCTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGAACTCCAGGCAACGAAGGCGACATCGTCAAGGt                                                                                                                 >  2:674245/1‑204 (MQ=255)
                                                                                                                                             aaCAGCACCACGGGCTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAg                                                                                                                   <  1:259351/201‑1 (MQ=255)
                                                                                                                                                        gggCTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAAc                                                                     >  1:206968/1‑236 (MQ=255)
                                                                                                                                                        gggCTCGGCCACGGCAACCGCACCAGACAACGTCGATGTCGGCACTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAgcggcg                                                                                                                                                                                                                           <  1:188802/86‑1 (MQ=255)
                                                                                                                                                          gCTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGCCGAGCCCGGCACTCCa                                                                                                                                           >  2:239869/1‑164 (MQ=255)
                                                                                                                                                          gCTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCaaa                                                                      >  2:307325/1‑233 (MQ=255)
                                                                                                                                                           cTCGACCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGc                                                                                                                                                      <  2:284611/152‑1 (MQ=255)
                                                                                                                                                            tCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGAc                                                                 <  1:16511/236‑1 (MQ=255)
                                                                                                                                                             cGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCa                                                                        >  2:561514/1‑228 (MQ=255)
                                                                                                                                                              ggCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCt                                                               <  1:307327/236‑1 (MQ=255)
                                                                                                                                                               gCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACTAATGCGCCGGTGGTCAATGCCATCGACGCGGTCAgcggcg                                                                                                                                                                                                                           >  2:719679/1‑79 (MQ=21)
                                                                                                                                                               gCCACTGCAACCACACCAGACAACGTCTATGTCGGCACTAATGCGCCGGTGGTCAATGCCATCGACGCGGTCAgcggcg                                                                                                                                                                                                                           <  1:719681/79‑1 (MQ=255)
                                                                                                                                                                ccACTGCAACCGCACCAGACAACGTCTATGTCGGCACTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCAc                                                                                                                                                              >  1:25688/1‑139 (MQ=255)
                                                                                                                                                                 cACTGCAACCGCACCAGACAACGTCTATGTCGGCACTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGcc                                                                                                    <  2:751676/196‑1 (MQ=255)
                                                                                                                                                                  aCTGCAACCGCACCAGACAACGTCTATGTCGGCACTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCg                                                                                                                              >  2:362589/1‑169 (MQ=255)
                                                                                                                                                                  aCTGCAACCGCACCAGACAACGTCTATGTCGGCACTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCt                                                               <  2:513783/232‑1 (MQ=255)
                                                                                                                                                                   cTGCAACCGCACCAGACAACGTCTATGTCGGCACTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACg                                                                           <  2:799461/219‑1 (MQ=255)
                                                                                                                                                                      cAACCGCACCAGACAACGTCTATGTCGGCACTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGAc                                                                 <  2:379208/226‑1 (MQ=255)
                                                                                                                                                                               cAGACAACGTCTATGTCGGCACTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAAcc                                                                    <  2:206966/214‑1 (MQ=39)
                                                                                                                                                                                aGACAACGTCTATGTCGGCACTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCaa                                               <  2:775503/234‑1 (MQ=255)
                                                                                                                                                                                                ggCACTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGa                             <  1:126544/236‑1 (MQ=255)
                                                                                                                                                                                                   aCTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAgg                                                                                                                                <  2:327195/134‑1 (MQ=37)
                                                                                                                                                                                                   aCTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAgg                                                                                                                                >  1:327197/1‑134 (MQ=38)
                                                                                                                                                                                                    cTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCCACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGCCCTCGGTAGCCgag                                                                               >  1:263/1‑182 (MQ=37)
                                                                                                                                                                                                    cTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCa                                                                        >  1:326777/1‑189 (MQ=37)
                                                                                                                                                                                                       aTGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAACACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGAt                       <  1:444779/235‑1 (MQ=255)
                                                                                                                                                                                                          cgcCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGAcccc                                                                                                                                                                                >  2:748580/1‑79 (MQ=34)
                                                                                                                                                                                                          cgcCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTc                                                             <  1:800705/194‑1 (MQ=37)
                                                                                                                                                                                                           gcCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCCACGAAGGCGACATCGTCAAGGTCACCATCACGGCCTACCACACCTCGGTAccc                                                                                  >  2:57018/1‑172 (MQ=255)
                                                                                                                                                                                                            ccGGTGGTCAATGCCATGGACGCCGGCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGAt                       <  2:620990/230‑1 (MQ=255)
                                                                                                                                                                                                               gtggtCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACgg                                                                                                      >  1:175271/1‑148 (MQ=32)
                                                                                                                                                                                                               gtggtCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGCCACCAt                                                                                                           >  1:190336/1‑143 (MQ=17)
                                                                                                                                                                                                                  gtCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACtgtg                                                       >  1:76543/1‑192 (MQ=38)
                                                                                                                                                                                                                          cATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACg                                                                                                                                                         >  1:352928/1‑86 (MQ=33)
                                                                                                                                                                                                                          cATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACg   >  1:527195/1‑236 (MQ=255)
                                                                                                                                                                                                                          cATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGGAGCCGAGAACGTCAAACCGACCTTCACTGTGCACg                                                   >  1:172238/1‑188 (MQ=39)
                                                                                                                                                                                                                           aTCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGc  <  1:475648/236‑1 (MQ=255)
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                   
CTGACCGCGACCCCTTCAGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGC  >  AE015451/200975‑201429

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