Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 757,142 T→C 5.2% F115F (TTT→TTC PP_0647 → conserved membrane protein of unknown function

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451757,1420TC 5.2% 267.4 / 3.5 77F115F (TTT→TTCPP_0647conserved membrane protein of unknown function
Reads supporting (aligned to +/- strand):  ref base T (42/31);  new base C (2/2);  total (44/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.95e-01

AATTTCACTGTTGCAACCGCGGCTATCCAAAGTCATCCCGGCCTGGAGCAAATGAAACTGTTTTTGGGCACGCGGCGCTTTGGCGCTCGGTGGCTGCTGATAACGTTTCTTTGTGGAATACTCATTTGGCCTGGCGTACATATGCGCCGAGGGCATCTGGACATCGAAGAACTGGAAAACTTCCCTGCAAAACTTAAGTTGGCTCTACAGGTATCCGCGTGGCTGACGGCCTTTGGCTTCGCTTGGATGTGGGTGACCTATGTATTCATCAAGGCGTGAACGTTGATGCCACGGAGATCGTAGCTCTCACTGGATCCTTCGCAAGTCCCGATTCGATTCCACTGAAAGCCATAAGTCTGGTTCAGTTACACCACTTCCCGGGAGGCAATCATAGTGTCAGGCTAAAGTACATAGATTTGGGCCTTTGTCTCGTGGATATCCAGGTTGCCTTCTCAGGCCAG  >  AE015451/756909‑757369
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                    
aaTTTCACTGTTGCAACCGCGGCTATCCAAAGTCATCCCGGCCTGGAGCAAATGAAACTGTTTTTGGGCACGCGGCGCTTTGGCGCTCGGTGGCTGCTGATAACGTTTCTTTGTGGAATACTCATTTGGCCTGGCGTACATATGCGCCGAGGGCATCTGGACATCGAAGAACTGGAAAACTTCCCTGCAAAACTTAAGTTGGCTCTACAGGTATCCGCGTGGCTGACGGCCTTTgg                                                                                                                                                                                                                                    >  1:1061195/1‑236 (MQ=255)
          ttGCAACCGCGGCTATCCAAAGTCATCCCGGCCTGGAGCAAATGAAACTGTTTTTGGGCACGCGGCGCTTTGGCGCTCGGTGGCTGCTGATAACGTTTCTTTGTGGAATACTCATTTGGCCTGGCGTACATATGCGCCGAGGGCATCTGGACATCGAAGAACTGGAAAACTTCCCTGCAAAACTTAAGTTGGCTCTACAGGTATCCGCGTGGCTGACGGCCTTTGGCTTCGCt                                                                                                                                                                                                                             <  2:555075/233‑1 (MQ=255)
             cAACCGCGGCTATCCAAAGTCATCCCGGCCTGGAGCAAATGAAACTGTTTTTGGGCACGCGGCGCTTTGGCGCTCGGTGGCTGCTGATAACGTTTCTTTGTGGAATACTCATTTGGCCTGGCGTACATATGCGCCGAGGGCATCTGGACATCGAAGAACTGGAAAACTTCCCTGCAAAACTTAAGTTGGCTCTACAGGTATCCGCGTGGCTGACGGCCTTTGGCTTCGCTggg                                                                                                                                                                                                                         >  2:1133279/1‑233 (MQ=255)
                  gcgGCTATCCAAAGTCATCCCGGCCTGGAGCAAATGAAACTGTTTTTGGGCACGCGGCGCTTTGGCGCTCGGTGGCTGCTGATAACGTTTCTTTGTGGAATACTCATTTGGCCTGGCGTACATATGCGCCGAGGGCATCTGGACATCGAAGAACTGGAAAACTTCCCTGCAAAACTTAAGTTGGCTCTACAGGTATCCGCGTGGCTGACGGCCTTTGGCTTCGCTTGGATGTGGGt                                                                                                                                                                                                                 >  1:158081/1‑236 (MQ=255)
                       tATCCAAAGTCATCCCGGCCTGGAGCAAATGAAACTGTTTTTGGGCACGCGGCGCTTTGGCGCTCGGTGGCTGCTGATAACGTTTCTTTGTGGAATACTCATTTGGCCTGGCGTACATATGCGCCGAGGGCATCTGGACATCGAAGAACTGGAAAACTTCCCTGCAAAACTTAAGTTGGCTCTACAGGTATCCGCGTGGCTGACGGCCTTTGGCTTCGCTTGGATGTGGGTGa                                                                                                                                                                                                               <  1:299952/233‑1 (MQ=255)
                                   tCCCGGCCTGGAGCAAATGAAACTGTTTTTGGGCACGCGGCGCTTTGGCGCTCGGTGGCTGCTGATAACGTTTCTTTGTGGAATACTCATTTGGCCTGGCGTACATATGCGCCGAGGGCATCTGGACATCGAAGAACTGGAAAACTTCCCTGCAAAACTTAAGTTGGCTCTACAGGTATCCGCGTGGCTGACGGCCTTTGGCTTCGCTTGGAtgt                                                                                                                                                                                                                     <  2:916405/215‑1 (MQ=255)
                                     ccGGCCTGGAGCAAATGAAACTGTTTTTGGGCACGCGGCGCTTTGGCGCTCGGTGGCTGCTGATACCGTTTCTTTGTGGAATACTCATTTGGCCTGGCGTACATATGCGCCGAGGGCATCTGGACATCGAAGAACTGGAAAACTTCCCTGCAAAACTTAAGTTGGCTCTACAGGTATCCGCGTGGCTGACGGCCTTTGGCtt                                                                                                                                                                                                                                 >  2:596843/1‑202 (MQ=255)
                                      cGGCCTGGAGCAAATGAAACTGTTTTTGGGCACGCGGCGCTTTGGCGCTCGGTGGCTGCTGATAACGTTTCTTTGTGGAATACTCATTTGGCCTGGCGTACATATGCGCCGAGGGCATCTGGACATCGAAGAACTGGAAAACTTCCCTGCAAAACTTAAGTTGGCTCTACAGGTATCCGCGTGGCTGACGGCCTTTGGCTTCGCTTGGATGTggg                                                                                                                                                                                                                  >  1:130008/1‑215 (MQ=255)
                                      cGGCCTGGAGCAAATGAAACTGTTTTTGGGCACGCGGCGCTTTGGCGCTCGGTGGCTGCTGATAACGTTTCTTTGTGGAATACTCATTTGGCCTGGCGTACATATGCGCCGAGGGCATCTGGACATCGAAGAACTGGAAAACTTCCCTGCAAAACTTAAGTTGGCTCTACAGGTATCCGCGTGGCTGACGGCCTTTGGCTTCGCTTGGATGTGGGTGACCTAt                                                                                                                                                                                                          <  2:1058942/223‑1 (MQ=255)
                                      cGGCCTGGAGCAAATGAAACTGTTTTTGGGCACGCGGCGCTTTGGCGCTCGGTGGCTGCTGATAACGTTTCTTTGTGGAATACTCATTTGGCCTGGCGTACATATGCGCCGAGGGCATCTGGACATCGAAGAACTGGAAAACTTCCCTGCAAAACTTAAGTTGGCTCTACAGGTATCCGCGTGGCTGACGGCCTTTGGCTTCGCTTGGATGTGGGTGACCTAt                                                                                                                                                                                                          >  1:1058947/1‑223 (MQ=255)
                                      cGGCCTGGAGCAAATGAAACTGTTTTTGGGCACGCGGCGCTTTGGCGCTCGGTGGCTGCTGATAACGTTTCTTTGTGGAATACTCATTTGGCCTGGCGTACATATGCGCCGAGGGCATCTGGACATCGAAGAACTGGAAAACTTCCCTGCAAAACTTAAGTTGGCTCTACAGGTATCCGCGTGGCTGACGGCCTTTGGCTTCGCTTGGATGTGGGTGACCTATGTATtcatca                                                                                                                                                                                                >  2:1081975/1‑233 (MQ=255)
                                      cGGCCTGGAGCAAATGAAACTGTTTTTGGGCACGCGGCGCTTTGGCGCTCGGTGGCTGCTGATAACGTTTCTTTGTGGAATACTCATTTGGCCTGGCGTACATATGCGCCGAGGGCATCTGGACATCGAAGAACTGGAAAACTTCCCTGCAAAACTTAAGTTGGCTCTACAGGTATCCGCGTGGCTGACGGCCTTTGGCTTCGCTTGGATGTGGGTGACCTATGTATTCATCAAgg                                                                                                                                                                                             >  1:481913/1‑236 (MQ=255)
                                      cGGCCTGGAGCAAATGAAACTGTTTTTGGGCACGCGGCGCTTTGGCGCTCGGTGGCTGCTGATAACGTTTCTTTGTGGAATACTCATTTGGCCTGGCGTACATATGCGCCGAGGGCATCTGGACATCGAAGAACTGGAAAACTTCCCTGCAAAACTTAAGTTGGCTCTACAGGTATCCGCGTGGCTGACGGCCTTTGGCTTCGCTTGGATGTGGGTGACCTATGTATTCATCAAg                                                                                                                                                                                              >  1:119857/1‑235 (MQ=255)
                                             gggCAAAGGAAACGGTTTTTGGGCACGCGGCGCTTTGGCGCTCGGTGGCTGCTGATAACGTTTCTTTGTGGAATACTCATTTGGCCTGGCGTACATATGCGCCGAGGGCATCTGGACATCGAAGAACTGGAAAACTTCCCTGCAAAACTTAAGTTGGCTCTACAGGTATCCGCGTGGCTGACGGCCTTTGGCTTCGCTTGGATGTggg                                                                                                                                                                                                                  <  1:648761/206‑1 (MQ=255)
                                             gAGCAAATGAAACTGTTTTTGGGCACGCGGCGCTTTGGCGCTCGGTGGCTGCTGATAACGTTTCTTTGTGGAATACTCATTTGGCCTGGCGTACATATGCGCCGAGGGCATCTGGACATCGAAGAACTGGAAAACTTCCCTGCAAAACTTAAGTTGGCTCTACAGGTATCCGCGTGGCTGACGGCCTTTGGCTTCGCTTGGATGTGGGTGACCTATGTATTCATCAAGGCGTGaa                                                                                                                                                                                       <  2:958931/235‑1 (MQ=255)
                                                           gTTTTTGGGCACGCGGCGCTTTGGCGCTCGGTGGCTGCTGATAACGTTTCTTTGTGGAATACTCATTTGGCCTGGCGTACATATGCGCCGAGGGCATCTGGACATCGACGAACTGGAAAACTTCCCTGCAAAACTTAAGTTGGCTCTACAGGTATCCGCGTGGCTGACGGCCTTTGGCTTCGCTTGGATGTGGGTGACCTATGTATTCATCaa                                                                                                                                                                                               >  1:19741/1‑213 (MQ=255)
                                                               ttGGGCACGCGGCGCTTTGGCGCTCGGTGGCTGCTGATAACGTTTCTTTGTGGAATACTCATTTGGCCTGGCGTACATATGCGCCGAGGGCATCTGGACATCGAAGAACTGGAAAACTTCCCTGCAAAACTTAAGTTGGCTCTACAGGTATCCGCGTGGCTGACGGCCTTTGGCTTCGCTTGGATGTGGGTGACCTATGTATTCATCAAGGCGTGAACGTTGATGCCACGGAGATc                                                                                                                                                                    >  1:222290/1‑236 (MQ=255)
                                                                  ggCACGCGGCGCTTTGGCGGTCGGTGGCTGCTGATAACGTTTCTTTGTGGAATACTCATTTGGCCTGGCGTACATATGCGCCGAGGGCATCTGGACATCGAAGAACTGGAGAACTTCCCTGCAAAACGTAAGTTGGCTCTACAGGTATCCGCGTGGCTGACGGCCTTTGGCTTCGCTTGGATGTGGGTGTCCTATGTATTCATCAAGGCGTGAACGTTGATGCCACGGAGATCGTa                                                                                                                                                                 <  1:1111357/236‑1 (MQ=255)
                                                                         ggcgCTTTGGCGCTCGGTGGCTGCTGATAACGTTTCTTTGTGGAATACTCATTTGGCCTGGCGTACATATGCGCCGAGGGCATCTGGACATCGAAGAACTGGAAAACTTCCCTGCAAAACTTAAGTTGGCTCTACAGGTATCCGCGTGGCTGACGGCCTTTGGCTTCGCTTGGATGTGGGTGACCTATGTATTCATCAAGGCGTGAACGTTGAt                                                                                                                                                                                >  2:457739/1‑214 (MQ=255)
                                                                          gcgcTTTGGCGCTCGGTGGCTGCTGATAACGTTTCTTTGTGGAATACTCATTTGGCCTGGCGTACATATGCGCCGAGGGCATCTGGACATCGAAGAACTGGAAAACTTCCCTGCAAAACTTAAGTTGGCTCTACAGGTATCCGCGTGGCTGACGGCCTTTGGCTTCGCTTGGAtgtg                                                                                                                                                                                                                    <  1:344419/177‑1 (MQ=255)
                                                                                         gTGGCTGCTGATAACGTTTCTTTGTGGAATACTCATTTGGCCTGGCGTACATATGCGCCGAGGGCATCTGGACATCGAAGAACTGGAAAACTTCCCTGCAAAACTTAAGTTGGCTCTACAGGTATCCGCGTGGCTGACGGCCTTTGGCTTCGCTTGGATGTGGGTGACCTATGTATTCATCAAGGCGTGAACGTTGATGCCACGGAGATCGTa                                                                                                                                                                 <  1:16937/213‑1 (MQ=255)
                                                                                          tGGCTGCTGATAACGTTTCTTTGTGGAATACTCATTTGGCCTGGCGTACATATGCGCCGAGGGCATCTGGACATCGAAGAACTGGAAAACTTCCCTGCAAAACTTAAGTTGGCTCTACAGGTATCCGCGTGGCTGACGGCCTTTGGCTTCGCTTGGATGTGGGTGACCTATGTATTCATCAAGGCGTGAACGTTGATGCCACGGAGATCGTAGCTCTCACTGGATCCTTCGCAAGt                                                                                                                                         >  1:1032116/1‑236 (MQ=255)
                                                                                                             tttGTGGAATACTCATTTGGCCTGGCGTACATATGCGCCGAGGGCATCTGGACATCGAAGAACTGGAAAACTTCCCTGCAAAACTTAAGTTGGCTCTACAGGTATCCGCGTGGCTGACGGCCTTTGGCTTCGCTTGGATGTGGGTGACCTATGTATTCATCAAGGCGTGAACGTTGATGcc                                                                                                                                                                             >  2:514622/1‑181 (MQ=255)
                                                                                                             tttGTGGAATACTCATTTGGCCTGGCGTACATATGCGCCGAGGGCATCTGGACATCGAAGAACTGGAAAACTTCCCTGCAAAACTTAAGTTGGCTCTACAGGTATCCGCGTGGCTGACGGCCTTTGGCTTCGCTTGGATGTGGGTGACCTATGTATTCATCAAGGCGTGAACGTTGATGCCACGGAGATCGTAGCTCTCACTGGATCCTTCGCAAGTCCCGATTCGATTCCACTGa                                                                                                                      >  1:635516/1‑236 (MQ=255)
                                                                                                              ttGTGGAATACTCATTTGGCCTGGCGTACATATGCGCCGAGGGCATCTGGACATCGAAGAACTGGAAAACTTCCCTGCAAAACTTAAGTTGGCTCTACAGGTATCCGCGTGGCTGACGGCCTTTGGCTTCGCTTGGATGTGGGt                                                                                                                                                                                                                 >  2:1007210/1‑144 (MQ=255)
                                                                                                              ttGTGGAATACTCATTTGGCCTGGCGTACATATGCGCCGAGGGCATCTGGACATCGAAGAACTGGAAAACTTCCCTGCAAAACTTAAGTTGGCTCTACAGGTATCCGCGTGGCTGACGGCCTTTGGCTTCGCTTGGATGTGGGt                                                                                                                                                                                                                 <  1:1007215/144‑1 (MQ=255)
                                                                                                              ttGTGGAATACTCATTTGGCCTGGCGTACATATGCGCCGAGGGCATCTGGACATCGAAGAACTGGAAAACTTCCCTGCAAAACTTAAGTTGGCTCTACAGGTATCCGCGTGGCTGACGGCCTTTGGCTTCGCTTGGATGTGGGTGACCTATGTATTCATCAAGGCGTGAACGTTGATGCCCCg                                                                                                                                                                          >  1:96791/1‑183 (MQ=255)
                                                                                                                  ggAATACTCATTTGGCCTGGCGTACATATGCGCCGAGGGCATCTGGACATCGAAGAACTGGAAAACTTCCCTGCAAAACTTAAGTTGGCTCTACAGGTATCCGCGTGGCTGACGGTCTTTGGCTTCGCTTGGATGTGGGTGACCTATGTATTCATCATGGCGTGAACGtt                                                                                                                                                                                   >  2:948440/1‑170 (MQ=255)
                                                                                                                  ggAATACTCATTTGGCCTGGCGTACATATGCGCCGAGGGCATCTGGACATCGAAGAACTGGAAAACTTCCCTGCAAAACTTAAGTTGGCTCTACAGGTATCCGCGTGGCTGACGGCCTTTGGCTTCGCTTGGATGTGGGTGACCTATGTATTCATCaa                                                                                                                                                                                               <  2:850314/158‑1 (MQ=255)
                                                                                                                  ggAATACTCATTTGGCCTGGCGTACATATGCGCCGAGGGCATCTGGACATCGAAGAACTGGAAAACTTCCCTGCAAAACTTAAGTTGGCTCTACAGGTATCCGCGTGGCTGACGGCCTTTGGCTTCGCTTGGATGTGGGTGACCTATGTATTCATCaa                                                                                                                                                                                               >  1:850319/1‑158 (MQ=255)
                                                                                                                   gAATACTCATTTGGCCTGGCGTACATATGCGCCGAGGGCATCTGGACATCGAAGAACTGGAAAACTTCCCTGCAAAACTTAAGTTGGCTCTACAGGTATCCGCGTGGCTGACGGCCTTTGGCTTCGCTTGGATGTGGGTGACCTATGTATTCATCAAGGCGTGAACGTTGATGCCACGGAGATCGTAGCTCTCACTGGATCCTTCGCAAGTCCcg                                                                                                                                     <  2:622004/215‑1 (MQ=255)
                                                                                                                       aCTCATTTGGCCTGGCGTACATATGCGCCGAGGGCATCTGGACATCGAAGAACTGGAAAACTTCCCTGCAAAACTTAAGTTGGCTCTACAGGTATCCGCGTGGCTGACGGCCTTTGGCTTCGCTTGGATGTGGGTGACCTATGTATTCATCAAGGCGTGAACGTTGATGCCACGGAGTTCGTAGCTCTCACTGGATCCTTCGCAAGTCCCGCTTCGATTCCACt                                                                                                                        >  2:837376/1‑224 (MQ=255)
                                                                                                                        cTCATTTGGCCTGGCGTACATATGCGCCGAGGGCATCTGGACATCGAAGAACTGGAAAACTTCCCTGCAAAACTTAAGTTGGCTCTACAGGTATCCGCGTGGCTGACGGCCTTTGGCTTCGCTTGGATGTGGGTGACCTATGTATTCATCAAGGCGTGAACGTTGATGCCACGGAGATCGTAGCTCTCACTGGATCCTTCGCAAGTCCcg                                                                                                                                     <  2:1061190/210‑1 (MQ=255)
                                                                                                                         tCATTTGGCCTGGCGTACATATGCGCCGAGGGCATCTGGACATCGAAGAACTGGAAAACTTCCCTGCAAAACTTAAGTTGGCTCTACAGGTATCCGCGTGGCTGACGGCCTTTGGCTTCGCTTGGATGTGGGTGACCTATGTATTCATCAAGGCGTGAACGTTGATGCCACGGAGATCGTAGCTCTCACTGGATCCTTCGCAAGTCCcgat                                                                                                                                   >  1:520450/1‑211 (MQ=255)
                                                                                                                         tCATTTGGCCTGGCGTACATATGCGCCGAGGGCATCTGGACATCGAAGAACTGGAAAACTTCCCTGCAAAACTTAAGTTGGCTCTACAGGTATCCGCGTGGCTGACGGCCTTTGGCTTCGCTTGGATGTGGGTGACCTATGTATTCATCAAGGCGTGAACGTTGATGCCACGGAGATCGTAGCTCTCACTGGATCCTTCGCAAGTCCCGATTCGATTCCACTGAAAGCCATAAGt                                                                                                           >  1:701359/1‑235 (MQ=255)
                                                                                                                         tCATTTGGCCTGGCGTACATATGCGCCGAGGGCATCTGGACATCGAAGAACTGGAAAACTTCCCTGCAAAACTTAAGTTGGCTCTACAGGTATCCGCGTGGCTGACGGCCTTTGGCTTCGCTTGGATGTGGGTGACCTATGTATTCATCAAGGCGTGAACGTTGATGCCACGGAGATCGTAGCTCTCACTGGATCCTTCGCAAGTCCCGATTCGATTCCACTGAAAGCCATAAGTc                                                                                                          >  1:853014/1‑236 (MQ=255)
                                                                                                                                   tGGCGTACATATGCGCCGAGGGCATCTGGACATCGAAGAACTGGAAAACTTCCCTGCAAAACTTAAGTTGGCTCTACAGGTATCCGCGTGGCTGACGGCCTTTGGCTTCGCTTGGATGTGGGTGACCTATGTATTCATCAAGGCGTGAACGTTGATGCCACGGAGATCGTAGCTCTCACTGGATCCTTCGCAAGTCCCGATTCGATTCCACTGAAAGCCATa                                                                                                              <  2:173916/222‑1 (MQ=255)
                                                                                                                                      cGTACATATGCGCCGAGGGCATCTGGACATCGAAGAACTGGAAAACTTCCCTGCAAAACTTAAGTTGGCTCTACAGGTATCCGCGTGGCTGACGGCCTTTGGCTTCGCTTGGATGTGGGTGACCTATGTATTCATCAAGGCGTGAACGTTGATGCCACGGAGATCGTAGCTCTCACTGGATCCTTCGCAAGTCCCGATTCGATTCCACTGAAAGCCATAAGTCTGGTTCAGTTaca                                                                                             <  2:222287/236‑1 (MQ=255)
                                                                                                                                       gTACATATGCGCCGAGGGCATCTGGACATCGAAGAACTGGAAAACTTCCCTGCAAAACTTAAGTTGGCTCTACAGGTATCCGCGTGGCTGACGGCCTTTGGCTTCGCTTGGATGTGGGTGACCTATGTATTCATCAAGGCGTGAACGTTGATGCCACGGAGATCGTAGCTCTCACTGGATCCTTCGCAAGTCCCGATTCGATTCCACTGAAAGCCATAAGTCTGGTTCAGTTacac                                                                                            >  1:1236316/1‑236 (MQ=255)
                                                                                                                                        tACATATGCGCCGAGGGCATCTGGACATCGAAGAACTGGAAAACTTCCCTGCAAAACTTAAGTTGGCTCTACAGGTATCCGCGTGGCTGACGGCCTTTGGCTTCGCTTGGGATGTGGGTGACCTATGTATTCATCAAGGCGTGAACGTTGATGCCACGGAGATCGTAGCTCTCACTGGATCCTTCGCAAGTCCCGATTCGATTCCACTGAAAGCCATAAGTCTGGTTCAGTTacac                                                                                            >  1:412187/1‑236 (MQ=255)
                                                                                                                                        tACATATGCGCCGAGGGCATCTGGACATCGAAGAACTGGAAAACTTCCCTGCAAAACTTAAGTTGGCTCTACAGGTATCCGCGTGGCTGACGGCCTTTGGCTTCGCTTGGATGTGGGTGACCTATGTATTCATCAAGGCGTGAACGTTGATGCCACGGAGATCGTAGCTCTCACTGGATCCTTCGCAAGTCCCGATTCGATTCCACTGAAAGCCa                                                                                                                >  2:291639/1‑215 (MQ=255)
                                                                                                                                                  ccGAGGGCATCTGGACATCGAAGAACTGGAAAACTTCCCTGCAAAACTTAAGTTGGCTCTACAGGTATCCGCGTGGCTGACGGCCTTTGGCTTCGCTTGGATGTGGGTGACCTATGTATTCATCAAGGCGTGAACGTTGATGCCACGGAGATCGTAGCTCTCACTGGATCCTTCGCAAGTCCCGATTCGATTCCACTGAAAGCCATAAGTCTGGTTCAGTTACACCACTTCCCggg                                                                                 <  1:1205086/236‑1 (MQ=255)
                                                                                                                                                  ccGAGGGCATCTGGACATCGAAGAACTGGAAAACTTCCCTGCAAAACTTAAGTTGGCTCTACAGGTATCCGCGTGGCTGACGGCCTTTGGCTTCGCTTGGATGTGGGTGACCTATGTATTCATCAAGGCGTGAACGTTGATGCCACGGAGATCGTAGCTCTCACTGGATCCTTCGCAAGTCCCGATTCGATTCCACTGAAAGCCATAAGTCTGGTTCAGTTACACCACTTCCCgg                                                                                  >  1:1160597/1‑235 (MQ=255)
                                                                                                                                                  ccGAGGGCATCTGGACATCGAAGAACTGGAAAACTTCCCTGCAAAACTTAAGTTGGCTCTACAGGTATCCGCGTGGCTGACGGCCTTTGGCTTCGCTTGGATGTGGGTGACCTATGTATTCATCAAGGCGTGAACGTTGATGCCACGGAGATCGTAGCTCTCACTGGATCCTTCGCAAGTCCCAATTCGATTCCACTAAAAGCCATaa                                                                                                             >  1:879379/1‑208 (MQ=255)
                                                                                                                                                  ccGAGGGCATCTGGACATCGAAGAACTGGAAAACTTCCCTGCAAAACTTAAGTTGGCTCTACAGGTATCCGCGTGGCTGACGGCCTTCGGCTTCGCTTGGATGTGGGTGACCTATGTATTCATCAAGGCGTGAACGTTGATGCCACGGAGATCGTAGCTCTCACTGGATCCTTCGCAAGTCCCGATTCGATTCCACTGAAAGCCATAAGTCTGGTTCAGTTACACCACTTCCCg                                                                                   >  2:516017/1‑234 (MQ=255)
                                                                                                                                                   cGAGGGCATCTGGACATCGAAGAACTGGAAAACTTCCCTGCAAAACTTAAGTTGGCTCTACAGGTATCCGCGTGGCTGACGGCCTTTGGCTTCGCTTGGATGTGGGTGACCTATGTATTCATCAAGGCGTGAACGTTGATGCCACGGAGATCGTAGCTCTCACTGGATCCTTCGCAAGTCCCGATTCGATTCCACTGAAAGCCATAAGTCTGGTTCAGTTACACCACTTCCCGGGa                                                                                <  1:114468/236‑1 (MQ=255)
                                                                                                                                                     aGGGCATCTGGACATCGAAGAACTGGAAAACTTCCCTGCAAAACTTAAGTTGGCTCTACAGGTATCCGCGTGGCTGACGGCCTTTGGCTTCGCTTGGATGTGGGTGACCTATGTATTCATCAAGGCGTGAACGTTGATGCCACGGAGATCGTAGCTCTCACTGGATCCTTCGCAAGTCCCGATTCGATTCCACTGAAAGCCATAAGTCTGGTTCAGTTACACCACTTCCCGGGAg                                                                               >  1:1086318/1‑235 (MQ=255)
                                                                                                                                                        gCATCTGGACATCGAAGAACTGGAAAACTTCCCTGCAAAACTTAAGTTGGCTCTACAGGTATCCGCGTGGCTGACGGCCTTTGGCTTCGCTTGGATGTGGGTGACCTATGTATTCATCAAGGCGTGAACGTTGATGCCACGGAGATCGTAGCTCTCACTGGATCCTTCGCAAGTCCCGATTCGATTCCACTGAAAGCCATACGTCTGGTTCAGTTACACCACTTcc                                                                                     >  2:9053/1‑226 (MQ=255)
                                                                                                                                                        gCATCTGGACATCGAAGAACTGGAAAACTTCCCTGCAAAACTTAAGTTGGCTCTACAGGTATCCGCGTGGCTGACGGCCTTTGGCTTCGCTTGGATGTGGGTGACCTATGTATTCATCAAGGCGTGAACGTTGATGCCACGGAGATCGTAGCTCTCACTGGATCCTTCGCAAGTCCCGATTCGATTCCACTGAAAGCCATAAGTCTGGTTCAGTTACACCACTTCCCGGGAGGCaa                                                                           >  1:30513/1‑236 (MQ=255)
                                                                                                                                                        gCATCTGGACATCGAAGAACTGGAAAACTTCCCTGCAAAACTTAAGTTGGCTCTACAGGTATCCGCGTGGCTGACGGCCTTCGGCTTc                                                                                                                                                                                                                                >  1:390240/1‑88 (MQ=255)
                                                                                                                                                        gCATCTGGACATCGAAGAACTGGAAAACTTCCCTGCAAAACTTAAGTTGGCTCTACAGGTATCCGCGTGGCTGACGGCCTTCGGCTTc                                                                                                                                                                                                                                <  2:390235/88‑1 (MQ=255)
                                                                                                                                                              ggACATCGAAGAACTGGAAAACTTCCCTGCAAAACTTAAGTTGGCTCTACAGGTATCCGCGTGGCTGACGGCCTTCGGCTTCGCTTGGATGTGGGTGACCTATGTATTCATCAAGGCGTGAACGTTGATGCCACGGAGATCGTAGCTCTCACTGGATCCTTCGCAAGTCCCGATTCGATTCCACTGAAAGCCATAAGTCTGGTTCAGTTACACCACTTCCCGGGAGGCAATCATAg                                                                     <  1:516022/236‑1 (MQ=255)
                                                                                                                                                                  aTCGAAGAACTGGAAAACTTCCCTGCAAAACTTAAGTTGGCTCTACAGGTATCCGCGTGGCTGACGGCCTTTGGCTTCGCTTGGATGTGGGTGACCTATGTATTCATCAAGGCGTGAACGTTGATGCCACGGAGATCGTAGCTCTCACTGGATCCTTCGCAAGTCCCGATTCGATTCCACTGAAAGCCATAAGTCTGGTTCAGTTACACCCCTTCCCgg                                                                                  >  2:1112381/1‑219 (MQ=255)
                                                                                                                                                                  aTCGAAGAACTGGAAAACTTCCCTGCAAAACTTAAGTTGGCTCTACAGGTATCCGCGTGGCTGACGGCCTTTGGCTTCGCTTGGATGTGGGTGACCTATGTATTCATCAAGGCGTGAACGTTGATGCCACGGAGATCGTAGCTCTCACTGGATCCTTCGCAAGTCCCGATTCGATTCCACTGAAAGCCATAAGTCTGGTTCAGTTACACCACTTc                                                                                      >  2:664818/1‑215 (MQ=255)
                                                                                                                                                                     gaagaaCTGGAAAACTTCCCTGCAAAACTTAAGTTGGCTCTACAGGTATCCGCGTGGCTGACGGCCTTTGGCTTCGCTTGGATGTGGGTGACCTATGTATTCATCAAGGCGTGAACGTTGATGCCACGGAGATCGTAGCTCTCACTGGATCCTTCGCAAGTCCCGATTCGATTCCACTGAAAGCCATAAGTCTGGTTCAGTTACACCACTTCCCGGGAGGCAATCATAGTGTCAgg                                                              <  1:664823/236‑1 (MQ=255)
                                                                                                                                                                       agaaCTGGAAAACTTCCCTGCAAAACTTAAGTTGGCTCTACAGGTATCCGCGTGGCTGACGGCCTTTGGCTTCGCTTGGATGTGGGTGACCTATGTATTCATCAAGGCGTGAACGTTGATGCCACGGAGATCGTAGCTCTCACTGGATCCTTCGCAAGTCCCGATTCGATTCCACTGAAAGCCATAAGTCTGGTTCAGTTACACCACTTCCCGGGAGGCAATCATAGTGTCAgg                                                              <  2:853009/234‑1 (MQ=255)
                                                                                                                                                                            tGGAAAACTTCCCTGCAAAACTTAAGTTGGCTCTACAGGTATCCGCGTGGCTGACGGCCTTTGGCTTCGCTTGGATGTGGGTGACCTATGTATTCATCAAGGCGTGAACGTTGATGCCACGGAGATCGTAGCTCTCACTGGATCCTTCGCACGTCCCGATTCGATTCCACTGAAAGCCATAAGTCTgg                                                                                                       >  2:42054/1‑188 (MQ=255)
                                                                                                                                                                                 aaCTTCCCTGCAAAACTTGAGTTGGCTCTACAGGTATCCGCGTGGCTGACGGCCTTTGGCTTCGCTTGGATGTggg                                                                                                                                                                                                                  >  2:1225095/1‑76 (MQ=255)
                                                                                                                                                                                 aaCTTCCCTGCAAAACTTAAGTTGGCTCTACAGGTATCCGCGTGGCTGACGGCCTTTGGCTTCGCTTGGATGTggg                                                                                                                                                                                                                  <  1:1225100/76‑1 (MQ=255)
                                                                                                                                                                                    ttCCCTGCAAAACTTAAGTTGGCTCTACAGGTATCCGCGTGGCTGACGGCCTTTGGCTTCGCTTGGATGTGGGTGACCTATGTATTCATCAAGGCGTGAACGTTGATGCCACGGAGATCGTAGCTCTCACTGGATCCTTCGCAAGTCCCGATTCGATTCCACTGAAAGCCATAAGTCTGGTTCAGTTACACCACTTCCCGGGAGGCAATCATAGTGTCAgg                                                              <  2:296578/221‑1 (MQ=255)
                                                                                                                                                                                      cccTGCAAAACTTAAGTTGGCTCTACAGGTATCCGCGTGGCTGACGGCCTTTGGCTTCGCTTGGATGTGGGTGACCTATGTATTCATCAAGGCGTGAACGTTGATGCCACGGAGATCGTAGCTCTCACTGGATCCTTCGCAAGTCCCGATTCGATTCCACTGAAAGCCATAAGTCTGGTTCAGTTACACCACTTCCCGGGAGGCAATCATAGTGTCAGGCTAAAGt                                                       <  2:496323/226‑1 (MQ=255)
                                                                                                                                                                                               aCTTAAGTTGGCTCTACAGGTATCCGCGTGGCTGACGGCCTTTGGCTTCGCTTGGATGTGGGTGACCTATGTATTCATCAAGGCGTGAACGTTGATGCCACGGAGATCGTAGCTCTCACTGGATCCTTCGCAAGTCCCGATTCGATTCCACTGAAAGCCATAAGTCTGGTTCAGTTACACCACTTCCCGGGAGGCAATCATAGTGTc                                                                 <  2:984070/207‑1 (MQ=255)
                                                                                                                                                                                                 ttAAGTGGGCTCTACAGGTATCCGCGTGGCTGACGGCCTTTGGCTTCGCTTGGATGTGGGTGACCTATGTATTCATCAAGGCGTGAACGTTGATGCCACGGAGATCGTAGCTCTCACTGGATCCTTCGCAAGTCCCGATTCGATTCCACt                                                                                                                        <  2:7455/150‑1 (MQ=255)
                                                                                                                                                                                                      ttGGCTCTTCAGGTATCCGCGTGGCTGACGGCCTTTGGCTTCGCTTGGATGTGGGTGACCTATGTATTCATCAAGGCGTGAACGTTGATGCCACGGAGATCGTAGCTCTCACTGGATCCTTCGCAAGTCCCGATTCGATTCCACTGAAAGCCATAAGTCTGGTTCAGTTACACCACTTCCCGGGAGGCAATCATAGTGTCAGGCTAAAGTACATAg                                                 <  2:1034906/216‑1 (MQ=255)
                                                                                                                                                                                                        ggCTCTACAGGTATCCGCGTGGCTGACGGCCTTTGGCTTCGCTTGGATGTGGGTGACCTATGTATTCATCAAGGCGTGAACGTTGATGCCACGGAGATCGTAGCTCTCACTGGATCCTTCGCAAGTCCCGATTCGATTCCACTGAAAGCCATAAGTCTGGTTCAGTTACACCACTTCCCGGGAGGCAATCATAGTGTCAGGCTAAAGt                                                       >  2:422514/1‑208 (MQ=255)
                                                                                                                                                                                                            ctACAGGTATCCGCGTGGCTGACGGCCTTTGGCTTCGCTTGGATGTGGGTGACCTATGTATTCATCAAGGCGTGAACGTTGATGCCACGGAGATCGTAGCTCTCACTGGATCCTTCGCAAGTCCCGATTCGATTCCACTGAAAGCCATAAGTCTGGTTCAGTTACACCACTTCCCGGGAGGCAATCATAGTGTCAGGCTAAAGTACATAGATTTGGGCCTTTGTCTCGTGGATATc                       <  1:81112/236‑1 (MQ=255)
                                                                                                                                                                                                             tACAGGTATCCGCGTGGCTGACGGCCTTTGGCTTCGCTTGGATGTGGGTGACCTTTGTATTCATCAAGGCGTGAACGTTGATGCCACGGAGATCGTAGCTCTCACTGGATCCTTCGCAAGTCCCGATTCGATTCCACTGAAAGCCATAAGTCTGGTTCAGTTACACCACTTcc                                                                                     >  2:768772/1‑173 (MQ=255)
                                                                                                                                                                                                             tACAGGTATCCGCGTGGCTGACGGCCTTTGGCTTCGCTTGGATGTGGGTGACCTATGTATTCATCAAGGCGTGAACGTTGATGCCACGGAGATCGTAGCTCTCACTGGATCCTTCGCAAGTCCCGATTCGATTCCACTGAAAGCCATAAGTCTGGTTCAGTTACACCACTTCCCGGGAGGCAATCATAGTGTCAGGCTAAAGTACATAGAt                                               <  2:956615/211‑1 (MQ=255)
                                                                                                                                                                                                                aGGTATCCGCGTGGCTGACGGCCTTTGGCTTCGCTTGGATGTGGGTGACCTAt                                                                                                                                                                                                          >  2:183578/1‑53 (MQ=255)
                                                                                                                                                                                                                aGGTATCCGCGTGGCTGACGGCCTTTGGCTTCGCTTGGATGTGGGTGACCTAt                                                                                                                                                                                                          <  1:183581/53‑1 (MQ=255)
                                                                                                                                                                                                                        gcgTGGCTGACGGCCTTTGGCTTCGCTTGGATGTGGGTGACCTATGTATTCATCAAGGCGTGAACGTTGATGCCACGGAGATCGTAGCTCTCACTGGATCCTTCGCAAGTCCCGATTCGATTCCACTGAAAGCCATAAGTCTGGTTCAGTTACACCACTTCCCGGGAGGCAATCATAGTGTCAGGCTAAAGTACATAGATTTGGGCCTTTGTCTCGTGGATATCCAGGTTGCCt             <  2:481908/234‑1 (MQ=255)
                                                                                                                                                                                                                          gTGGCTGACGGCCTTTGGCTTCGCTTGGATGTGGGTGACCTATGTATTCATCAAGGCGTGAACGTTGATGCCACGGAGATCGTAGCTCTCACTGGATCCTTCGCAAGTCCCGATTCGATTCCACTGAAAGCCATAAGTCTGGTTCAGTTACACCACTTCCCGGGAGGCAATCATAGTGTCAGGCTAAAGTACATAGATTTg                                            >  2:1245952/1‑201 (MQ=255)
                                                                                                                                                                                                                          gTGGCTGACGGCCTTTGGCTTCGCTTGGATGTGGGTGACCTATGTATTCATCAAGGCGTGAACGTTGATGCCACGGAGATCGTAGCTCTCACTGGATCCTTCGCAAGTCCCGATTCGATTCCACTGAAAGCCATAAGTCTGGTTCAGTTACACCACTTCCCGGGAGGCAATCATAGTGTCAGGCTAAAGTACATAGATTTGGGCCTTTGTCTCGTGGATATCCAGGTTGCCTtct          >  2:613534/1‑235 (MQ=255)
                                                                                                                                                                                                                            ggCTGACGGCCTTTGGCTTCGCTTGGATGTGGGTGACCTATGTATTCATCAAGGCGTGAACGTTGATGCCACGGAGATCGTAGCTCTCACTGGATCCTTCGCAAGTCCCGATTCGATTCCACTGAAAGCCATAAGTCTGGTTCAGTTACACCACTTCCCGGGAGGCAATCATAGTGTCAGGCTAAAGTACATAGATTTGGGCCTTTGTCTCGTGGATATCCAGGTTGCCTTCTCac       >  2:903587/1‑235 (MQ=255)
                                                                                                                                                                                                                            ggCTGACGGCCTTTGGCTTCGCTTGGATGTGGGTGACCTATGTATTCATCAAGGCGTGAACGTTGATGCCACGGAGATCGTAGCTCTCACTGGATCCTTCGCAAGTCCCGATTCGATTCCACTGAAAGCCATAAGTCTGGTTCAGTTACACCACTTCCCGGGAGGCAATCATAGTGTCAGGCTAAAGTACATAGATTTGGGCCTTTGTCTCGTGGATATCCAGGTTGCCTTCTCAg       >  1:341595/1‑236 (MQ=255)
                                                                                                                                                                                                                            ggCTGACGGCCGTTGGCTTCGCTTGGATGTGGGTGACCTATGTATTCATCAAGGCGTGAACGTTGATGCCACGGAGATCGTAGCTCTCACTGGATCCTTCGCAAGTCCCGATTCGATTCCACTGAAAGCCATAAGTCTGGTTCAGTTACACCACTTCCCGGGAGGCAATCATAGTGTCAGGCTAAAGTACATAGATTTGGGCCTTTGTCTCGTGGata                         <  2:320755/218‑1 (MQ=255)
                                                                                                                                                                                                                                   ggCCTTTGGCTTCGCTTGGATGTGGGTGCCCTATGTATTCATCAAGGCGTGAACGTTGATGCCACGGAGATCGTAGCTCTCACTGGATCCTTCGCAAGTCCCGATTCGATTCCACTGAAAGCCATAAGTCTGGTTCAGTTACACCACTTCCCGGGAGGCAATCATAGTGTCAGGCTAAAGTACATAGATTTGGGCCTTTGTCTCGTGGATATCCAGGTTGCCTTCTCAGGCCAg  <  1:510967/234‑1 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                    
AATTTCACTGTTGCAACCGCGGCTATCCAAAGTCATCCCGGCCTGGAGCAAATGAAACTGTTTTTGGGCACGCGGCGCTTTGGCGCTCGGTGGCTGCTGATAACGTTTCTTTGTGGAATACTCATTTGGCCTGGCGTACATATGCGCCGAGGGCATCTGGACATCGAAGAACTGGAAAACTTCCCTGCAAAACTTAAGTTGGCTCTACAGGTATCCGCGTGGCTGACGGCCTTTGGCTTCGCTTGGATGTGGGTGACCTATGTATTCATCAAGGCGTGAACGTTGATGCCACGGAGATCGTAGCTCTCACTGGATCCTTCGCAAGTCCCGATTCGATTCCACTGAAAGCCATAAGTCTGGTTCAGTTACACCACTTCCCGGGAGGCAATCATAGTGTCAGGCTAAAGTACATAGATTTGGGCCTTTGTCTCGTGGATATCCAGGTTGCCTTCTCAGGCCAG  >  AE015451/756909‑757369

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