Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 2,550,125 G→A 48.1% noncoding (1436/1522 nt) PP_16SF → 16S ribosomal RNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154512,550,1250GA48.1% 49.4 / 122.4 79noncoding (1436/1522 nt)PP_16SF16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base G (21/20);  new base A (18/20);  total (39/40)
Fisher's exact test for biased strand distribution p-value = 8.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCGGGCCTTGTACACACCGCCCGTCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGCAGTTCAAATCT  >  AE015451/2550051‑2550359
                                                                          |                                                                                                                                                                                                                                          
cccGGGCCTTGTACACACCGCCCGTCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAAcgcg                                                                           <  1:109670/236‑1 (MQ=11)
 ccGGGCCTTGTACCCACCGCCCGTCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGa                                                                          <  1:277425/236‑1 (MQ=11)
 ccGGGCCTTGTACACACCGCCCGTCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGa                                                                          <  1:1173579/236‑1 (MQ=11)
 ccGGGCCTTGTACACACCGCCCGTCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGa                                                                          <  1:1132932/236‑1 (MQ=11)
 ccGGGCCTTGTACACACCGCCCGTCACACCATGGGAGTGGGGTGCACCAGAAGTAGCTAGTCTAACCTTCGGGAGGACGGTTACCACGGAGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTACCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGa                                                                          >  1:242725/1‑236 (MQ=2)
   gggCCTTGTACACACCGCCCGTCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGAcc                                                                        >  2:287936/1‑236 (MQ=31)
   gggCCTTGTACACACCGCCCGTCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGAcc                                                                        >  1:941187/1‑236 (MQ=14)
    ggCCTTGTACACACCGCCCGTCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGAccc                                                                       <  1:548392/236‑1 (MQ=31)
     gCCTTGTACACACCGCCCGTCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGAccc                                                                       >  1:870397/1‑235 (MQ=31)
     gCCTTGTACACACCGCCCGTCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCt                                                                      >  1:144871/1‑236 (MQ=14)
         tGTACACACCGCCCGTCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGa                                                                          <  1:591061/228‑1 (MQ=11)
            acacacCGCCCGTCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTtata                                                                >  1:1128348/1‑235 (MQ=31)
            acacacCGCCCGTCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAg                                                               >  1:134367/1‑236 (MQ=31)
             cacacCGCCCGTCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTtata                                                                >  2:775590/1‑234 (MQ=14)
              acacCGCCCGTCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAAcgcg                                                                           >  2:294973/1‑222 (MQ=11)
                 ccGCCCGTCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTg                                                          <  2:258895/236‑1 (MQ=14)
                  cGCCCGTCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTg                                                          <  1:429663/235‑1 (MQ=31)
                   gCCCGTCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTAGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTa                                                        <  1:784580/236‑1 (MQ=14)
                     ccGTCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGc                                                      <  1:356697/236‑1 (MQ=31)
                     ccGTCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGc                                                      >  1:1206704/1‑236 (MQ=14)
                       gTCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGAcc                                                                        >  1:6077/1‑216 (MQ=31)
                       gTCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTc                                                            >  1:1025417/1‑228 (MQ=31)
                       gTCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTa                                                        >  2:1074632/1‑232 (MQ=31)
                       gTCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTa                                                        <  1:1074637/232‑1 (MQ=31)
                       gTCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGc                                                      >  1:627214/1‑234 (MQ=31)
                       gTCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTc                                                    >  1:27051/1‑236 (MQ=31)
                       gTCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCAAAGCCGATGCTGTAACGCGACCCTg                                                                     >  2:536276/1‑219 (MQ=14)
                       gTCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGa                                                                          <  1:1155724/214‑1 (MQ=11)
                       gTCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGc                                                      >  2:859753/1‑234 (MQ=14)
                       gTCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCt                                                     >  1:708483/1‑235 (MQ=14)
                       gTCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTc                                                    >  1:168167/1‑236 (MQ=14)
                       gTCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCATGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCt                                                     >  1:192034/1‑235 (MQ=14)
                          acacCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGt                                                 >  1:827013/1‑236 (MQ=31)
                              cATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTgg                                              <  1:1274951/235‑1 (MQ=31)
                              cATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTgg                                              <  1:465666/235‑1 (MQ=14)
                               aTGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAg                                                               >  2:1099593/1‑217 (MQ=14)
                                  ggAGTGGGTTGCACCAGAAGTAGCTAGTCTATCCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTaga                                         <  1:962404/236‑1 (MQ=14)
                                  ggAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTaga                                         <  1:294977/236‑1 (MQ=31)
                                  ggAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGCAGCTCAGTTGGTTaga                                         <  2:144868/236‑1 (MQ=14)
                                   gAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTagag                                        >  1:1083943/1‑236 (MQ=31)
                                   gAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTaga                                         <  1:601684/235‑1 (MQ=14)
                                    aGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTaga                                         <  2:1044732/234‑1 (MQ=31)
                                    aGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAgc                                       <  1:471175/236‑1 (MQ=14)
                                     gTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTaga                                         <  2:870392/233‑1 (MQ=31)
                                     gTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTaga                                         >  2:1257762/1‑233 (MQ=14)
                                     gTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTaga                                         <  1:1257767/233‑1 (MQ=14)
                                      tGGGTTGCTCCAGAAGTAGCTAGTCTAACCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAgcgc                                     <  1:347702/236‑1 (MQ=14)
                                      tGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGAGGACGGTTACCACGGTGCGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAgcgc                                     <  1:1218357/236‑1 (MQ=14)
                                      tGGGTTGCACCAGAAGTAGATAGTCTAACCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAgcgc                                     <  1:295205/236‑1 (MQ=14)
                                         gTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTaga                                         <  1:854777/229‑1 (MQ=31)
                                         gTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTg                                                                     >  2:351259/1‑201 (MQ=14)
                                           tGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAgc                                       >  2:1016747/1‑229 (MQ=31)
                                           tGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAGCGCAccc                                 >  1:544210/1‑235 (MQ=31)
                                           tGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAGCGCAcc                                  >  2:897199/1‑234 (MQ=31)
                                           tGCACCAGAAGTAGCTAGTCTAACCTTCGGGGCGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCAGTAGGGCAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCATGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAg                                                       >  2:1213304/1‑213 (MQ=9)
                                           tGCACCAGAAGTAGCTAGTCTAACCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAgc                                       >  2:1225715/1‑229 (MQ=14)
                                           tGCACCAGAAGTAGCTAGTCTAACCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAGCGCAccc                                 >  1:1071482/1‑235 (MQ=14)
                                             cACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAGCGCACCCCTg                              <  2:889487/236‑1 (MQ=31)
                                               ccAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAGCGCAcccc                                <  1:758215/232‑1 (MQ=31)
                                                cAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGc                                                      >  2:51330/1‑209 (MQ=31)
                                                 aGAAGTAGCTAGTCTAACCGTCGGGGGGACGGTTACCACGGTGTTATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGTGACCCTGTTATAGGTCTGTa                                                        <  2:390391/206‑1 (MQ=14)
                                                     gTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAGCGCACCCCTGATAAggg                       >  1:1233773/1‑235 (MQ=31)
                                                     gTAGCTAGGCGAACCGTCGGGAGGACGGTTACCACGGTGTGATTCAGGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTaga                                         <  2:57582/217‑1 (MQ=9)
                                                      tAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAGCGCACCCCTGATAAgg                        >  2:1102732/1‑233 (MQ=31)
                                                         cTAGTCTAACCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTgg                                              <  2:304241/208‑1 (MQ=14)
                                                         cTAGTCTAAACTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAgc                                       <  1:151149/215‑1 (MQ=14)
                                                          tAGTCTAACCTACGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAGCGCACCCCTGATAAGGGTGAGGt                 <  1:270732/236‑1 (MQ=14)
                                                           aGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAGCGCAcccc                                <  1:1276447/220‑1 (MQ=31)
                                                           aGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAGCGCACCCCTGATAAGGGTGAgg                  >  2:1126159/1‑234 (MQ=31)
                                                           aGTCTAACCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACTTCCTTAATCGCCGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCCTTCTAAAAGACGATGCTGTAACGCGAcc                                                                        >  1:1037009/1‑180 (MQ=255)
                                                            gTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTaga                                         <  1:1031298/210‑1 (MQ=31)
                                                            gTCTAACCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAGCGCACCCCTGATAAGGGTGAGGt                 <  2:714873/234‑1 (MQ=14)
                                                                aaCCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAg                                                       >  2:1080182/1‑192 (MQ=14)
                                                                  ccTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAgcgc                                     <  2:1252189/208‑1 (MQ=31)
                                                                  ccTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAGCGCAc                                   <  1:784363/210‑1 (MQ=14)
                                                                    ttCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGCAGTTCa       <  1:1092040/236‑1 (MQ=31)
                                                                    ttCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATAAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGCAGTTCa       <  1:584426/236‑1 (MQ=14)
                                                                    ttCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGCAGTTCa       >  1:278942/1‑236 (MQ=14)
                                                                     tCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCCGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAGCGCa                                    >  1:774108/1‑206 (MQ=14)
                                                                       ggggggACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGa                                                                          <  2:149204/166‑1 (MQ=31)
                                                                       gggaggACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGa                                                                          <  2:856121/166‑1 (MQ=14)
                                                                         ggggACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGGATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGCAGTTCAAATCt  <  1:343186/236‑1 (MQ=14)
                                                                         ggggACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGtt                                                                   >  2:1053345/1‑171 (MQ=31)
                                                                         ggggACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTa                                           >  2:87413/1‑195 (MQ=31)
                                                                         ggggACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGCAGTTCAAATCt  >  1:196521/1‑236 (MQ=31)
                                                                          |                                                                                                                                                                                                                                          
CCCGGGCCTTGTACACACCGCCCGTCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGCAGTTCAAATCT  >  AE015451/2550051‑2550359

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