Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 210,319 C→G 7.1% T5275T (ACC→ACG PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451210,3190CG 7.1% 236.6 / 5.7 69T5275T (ACC→ACGPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base C (28/36);  new base G (3/2);  total (31/38)
Fisher's exact test for biased strand distribution p-value = 6.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

CGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACG  >  AE015451/210085‑210547
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                    
cGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCAccc                                                                                                                                                                                                                                     <  1:249879/236‑1 (MQ=255)
              acagcacCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAg                                                                                                                                                                                                                         <  1:1347635/230‑1 (MQ=255)
                   acCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCt                                                                                                                                                                                                                  <  1:891135/236‑1 (MQ=255)
                     ccACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGc                                                                                                                                                                                                                <  1:1312487/236‑1 (MQ=255)
                                 gggTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAAcc                                                                                                                                                                                                    <  1:283361/236‑1 (MQ=255)
                                         aCGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGc                                                                                                                                                                                             <  1:686489/235‑1 (MQ=255)
                                            tttACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCt                                                                                                                                                                                                                                    <  2:896536/193‑1 (MQ=255)
                                                 aaCGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCtt                                                                                                                                                                                    >  1:100681/1‑236 (MQ=255)
                                                 aaCGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCtt                                                                                                                                                                                    >  1:1372784/1‑236 (MQ=255)
                                                           gCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACACAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAg                                                                                                                                                                           <  1:840500/235‑1 (MQ=255)
                                                             cAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAg                                                                                                                                                                           <  1:639775/233‑1 (MQ=255)
                                                              aGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGa                                                                                                                                                                        >  1:1169242/1‑235 (MQ=255)
                                                                  cAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACTCCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATc                                                                                                                                                                                                           <  1:1441189/196‑1 (MQ=255)
                                                                     ccTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCaacaac                                                                                                                                                                                                     <  2:1025299/199‑1 (MQ=255)
                                                                        gggCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGAcc                                                                                                                                                              >  1:850247/1‑235 (MQ=255)
                                                                        gggCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCa                                                                                                                                                                                  >  2:889801/1‑215 (MQ=255)
                                                                               aCCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCAc                                                                                                                                                                                 <  2:1033676/209‑1 (MQ=255)
                                                                                 cTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCg                                                                                                                                                    <  1:805394/236‑1 (MQ=255)
                                                                                  tCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCAt                                                                                                                                                                                           >  2:1134609/1‑196 (MQ=255)
                                                                                     gcggTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCCGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTTGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCaacaac                                                                                                                                                                                                     <  1:807048/183‑1 (MQ=255)
                                                                                           gACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGtcat                                                                                                                                                         >  2:573856/1‑221 (MQ=255)
                                                                                            aCGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCAACACCCTGACcaaca                                                                                                                                                                                                                            >  2:435648/1‑153 (MQ=255)
                                                                                                    ggTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTTAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGc                                                                                                                                                                                             <  2:800830/176‑1 (MQ=255)
                                                                                                        cGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGc                                                                                                                                                                                                                <  1:952080/153‑1 (MQ=11)
                                                                                                        cGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGc                                                                                                                                                                                                                >  2:952078/1‑153 (MQ=14)
                                                                                                           ccTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACTGcc                                                                                                                                                                                                                   >  1:920274/1‑147 (MQ=255)
                                                                                                           ccTTCGAGAACCTGGAGCTGGGCGGCGAGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGAc                                                                                                                                                                                                                       <  2:866656/143‑1 (MQ=12)
                                                                                                             ttCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGtca                                                                                                                                                          <  2:947497/202‑1 (MQ=255)
                                                                                                                  gaACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTg                                                                                                                                                  >  2:600273/1‑205 (MQ=255)
                                                                                                                  gaACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTg                                                                                                                                                  <  1:600275/205‑1 (MQ=255)
                                                                                                                  gaACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCg                                                                                                                   <  1:791553/236‑1 (MQ=255)
                                                                                                                  gaACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATGGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAAcc                                                                                                                    >  1:835898/1‑235 (MQ=255)
                                                                                                                   aaCCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCg                                                                                                                   <  1:1078919/235‑1 (MQ=255)
                                                                                                                   aaCCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCg                                                                                                                   <  1:1504327/235‑1 (MQ=255)
                                                                                                                    aCCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCg                                                                                                                   <  1:1122531/234‑1 (MQ=255)
                                                                                                                        ggAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACcaaca                                                                                                                                                                                                                            >  2:1422182/1‑125 (MQ=255)
                                                                                                                            cTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGcccc                                                                                                         <  1:573858/236‑1 (MQ=255)
                                                                                                                              gggCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCg                                                                                                        <  1:407167/235‑1 (MQ=255)
                                                                                                                                       ggCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGtcat                                                                                                                                                         >  2:700139/1‑177 (MQ=38)
                                                                                                                                          ttCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAAc                                                                                                                                         <  1:259634/190‑1 (MQ=37)
                                                                                                                                           tCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCg                                                                                                        <  1:1238175/222‑1 (MQ=40)
                                                                                                                                                 cAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCAATCAACAAcc                                                                                                                                                                                                    >  1:1204511/1‑124 (MQ=17)
                                                                                                                                                       aCCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCAcc                                                                                                                                                                                                                                      <  1:644050/84‑1 (MQ=31)
                                                                                                                                                        ccGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGt                                                                                                                                 <  1:1261297/184‑1 (MQ=38)
                                                                                                                                                        ccGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACg                                                                               <  2:1456990/234‑1 (MQ=255)
                                                                                                                                                            caccacCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGAc                                                                                                               >  2:1177244/1‑198 (MQ=40)
                                                                                                                                                                       aaGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCg                                                              >  1:377405/1‑236 (MQ=40)
                                                                                                                                                                         gtggtgGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCg                                                                                                                   >  2:1374778/1‑181 (MQ=14)
                                                                                                                                                                            gtggCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGCCCGTCATCACCGTGCCGGCCCACgg                                                                                                                                       >  2:32730/1‑158 (MQ=255)
                                                                                                                                                                                 cAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACAGCAAGACCGTCACCACCGTGCCGGCCaa                                                                                                                                          >  1:1268666/1‑150 (MQ=255)
                                                                                                                                                                                      tcaCTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGTCGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACt                                                                                                                          <  1:28611/159‑1 (MQ=255)
                                                                                                                                                                                       gACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCg                                                    >  1:1086639/1‑230 (MQ=255)
                                                                                                                                                                                       gACTGCGACCCCGCCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGgcagca                                                                                  >  2:811826/1‑200 (MQ=255)
                                                                                                                                                                                            cGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCCCCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCCGCAACGACCCCGt                                                                        >  2:112382/1‑205 (MQ=255)
                                                                                                                                                                                             gACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTc                                                                       <  2:1201258/205‑1 (MQ=255)
                                                                                                                                                                                                 ccGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGGGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCgct                                                            >  1:281282/1‑210 (MQ=255)
                                                                                                                                                                                                     cGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTc                                                               >  1:848244/1‑205 (MQ=255)
                                                                                                                                                                                                     cGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTc                                                               <  2:848242/205‑1 (MQ=255)
                                                                                                                                                                                                       gTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGACGTCGGTaa                              >  1:569262/1‑236 (MQ=255)
                                                                                                                                                                                                         cACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTg                                           >  1:1169847/1‑221 (MQ=255)
                                                                                                                                                                                                            cGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGACGTCGGTAAGTTc                          >  1:1206570/1‑235 (MQ=255)
                                                                                                                                                                                                              agggcggcGAGATCACCTCCCCCAGCCCCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGt                                                                                                                                 <  2:131547/128‑1 (MQ=11)
                                                                                                                                                                                                              aaGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAAcc                                                                                                                                                                                                    <  2:1287048/63‑1 (MQ=40)
                                                                                                                                                                                                                   ggcgAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACgg                                                                                                                                       <  2:1158244/119‑1 (MQ=39)
                                                                                                                                                                                                                   ggcgAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACgg                                                                                                                                       >  1:1158247/1‑119 (MQ=39)
                                                                                                                                                                                                                     tgagaTCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGc                                                                                                                        <  2:289962/131‑1 (MQ=12)
                                                                                                                                                                                                                       agaTCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGACGTCGGTAAGTTCGAGCAACTg                 <  2:1372781/233‑1 (MQ=255)
                                                                                                                                                                                                                          tCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGACGTCGGTAAGTTCGAGCAACTGACg              <  2:1206567/233‑1 (MQ=255)
                                                                                                                                                                                                                                acacCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACACCCACAGTGCAGTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTc                          <  1:114333/215‑1 (MQ=255)
                                                                                                                                                                                                                                acacCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCCACCGTTGGCGGTGc                                          >  2:181620/1‑199 (MQ=255)
                                                                                                                                                                                                                                   ccATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGACGTCGGTAAGGTCGAGCAACTGACGCTGGACAAGACg  >  1:325950/1‑236 (MQ=255)
                                                                                                                                                                                                                                          ccTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGGCATCACCGAGCCCGCCAACGGCACCGTGGGGCCTGCAACAGTGACTGCCCCGGACa                                                                                                    >  2:1235930/1‑131 (MQ=11)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                    
CGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACG  >  AE015451/210085‑210547

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