Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 4,022,394 C→T 5.9% intergenic (‑89/+579) PP_3547 ← / ← emrB putative NAD(P)‑binding oxidoreductase with NAD(P)‑binding Rossmann‑fold domain/multidrug efflux transport system ‑ permease subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154514,022,3940CT 5.9% 188.3 / 5.0 68intergenic (‑89/+579)PP_3547/emrBputative NAD(P)‑binding oxidoreductase with NAD(P)‑binding Rossmann‑fold domain/multidrug efflux transport system ‑ permease subunit
Reads supporting (aligned to +/- strand):  ref base C (27/37);  new base T (2/2);  total (29/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.35e-01

AACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCGAAAAGGTGAGTCAGCGTCGCCATCGCAAACGGACGGTTCGCTGAGTCAAAACAGACGAGTTCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCACGCAGGTGGCGCATA  >  AE015451/4022178‑4022626
                                                                                                                                                                                                                        |                                                                                                                                                                                                                                        
aaCATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGcc                                                                                                                                                                                                                              >  2:1153468/1‑229 (MQ=255)
aaCATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGgatgat                                                                                                                                                                                                                       >  1:1306101/1‑236 (MQ=255)
aaCATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGgat                                                                                                                                                                                                                          >  2:1033599/1‑233 (MQ=255)
aaCATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCAGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGgatgat                                                                                                                                                                                                                       >  1:1336338/1‑236 (MQ=255)
aaCATCAGCACAGCCAACCCGCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCACGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGgat                                                                                                                                                                                                                          >  2:918720/1‑233 (MQ=255)
 aCATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCAGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGgatga                                                                                                                                                                                                                        <  2:1336335/234‑1 (MQ=255)
     cAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAgg                                                                                                                                                                                                                                  <  2:1058499/220‑1 (MQ=255)
                  ccTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGtttttt                                                                                                                                                                                                     <  1:560124/236‑1 (MQ=255)
                          tttgatctggcttttgatctggctTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACtt                                                                                                                                                                                             >  1:1410651/1‑236 (MQ=255)
                            tgatctggcttttgatctggctttTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATCCATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGcc                                                                                                                                                                                                                              <  1:1420497/200‑1 (MQ=14)
                                    cttttgatctggcttttgatctTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGt                                                                                                                                                                                                          <  1:922593/213‑1 (MQ=38)
                                      tttgatctggcttttgatctTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGtt                                                                                                                                                                                  <  1:918722/235‑1 (MQ=37)
                                                                tgatcttgatcttCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAggg                                                                                                                                                        <  1:1148048/235‑1 (MQ=31)
                                                                    cttgatcttCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGtgt                                                                                                                                                                                                        >  1:90214/1‑181 (MQ=14)
                                                                           ttCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACaaaa                                                                                                                                                                          <  2:844089/206‑1 (MQ=31)
                                                                               cgACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCGAAAAGGTGAGTCa                                                                                                                                         >  1:498867/1‑235 (MQ=33)
                                                                                gACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCGAAAAGGTGAGTCAGc                                                                                                                                       >  1:778832/1‑236 (MQ=33)
                                                                                gACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCGAAAAGGTGAGTCAGc                                                                                                                                       >  1:572882/1‑236 (MQ=33)
                                                                                                  gCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATCCTGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTcg                                                                                                                                                                                      <  1:1068626/171‑1 (MQ=9)
                                                                                                     ggCCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACaaaa                                                                                                                                                                          <  1:1043815/180‑1 (MQ=17)
                                                                                                                       gACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCGAAAAGGTGAGTCAGCGTCGCCATCGCAAACGGACGGTTCGCTGAGTCaaaa                                                                                                   >  2:772697/1‑233 (MQ=38)
                                                                                                                           ggCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCGAAAAGGTGAGTCAGCGTCGCCATCGCAAACGGACGGTTCg                                                                                                              >  2:1008121/1‑218 (MQ=38)
                                                                                                                            gCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGcc                                                                                                                                                                                                                              <  2:1510638/105‑1 (MQ=11)
                                                                                                                            gCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGcc                                                                                                                                                                                                                              >  1:1510641/1‑105 (MQ=11)
                                                                                                                                  aTGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCGAAAAGGTGAGTCAGCGTCGCCATCGCAAACGGACGGTTCGCTGAGTCAAAACAGACGAGTTCg                                                                                       >  2:1104194/1‑234 (MQ=40)
                                                                                                                                   tGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTCGCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACa                                                                                                                                                                             <  2:350145/146‑1 (MQ=12)
                                                                                                                                          cAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCAGTAAGGGCGAAAAGGTg                                                                                                                                              >  2:1395779/1‑171 (MQ=17)
                                                                                                                                             gcgcTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACttttt                                                                                                                                                                                          >  2:1140908/1‑124 (MQ=32)
                                                                                                                                                      ccccAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCGAAAAGGTGAGTCAGCGTCGCCATCGCAAACGGACGGTTCGCTGAGTCAAAACAGACGAGTTCGCAtcttgcttcttgcttcttg                                                                  >  2:221497/1‑235 (MQ=255)
                                                                                                                                                         cAGGATGGGGCCGTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCGAAAAGGTGAGTCAGCGTCGCCATCGCAAACGGACGGTTCGCTGCGTCAAAACAGACGAGTTCGc                                                                                      <  2:1512151/212‑1 (MQ=255)
                                                                                                                                                           ggATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCGAAAAGGTGAGTCAGCGTTGCCATCGCAAACGGACGGTTCGCTGAGTCAAAACAGACGAGTTCGCAtcttgcttcttgcttcttgcttc                                                              <  2:778830/234‑1 (MQ=255)
                                                                                                                                                            gATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCGAAAAGGTGAGTCAGCGTCGCCATCGCAAACGGACGGTTCGCTGAGTCAAAACAGACGAGTTCGCAtcttgcttcttgcttcttgcttcttg                                                           <  1:385255/236‑1 (MQ=255)
                                                                                                                                                                   ccTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGGTACTTTTTGTCGCGTGTGACAAACAgg                                                                                                                                                                       >  2:557418/1‑120 (MQ=255)
                                                                                                                                                                     ttCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCGGAAGGGTGactatgc                                                                                                                                       <  2:201539/151‑7 (MQ=9)
                                                                                                                                                                       cggcggGGTCCTGCCGGGAGCAAGGGCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGACGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCGAAAAGGTGAGTCAGCGTCGCCATCGCAAACGGACGGTTCGCTGAGTCAAAACAGACGAGTTCGCATCTTgcttcttgcttcttgcttcttgcttctt                                                     <  1:450790/231‑1 (MQ=255)
                                                                                                                                                                        ggCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCGAAAAGGTGAGTCAGCGTCGCCATCGCAAACGGACGGTTCGCTGAGTCAAAACAGACGAGTTCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGGAt                                                 <  1:1033602/234‑1 (MQ=255)
                                                                                                                                                                         gcgcGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGt                                                                                                                                                                                                          >  2:1422014/1‑80 (MQ=34)
                                                                                                                                                                         gcgcGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGGCAAAAAGTGACTccc                                                                                                                                                                >  1:110772/1‑120 (MQ=17)
                                                                                                                                                                         gcgcGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGAcac                                                                                                                                                                  >  2:264052/1‑118 (MQ=17)
                                                                                                                                                                         gcgcGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCGAAAAgg                                                                                                                                                >  2:1322594/1‑138 (MQ=32)
                                                                                                                                                                         gcgcGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCGAAAAGGTGAGTCAGCGTCGCCATCGCAAACGGACGGTTCGCTGAGTCAAAACAGACGAGTTCGCAtcttgct                                                                              <  1:1008123/204‑1 (MQ=255)
                                                                                                                                                                            cgGTGCTGGTGGGAGCAAGGCCGGAGTGAGGGAAACCCGGGGCGCAGCGAAGGGGCCGGATGAGGGGAGCGCAGCGTTTTGTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCg                                                                                                                                                                 <  2:104800/118‑1 (MQ=11)
                                                                                                                                                                              gTCCTACCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTg                                                                                                                                                                      <  2:1187968/111‑1 (MQ=17)
                                                                                                                                                                                 cTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCGAAAAGGTGAGTCAGCGTCGCCATCGCAAACGGACGGTTCGCTGAGTCAAAACAGACGAGTTCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGGATGTGGGCATTc                                       <  1:1104197/235‑1 (MQ=255)
                                                                                                                                                                                 cTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCGAAAAGGTGAGTCAGCGTCGCCATCGCAAACGGACGGTTCGCTGAGTCAAAACAGACGAGTTCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGGATGTGGGCATTc                                       <  1:1512561/235‑1 (MQ=255)
                                                                                                                                                                                  tGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCGaca                                                                                                                                                   >  2:1108312/1‑124 (MQ=18)
                                                                                                                                                                                    cggggAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCGAAAAGGTGAGTCAGCGTCGCCATCGCAAACGGACGGTTCGCTGAGTCAAAACAGACGAGTTCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGGATGTGGGCATTca                                      <  1:488794/231‑1 (MQ=255)
                                                                                                                                                                                    ccGGGAGTAAGGCCGGAGGGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCGAAAAGGTGAGTCAGCGTCGCCATCGCAAACGGACGGTTCGCTGAGTCAAAACAGACGAGTTCGCAtcttgcttcttgc                                                                        <  1:481289/199‑1 (MQ=255)
                                                                                                                                                                                     cGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCGAAAAGGTGAGTCAGCGTCGCCATCGCAAACGGACGGTTCGCTGAGTCAAAACAGACGAGTTCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGGATGTGGGCATTcac                                     <  2:498865/233‑1 (MQ=255)
                                                                                                                                                                                     cGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCGAAAAGGTGAGTCAGCGTCGCCATCGCAAACGGACGGTTCGCTGAGTCAAAACAGACGAGTTCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGGATGTGGGCATTCACATc                                  <  1:221498/236‑1 (MQ=255)
                                                                                                                                                                                      gggAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCGAAAAGGTGAGTCAGCGTCGCCATCGCAAACGGACGGTTCGCTGAGTCAAAACAGACGAGTTCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGGATGTGGGCATTCACATc                                  <  1:772699/235‑1 (MQ=255)
                                                                                                                                                                                       ggggCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCGAAAAGGTGAGTCAGCGTCGCCATCGCAAACGGACGGTTCGCTGAGTCAAAACAGACGAGTTCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGGATGTGGGCAtt                                        <  1:384442/225‑1 (MQ=255)
                                                                                                                                                                                       ggAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCGAAAAGGTGAGTCAGCGTCGCCATCGCAAACGGACGGTTCGCTGAGTCAAAACAGACGAGTTCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGGATGTGGGCATTCACATCg                                 <  1:1160179/235‑1 (MQ=255)
                                                                                                                                                                                            aaGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGCCAAAAAGTGACTCGCCGTAAGGGCGAAAAggc                                                                                                                                               >  2:424613/1‑119 (MQ=17)
                                                                                                                                                                                            aaGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCGAAAAGGTGAGt                                                                                                                                           >  1:566021/1‑124 (MQ=33)
                                                                                                                                                                                            aaGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCGAAAAGGTGAGTCAGCGTCGCCATCGCAAACGGACGGTTCGCTGAGTCAAAACAGACGAGTTCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGGATGTGGGCATTcaca                                    <  1:1439308/227‑1 (MQ=255)
                                                                                                                                                                                                 cGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGGTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTa                                                                                                                                                            >  1:276330/1‑102 (MQ=17)
                                                                                                                                                                                                  ggAGTGAGGGAACCCCGGGGCGCGGCGCAGGGGCCGGAGGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCGAAAAGGTGAGTCAGCGTCGCCATCGCAAACGGACGGTTCGc                                                                                                             <  2:32397/148‑1 (MQ=255)
                                                                                                                                                                                                       gAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCGAAAAGGTGAGTCAGCGTCGCCATCGCAAACGGACGGTTCGCTGAGTCAAAACAGACGAGTTCGCAtcttgcttcttgcttcttgcttc                                                              >  1:1406309/1‑189 (MQ=255)
                                                                                                                                                                                                       gAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCGAAAAGGTGAGTCAGCGTCGCCATCGCAAACGGACGGTTCGCTGAGTCAAAACAGACGAGTTCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTAttcattca                 <  1:1374954/235‑1 (MQ=255)
                                                                                                                                                                                                        aGGGAAACCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCGAAAAGGTGAGTCAGCGTCGCCATCGCAAACGGACGGTTCGCTGAGTCAAAACAGACGAGTTCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTgg                                                   <  2:1036981/200‑1 (MQ=255)
                                                                                                                                                                                                          ggAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCGAAAAGGTGAGTCAGCGTCGCCATCGCAAACGGACGGTTCGCTGAGTCAAAACAGACGAGTTCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGGATGTGGGCATTCACATCGAg                               <  2:572880/218‑1 (MQ=255)
                                                                                                                                                                                                          ggAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCGAAAAGGTGAGTCAGCGTCGCCATCGCAAACGGACGGTTCGCTGAGTCAAAACAGACGAGTTCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCACg               <  2:1306098/234‑1 (MQ=255)
                                                                                                                                                                                                                ccGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCGAAAAGGTGAGTCAGCGTCGCCATCGCAAACGGACGGTTCGCTGAGTCAAAACAGACGAGTTCGCAtcttgcttcttgcttcttgcttc                                                              <  2:1371534/181‑1 (MQ=255)
                                                                                                                                                                                                                 cGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCGAAAAGGTGAGTCAGCGTCGCCATCGCAAACGGACGGTTGGCTGAGTCAAAACAGACGAGTTCGCAtcttgcgtcttgcttcttgcttcttgct                                                         >  1:381993/1‑185 (MQ=255)
                                                                                                                                                                                                                 cGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCGAAAAGGTGAGTCAGCGTCGCCATCGCAAACGGACGGTTCGCTGAGTCAAAACAGACGAGTTCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGGAt                                                 <  2:140008/193‑1 (MQ=255)
                                                                                                                                                                                                                 cGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTCTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCAGTAAGGGCGAAGAGGTGAGTCAGCGTc                                                                                                                                    >  2:1383696/1‑110 (MQ=255)
                                                                                                                                                                                                                  ggAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTGTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCGAAAAGGTGAGTCAGCGTCGCCATCGCAAACGGACGGTTCGCTGAGTCAAAACAGACGAGTTCGCAtcttgcttcttgcttcttgcttc                                                              <  2:407659/179‑1 (MQ=255)
                                                                                                                                                                                                                       gcagcgcagGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCGAAAAGGTGAGTCAGCGTCGCCATCGCAAACGGACGGTTCGCTGAGTCAAAACAGACGAGTTCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCACGCAGGTGGCGCATa  >  2:899806/1‑234 (MQ=255)
                                                                                                                                                                                                                        cagcgcagGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCGAAAAGGTGAGTCAGCGTCGCCATCGCAAACGGACGGTTCGCTGAGTCAAAACAGACGAGTTCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGGATGTGGGCATTCACATCGAGTAg                            <  1:147339/207‑1 (MQ=255)
                                                                                                                                                                                                                        |                                                                                                                                                                                                                                        
AACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCGAAAAGGTGAGTCAGCGTCGCCATCGCAAACGGACGGTTCGCTGAGTCAAAACAGACGAGTTCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCACGCAGGTGGCGCATA  >  AE015451/4022178‑4022626

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