Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 4,635,590 T→C 46.3% K56E (AAA→GAA)  uvrY ← two‑componenent system BarA/UvrY ‑ regulatory subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154514,635,5900TC46.3% 9.5 / 77.7 54K56E (AAA→GAA) uvrYtwo‑componenent system BarA/UvrY ‑ regulatory subunit
Reads supporting (aligned to +/- strand):  ref base T (16/13);  new base C (13/12);  total (29/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.50e-01

GCGGGCTGATATAGCGCTGGCCGGCAAAGGCCAGACGAATGGCCTGGACCATTTCGTCAAGGCCTGCACCCTTGGTCAGGTAGCCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATACCGGTTCGTACCAGATCGTGATCATCGACCACTAAGACCCTAATCAAGCACACACCTCGTCAGC  >  AE015451/4635374‑4635772
                                                                                                                                                                                                                        |                                                                                                                                                                                      
gCGGGCTGATATAGCGCTGGCCGGCAAAGGCCAGACGAATGGCCTGGACCATTTCGTCAAGGCCTGCACCCTTGGTCAGGTAGCCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACAGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTCCACATCCATCAGcacac                                                                                                                                                                       >  2:813570/1‑233 (MQ=255)
                     cGGCAAAGGCCAGACGAATGGCCTGGACCATTTCGTCAAGGCCTGCACCCTTGGTCAGGTAGCCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCggg                                                                                                                                                  >  2:450601/1‑234 (MQ=255)
                      ggCAAAGGCCAGACGAATGGCCTGGACGATTTAGTCAAGGCCTGCACCCTTGGTCAGGTAGCCGGCAGCGCCGGCCTGCAACAAGCGTGTGGAGAACGGGTCTTCCTCGCACACCGTCACCGCCACTCCCTTCATGTCAGGGTGGCTGCGCAGCAGTGTGCTGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGcacatc                                                                                                                                                                      >  2:577446/1‑211 (MQ=255)
                      ggCAAAGGCCAGACGAATGGCCTGGACCATTTCGTCAAGGCCTGCACCCTTGGTCAGGTAGCCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAg                                                                                                                                                                            >  2:805003/1‑207 (MQ=255)
                      ggCAAAGGCCAGACGAATGGCCTGGACCATTTCGTCAAGGCCTGCACCCTTGGTCAGGTAGCCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTCCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCggg                                                                                                                                                  >  2:474579/1‑233 (MQ=255)
                                   cGAATGGCCTGGACCATTTCGTCAAGGCCTGCACCCTTGGTCAGGTAGCCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCg                                                                                                                                   <  1:862620/235‑1 (MQ=255)
                                    gAATGGCCTGGACCATTTCGTCAAGGCCTGCACCCTTGGTCAGGTAGCCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTCCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGCTTCAGTGCCgc                                                                                                                                  >  1:304993/1‑234 (MQ=255)
                                          ccTGGACCATTTCGTCAAGGCCTGCACCCTTGGTCAGGTAGCCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGGTTCAGCTCCAggg                                                                                                                                                  >  2:108376/1‑213 (MQ=255)
                                             ggACCATTTCGTCAAGGCCTGCACCCTTGGTCAGGTAGCCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTTc                                                                                                                                                                                       <  1:698495/173‑1 (MQ=255)
                                                      cGTCAAGGCCTGCACCCTTGGTCAGGTAGCCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTca                                                                                                               >  1:417945/1‑236 (MQ=255)
                                                      cGTCAAGGCCTGCACCCTTGGTCAGGTAGCCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTCCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTca                                                                                                               >  1:533393/1‑236 (MQ=255)
                                                        tCAAGGCCTGCACCCTTGGTCAGGTAGCCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTcacc                                                                                                             <  1:450604/236‑1 (MQ=255)
                                                                  cACCCTTGGTCAGGTAGCCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTCCACATCCATCAGCACAACGTCCGGtt                                                                                                                                                              <  1:533851/177‑1 (MQ=255)
                                                                   aCCCTTGGTCAGGTAGCCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTCCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCt                                                                                                     >  2:902402/1‑233 (MQ=255)
                                                                     ccTTGGTCAGGTAGCCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCa                                                                                                                                               >  2:684532/1‑189 (MQ=255)
                                                                     ccTTGGTCAGGTAGCCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACccacca                                                                                                        >  2:907472/1‑228 (MQ=255)
                                                                     ccTTGGTCAGGTAGCCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTCCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTcaccc                                                                                                            >  1:254886/1‑224 (MQ=255)
                                                                          gTCAGGTAGCCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTCCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGt                                                                                                                                             >  2:341455/1‑186 (MQ=255)
                                                                          gTCAGGTAGCCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTCCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCCGcc                                                                                                  >  2:843298/1‑228 (MQ=255)
                                                                                      gCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGc                                                                                                                             >  2:599875/1‑190 (MQ=255)
                                                                                                     aacaaGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATAcc                                                                <  1:446055/236‑1 (MQ=255)
                                                                                                      acaaGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATAcc                                                                <  1:907475/235‑1 (MQ=255)
                                                                                                       caaGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTCCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATccc                                                                >  2:643005/1‑234 (MQ=255)
                                                                                                            gtgtGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCAAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCac                                                                               >  2:582002/1‑213 (MQ=255)
                                                                                                                  ggAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTCCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTca                                                                                                               <  2:529965/176‑1 (MQ=255)
                                                                                                                        ggTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATACCGGTTCGTACCAg                                                    >  2:619243/1‑229 (MQ=255)
                                                                                                                            ttCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTCCACATCCATCCGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGc                                                                                                                             >  2:85421/1‑152 (MQ=255)
                                                                                                                               cccGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTCCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATACCGGTTCGTa                                                        <  1:703492/216‑1 (MQ=255)
                                                                                                                                tCGAACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTc                                                                                           <  1:857058/182‑1 (MQ=255)
                                                                                                                                 cGCACCCCGCCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTCCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAAtctc                                                                                                                    <  1:181781/156‑1 (MQ=255)
                                                                                                                                  gCACACTGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCCTTTCCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATACCGGTTCGTACCAGATCGTGATCATCGAcc                                    <  1:886667/235‑1 (MQ=255)
                                                                                                                                   cacaCCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTCCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATACCGGTTCGt                                                         <  1:927207/213‑1 (MQ=255)
                                                                                                                                   cacaCCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTCCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATACCGGTTCGTACCAGATCGTGATCATCGACCAc                                  >  1:250054/1‑236 (MQ=255)
                                                                                                                                        aggcACCGCCACGACCTGGATGGCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGGCTCCAGGCCGCCGATGCCGGGCATTTTc                                                                                                                                                                                       <  2:53730/79‑1 (MQ=255)
                                                                                                                                          tCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTCCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATACCGGTTCGTACCAGATCGTGATCATCGACCa                                   <  1:630790/228‑1 (MQ=255)
                                                                                                                                           cACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTCCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATACCGGTTCGTACCAGATCGTGATCATCGACCACTAAGAcc                           >  1:380066/1‑235 (MQ=255)
                                                                                                                                             ccGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATACCGGTTCGTACCAGATCGTGATCATCGACCACTAAGACCCt                         >  1:415101/1‑235 (MQ=255)
                                                                                                                                                   ttACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTCCACATCCATCAGCACACCGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTcaccc                                                                                                            <  2:451639/145‑1 (MQ=255)
                                                                                                                                                               cAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATACCGGTTCGTACCAGATCGTGATCATCGACCACTAAGACCCTAATCAAGCACACACCTCGt      >  1:896521/1‑236 (MQ=255)
                                                                                                                                                                agcgTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTCCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATACCGGTTCGTACCAGATCGTGatca                                          <  1:913801/196‑1 (MQ=255)
                                                                                                                                                                aGGGTGGCTGCGCAGCAGTTGGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATACCGGTTCGTACCAGATCGTGATCATCGACCACTAAGACCCTAATCAAGca                <  1:272896/225‑1 (MQ=255)
                                                                                                                                                                 gggTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATACCGGTTCGTACCAGATCGTGATCATCGACCACTAAGACCCTAATCAAGCACACACCTCGTc     <  1:545817/235‑1 (MQ=255)
                                                                                                                                                                   gTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATACCGGTTCGTACCAGATCGTGATCATCGACCACTAAGACCCTAATCAAGCACACACCTc        >  2:588902/1‑230 (MQ=255)
                                                                                                                                                                   gTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATACCGGTTCGTACCAGATCGTGATCATCGACCACTAAGACCCTAATCAAGCACACACCTCGTcagc  >  1:145062/1‑236 (MQ=255)
                                                                                                                                                                   gTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATACCGGTTCGTACCAGATCGTGATCATCGACCACTAAGACCCTAATCAAGCACACACCTCGTcagc  <  1:824472/236‑1 (MQ=255)
                                                                                                                                                                    tGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTCCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATACCGGTTCGTACCAGATCGTGatcat                                         >  2:846673/1‑196 (MQ=255)
                                                                                                                                                                       cTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATACCGGTTCGTACCAGATCGTGATCATCGACCACTa                                >  1:573796/1‑202 (MQ=255)
                                                                                                                                                                          cgcAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTCCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATAcc                                                                <  1:177635/167‑1 (MQ=255)
                                                                                                                                                                                  cttGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCCCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATACCGGTTCGTACCAGATCGTGATCATCGACCACTAAGACCCTAATCAAGca                <  2:145061/206‑1 (MQ=255)
                                                                                                                                                                                        ggTGGCCTCCAGGCCGCCGATGCCGGGCATTTCCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTcacc                                                                                                             <  2:118919/108‑1 (MQ=255)
                                                                                                                                                                                          tGGCCTCCAGGCCGCCGATGCCGGGCATTTCCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGTAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATACCGGTTCGTACCAGATCGTGATCATCGACCAc                                  <  1:813573/181‑1 (MQ=255)
                                                                                                                                                                                                       gccgATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTcacc                                                                                                             <  2:264053/93‑1 (MQ=255)
                                                                                                                                                                                                           aTGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATcc                                                                          <  2:633301/124‑1 (MQ=255)
                                                                                                                                                                                                           aTGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATcc                                                                          >  1:633304/1‑124 (MQ=255)
                                                                                                                                                                                                                      tttCACACCCATCAGCTCAACGTCCGCTTTCAGCTCCCGGGCCAGTTTCAGGGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATAcc                                                                <  2:119837/123‑1 (MQ=255)
                                                                                                                                                                                                                        |                                                                                                                                                                                      
GCGGGCTGATATAGCGCTGGCCGGCAAAGGCCAGACGAATGGCCTGGACCATTTCGTCAAGGCCTGCACCCTTGGTCAGGTAGCCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATACCGGTTCGTACCAGATCGTGATCATCGACCACTAAGACCCTAATCAAGCACACACCTCGTCAGC  >  AE015451/4635374‑4635772

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