Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 520,534 A→G 24.8% E40G (GAG→GGG)  ybbP → putative ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913520,5340AG24.8% 101.1 / 5.3 56E40G (GAG→GGG) ybbPputative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (8/32);  new base G (13/2);  total (22/34)
Fisher's exact test for biased strand distribution p-value = 9.93e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTTTTCCCTCAACCGTGAACATGGCACCACGTTGATTATGGTGACCCACGACCTGCAACTGGCGGCACGCTGCGACCGCTGCTTACGGCTGGTGAACGGGCAGTTGCAGGAGGAAGCATGATTGCACGTTGGTTCTGGCGCGAATGGCGTTCGCCGTCGCTATTAATTGTCTGGCTGGCGCTAAGCCTGGCGGTGGCCTGCGTGCTGGCGCTGGGCAATATCAGCGATCGCA‑TGGAGAAGGGCTTAAGCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGTTCACGCGAAGTGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAGCGCGGCCTGAAAGTCGGCAAGCAGCTGACTTTCG‑CCACAATGACCTTTGCAGGCGACACACCGCAGCTGGCGAACGTCAAAGCGGTGGATGATATCTACCCGATGTATGGCGATCTGCAAACTAATCCCCCTGGCCT  >  NC_000913/520298‑520771
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                              
tgTTTTCCCTCAACCGTGAACATGGCCCCCGGTTGATTATGGTGACCCACGACCTGCAACGGGCGGCACGCTGCGACCGCTGCTTACGGCTGGTGAACGGGCAGTTGCAGGAGGAAGCATGATTGCACGTTGGTTCTGGCGCGAATGGCGTTCGCCGTCGCTATTAATTGTCTGGCTGGCGCTAAGCCTGGCGGTGGCCTGCGTGCTGGCGCTGGGCAATATCAGCGATCGCA‑TGgaga                                                                                                                                                                                                                                              <  1:78032/239‑1 (MQ=255)
      ccccCAAACGTGACCATGGCACCCCGTTGATTATGGTGACCCACCGCCTGCAACTGGCGGCACGCTGCGACCGCTGCTTACGGCTGGTGAACGGGCAGTTGCAGGAGGAAGCATGATTTCACGTTGGTTCTGGCGCGAATGGCGTTCGCCGTCGCTATTAATTGTCTGGCTGGCGCTAAGCCTGGCGGTGGCCTGCGTGCTGGCGCTGGGCAATATCAGCGATCGCA‑TGGAGAAGGGCt                                                                                                                                                                                                                                        <  1:602222/239‑1 (MQ=255)
         tCAACCGTGAACATGGCACCACGTTGATTATGGTGACCCACGACCTGCAACTGGCGGCACGCTGCGACCGCTGCTTACGGCTGGTGAACGGGCAGTTGCAGGAGGAAGCATGATTGCACGTTGGTTCTGGCGCGAATGGCGTTCGCCGTCCCTATTAATTGTCTGGCTGGCGCTAAGCCTGGCGGTGGCCTGCGGGCTGGCGCTGGGCAATATCAGCGATCGCA‑TGGGGACGGGCttt                                                                                                                                                                                                                                      >  1:659607/1‑237 (MQ=255)
          cAACCGTGAACATGGCACCACGTTGATTATGGTGACCCACGACCTGCAACTGGCGGCACGCTGCGACCGCTGCTTACGGCTGGTGAAAGGGCAGTTGCAGGAGGAAGCATGATGGCACGTTGGTTCTGGCGCGAATGGCGTTCGCCGTCGCTATAAATTTTCTGGCTGGCCCTAAGCCTGGCGGTGGCCTGCGTGCTGGCGCTGGGCAATATCAGCAAGCACA‑AGGGGAAGGGCTAAAg                                                                                                                                                                                                                                    >  2:145554/1‑239 (MQ=255)
                tGAACATGGCACCACGTTGATTATGGTGACCCACGACCTGCAACTGGCGGCACGCTGCGACCGCTGCTTACGGCTGGTGAACGGGCAGTTGCAGGAGGAAGCATGATTGCACGGTGGTTCTGGCGCGAATGGCGTTCGCCGTCGCTATTAATTGTCTGGCTGGCGCTAAGCCTGGCGGTGGGCTGCGTGCTTGGACTGGGGAAATACAGCGATCGCA‑TGGGGAAGGGGTTAAGCCAgac                                                                                                                                                                                                                              >  2:684134/1‑237 (MQ=255)
                 gAACATGGCACCACGTTGATTATGGTTACCCACGACCTGCAACTGGCGGCACGCTGCGACCGCTGCTTACGGCTGGTGAACGGGCAGTTGCAGGAGGAAGCATGATTGCACGTTGGTTCTGGCGCGAATGGCGTTCGCCGTCGCTATTAATTGTCTGGCTGGCGCTAAGCCTGGCGGTGGCCTGCGTGCTGGCGCTGGGCAATATCAGCGATCGCA‑TGGAGAAGGGCTTAAGCCAGCa                                                                                                                                                                                                                              <  1:768049/238‑1 (MQ=255)
                       ggCACCACGTTGATTATGGTGACCCACGACCTGCAACTGGCGGCACGCTGCGACCGCTGCTTACGGCTGGTGAACGGGCAGTTGCAGGAGGAAGCATGATTGCACGTTGGTTCTGGCGCGAATGGCGTTCGCCGTCGCTATTAATTTTCTGGCTGGCGCTAAGCCTGGCGGGGGCCTGCGTGCTGGGGCTGGGCCATATCCGCGATCGCA‑TGGGGAAGGGCTTAAGCCAGCaaaagcgg                                                                                                                                                                                                                       >  2:96699/1‑234 (MQ=255)
                       ggCACCACGTTGATTATGGTGACCCACGACCTGCAACTGGCGGCACGCTGCGACCGCTGCTTACGGCTGGTGAACGGGCAGTTGCAGGAGGAAGCATGATTGCACGTTGGTTCTGGCGCGAATGGCGTTCGCCGTCGCTATTAATTGTCTGGCTGGCGCTAAGCCTGGCGGTGGCCTGCGGGCTGGCGCTGGGCAATATCAGGGATCGCA‑TGGCGAAGGGGTTAAGCCCGCAAAGCCGt                                                                                                                                                                                                                       >  1:480769/1‑239 (MQ=255)
                         ccccccGTTGATTATGGTGACCCCCGGCCTGCAACTGGCGGCCCGCTGCGACCGCTGCTTACGGCTGGGGAACGGGCAGTTGCAGGAGGAAGCATGATTGCACGTTGGTTCTGGCGCGAATGGCGTTCGCCGTCGCTATTAATTGTCTGGCTGGCGCTAAGCCTGGCGGTGGCCTGCGTGCTGGCGCTGGGCAATATCAGCGATCGCA‑TGGAGAAGGGCTTAAGCCAGCAAAGCCg                                                                                                                                                                                                                        <  2:688409/231‑1 (MQ=255)
                              cGTTGATTATGGTGACCCACGACCTGCAACTGGCGGCACGCTGCGACCGCTGCTTACGGCTGGTGAACGGGCAGTTGCAGGAGGAAGCATGATTGCACGTTGGTTCTGGCGCGAATGGCGTTCGCCGTCGCTATTAATTGTCTGGCTGGCGCTAAGCCTGGCGGTGGCCTGCGTGCTGGCGCTGGGCAATATCAGCGATCGCA‑TGGGGAAGGGCTTAAGCCAGCAAAGCCGTGAGTTTa                                                                                                                                                                                                                >  1:894890/1‑239 (MQ=255)
                                 tGATTATGGTGACCCACGACCTGCAACTGGCGGCACGCTGCGACCGCTGCTTACGGCTGGTGCACGGGCAGTTGCAGGGGGAAGCATGATTGCACGTTGGTTCTGGCGCGAATGGCCTTCGCCGTCGCTATTAATTTTCTGGCTGGCGCTAAGCCTGGCGGTGGCCTGCGGGCTGGCCGTGGGCGATATCAGCGCTCGCA‑GGGGGAAGGGGTTAAGCCCGCaaaagccggggtttaagg                                                                                                                                                                                                             >  2:788725/1‑224 (MQ=255)
                                  gATTATGGTGACCCACGACCTGCAACTGGCGGCACGCTGCGACCGCTGCTTACGGCTGGTGAACGGGCAGTTGCAGGAGGAAGCATGATTGCACGTTGGTTCTGGCGCGAATGGCGTTCGCCGTCGCTATTAATTGTCTGGCTGGCGCTAAGCCTGGCGGTGGCCTGCGTGCTGGCGCTGGGCAATATCAGCGATCGCA‑TGGAGAAGGGCTTAAGCCAGCAAAGCCGTGAGTTTATggc                                                                                                                                                                                                            <  1:276603/239‑1 (MQ=255)
                                   aTTATGGTGACCCACGACCTGCAACTGGCGGCACGCTGCGACCGCTGCTTACGGCTGGTGAACGGGCAGTTGCAGGAGGAAGCATGATTGCACGTTGGTTCTGGCGCGAATGGCGTTCGCCGTCGCTATTAATTGTCTGGCTGGCGCTAAGCCTGGCGGTGGCCTGCGTGCTGGGGCTGGGCAATATCAGCGATCGCA‑TGGGGAAGGGCTTAAGCCAGCCAAGCCGGGAGTTTATg                                                                                                                                                                                                              >  1:688409/1‑236 (MQ=255)
                                   aTTATGGTGACCCACGACCTGCAACTGGCGGCACGCTGCGACCGCTGCTTACGGCTGGTGAACGGGCAGTTGCAGGAGGAAGCATGATTGCACGTTGGTTCTGGCGCGAATGGCGTTCGCCGTCGCTATTAATTGTCTGGCTGGCGCTAAGCCTGGCGGTGGCCTGCGTGCTGGCGCTGGGCAATATCAGCGATTGCA‑TGGAGAAGGGCTTAAGCCAGCAAAGCCGTGAGTTTATggcg                                                                                                                                                                                                           >  1:851442/1‑239 (MQ=255)
                                   aTTATGGTGACCCACGACCTGCAACTGGCGGCACGCTGCGACCGCTGCTTACGGCTGGTGAACGGGCAGTTGCAGGAGGAAGCATGATTGCACGTTGGTTCTGGCGCGAATGGCGTTCGCCGTCGCTATTAATTGTCTGGCTGGCGCTAAGCCTGGCGGTGGCCTGCGTGCTGGCGCTGGGCAATATCAGCGATCGCA‑TGGAGAAGGGCTTAAGc                                                                                                                                                                                                                                   <  1:231610/215‑1 (MQ=255)
                                    ttttggtccccccCGCCCTGCAACGGGCGGGCCCCTGCGGCCGCTGCTGACGGGTGGGGAACGGGCAGTTGCAGGAGGAAGCATGATTGCACGTTGGTTCGGGCGCGAATGGCTTTCGCCGTCGCTATTAATTGTCTGGCTGGCGCTAAGCCTGGCGGTGGCCTGCGTGCTGGCGCTGGGCAATATCAGCGATCGCA‑TGGAGAAGGGCTTAAGCCAGCAAAGCCGTGAGTTTATggcg                                                                                                                                                                                                           <  2:737008/229‑1 (MQ=255)
                                       tGGTGACCCACGACCTGCAACTGGCGGCACGCTGCGACCGCTGCTTACGGCTGGTGAACGGGCAGTTGCAGGGGGAAGCATGATTGCACGTTGGTTCTGGCGCGAATGGCGTTCGCCGTCGCTATTAATTGTCTGGCTGGCGCTAAACCTGGCGGGGGCCGTCGGTGTTGCCCTGGGCAAAATGCGCGATCGCC‑TGGGGAAGGGGGTAAGCCCGCAAACACGGGGGTTTagtggggggg                                                                                                                                                                                                       >  2:276603/1‑230 (MQ=255)
                                             cccACGACCTGCAACTGGCGGCACGCTGCGACCGCTGCTTACGGCTGGTGAACGGGCAGTTGCAGGAGGAAGCATGATTGCACGTTGGTTCTGGCGCGAATGGCGTTCGCCGTCGCTATTAATTGTCTGGCTGGCGCTAAGCCTGGCGGTGGCCTCCGTGCTGGCGCTGGGCACAATACGCGATTCCA‑TGGGGAAGGGCTTTAGCCACCAAAACCCGGGGTTTTTgggggggggggggg                                                                                                                                                                                                 >  2:911200/1‑231 (MQ=255)
                                             cccACGACCTGCAACTGGCGGCACGCTGCGACCGCTGCTTACGGCTGGTGAACGGGCAGTTGCAGGAGGAAGCATGATTGCACGTTGGTTCTGGCGCGAATGGCGTTCGCCGTCGCTATTAATTGTCTGGCTGGCGCTAAGCCTGGCGGGGGGCTCCGTGATGGCGCTGGGCAAAAATAGCGGTCCCAGTGGA‑AAGGGCTTAACCCAGCAAAACCGTGGGGTTCTggggcgggggcagg                                                                                                                                                                                                 >  2:748768/1‑227 (MQ=255)
                                             cccACGACCTGCAACTGGCGGCACGCTGCGACCGCTGCTTACGGCTGGTGAACGGGCAGTTGCAGGAGGAAGCATGATTGCACGTTGGTTCTGGCGCGAATGGCGTTCGCCGTCGCTATTAATTGTCTGGCTGGCGCTAAGCCTGGCGGGGGCCTGCGGGCGGGCGCTGGGCAATATCCGCGATCGCA‑TGGGGAAGGGCTTAAGCCAGCAAAACCGGGGGTTTATggggggggggcggg                                                                                                                                                                                                 >  2:847089/1‑231 (MQ=255)
                                             cccACGACCTGCAACTGGCGGCACGCTGCGACCGCTGCTTACGGCTGGTGAACGGGCAGTTGCAGGAGGAAGCATGATTGCACGGTGGTTCTGGCGCGAATGGCGTTCGCCGTCGCTATTAATTGTCTGGCTGGCGCTCAGACTGGCGGTGGCCTGCGTGCTGGCGCTGGGCCATATCAGCGAACGCA‑TGGAGAAGGGCTTAAACCAGCACAACCgccggcccttatacccaactcac                                                                                                                                                                                                  >  2:231610/1‑216 (MQ=255)
                                                     cTGCAACTGGCGGCACGCTGCGACCGCTGCTTACGGCTGGTGAACGGGCAGTTGCAGGAGGAAGCATGATTGCACGTTGGTTCTGGCGCGAATGGCGTTCGCCGTCGCTATTAATTGTCTGGCTGGCGCTAAGCCTGGCGGTGGCCTGCGTGCTGGCGCTGGGCAATATCAGCGATCGCA‑TGGAGAAGGGCTTAAGCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTgcg                                                                                                                                                                                          <  1:781401/238‑1 (MQ=255)
                                                              gcggcACGCTGCGACCGCTGCTTACGGCTGGTGAACGGGCAGTTGCAGGAGGAAGCATGATTGCACGTTGGTTCTGGCGCGAATGGCGTTCGCCGTCGCTATTAATTGTCTGGCTGGCGCTAAGCCTGGCGGTGGCCTGCGTGCTGGCGCTGGGCAATATCAGCGATCGCA‑TGGGGAAGGGCTTAAGCCAGCAAAGCCGTGGGTTTATggcgggcggcaggggggtggcccgtttccc                                                                                                                                                                                 >  1:453973/1‑216 (MQ=255)
                                                                             cGCTTCTTCCGGGTGGGGAACGGGCAGTTGCAGGAGGAAGTATTATTGCACGTTGGTTCTGGCGCGAATGGCGTTCGCCGTCGCTATTAATTGTCTGGCTGGCGCTAAGCCTGGCGGTGGCCTGCGTGCTGGCGCTGGGCAATATCAGCGATCGCA‑TGGAGAAGGGCTTAAGCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGTTCACGCGAAGTGCCGCAAGCg                                                                                                                                                                 <  2:1044867/239‑1 (MQ=255)
                                                                             cGCTGTTTACGGCTGGTGAACGGGCAGTTGCAGGAGGAAGCATGATTGCACTTTGGTTCTGGCGCGAATGGCGTTCGCCGTCGCTATTAATTGTCTGGCTGGCGCTAAGCCTGGCGGTGGCCTGCGTGCTGGCGCTGGGCAATATCAGCGATCGCA‑TGGAGAAGGGCTTAAGCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGTTCACGCGAAGTGCCGCAAGCg                                                                                                                                                                 <  2:833376/239‑1 (MQ=255)
                                                                             cGCTGGTTACGGCTGGTGAACGGGCAGTTGCAGGAGGAAGCATGATTGCACGTTGGTTCTGGCGCGAATGGCGTTCGCCGTCGCTATTAATTGTCTGGCTGGCGCTAAGCCTGGCGGTGGCCTGCGTGCTGGCGCTGGGCAATATCAGCGATCGCA‑TGGAGAAGGGCTTAAGCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGTTCACGCGAAGTGCCGCAAGCg                                                                                                                                                                 <  1:96699/239‑1 (MQ=255)
                                                                               ctgctTACGGCTGGTGAACGGGCAGTTGCAGGAGGAAGCATGATTGCACGTTGGTTCTGGCGCGAATGGCGTTCGCCGTCGCTATTAATTGTCTGGCTGGCGCTAAGCCTGGCGGTGGCCTGCGGGCTGGCGCTGGGCAATGTCCGCGAGCGCA‑TGGGGAAGGGCTTAAGCCACACAGGCCGGGGGTTGTTGGCGGGCGGTCGGGGGGGGCGCGGCTGACGCGGAGGGGGGCAAggggg                                                                                                                                                               >  2:1078776/1‑235 (MQ=255)
                                                                                 ggtttgggCTGGTGACGGGGCAGGTGCAGGAGGAGGCTTGATTGCACGTTGGTTCTGGCGCGAATGGCGTTCGCCGTCGCTATTAATTGTCTGGCTGGCGCTAAGCCTGGCGGTGGCCTGCGTGCTGGCGCTGGGCAATATCAGCGATCGCA‑TGGAGAAGGGCTTAAGCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGTTCACGCGAAGTGCCGCAAGCGTGGc                                                                                                                                                             <  2:851442/233‑1 (MQ=255)
                                                                                        gCTGGTGAACGGGCAGTTGGAGGAGGAAGCATGAGTGCACGGTGGTTCTGGCGCGAATGGCGGTCGCCGTCGCTATTAATTGTCTGGCTGGCGCTAAGCCTGGCGGTGGCCTGCGTGCTGGCGCTGGGCAATATCAGCGATCGCA‑TGGAGAAGGGCTTAAGCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGTTCACGCGAAGTGCCGCAAGCGTGGCTggagg                                                                                                                                                       <  1:374173/238‑1 (MQ=255)
                                                                                                     cAGTTGCAGGAGGAAGCATGATTGCACGTTGGTTCTGGCGCGAATGGCGTTCGCCGGCGCTATTAATTGTCTGGCTGGCGCTAAGCCTGGCGGTGGCCTGCGTGCTGGCGCTGGGCAATATCAGCGATCGCA‑TGGAGAAGGGCTTAAGCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGTTCACGCGAAGTGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAgcg                                                                                                                                          <  1:263997/238‑1 (MQ=255)
                                                                                                      cgTTGCAGGAGGAAGCATGATTGCACGTTGGTTCTGGCGCGAATGGCGTTCGCCGTCGCTATTAATTGTCTGGCTGGCGCTAAGCCTGGCGGTGGCCTGCGTGCTGGCGCTGGGCAATATCAGCGATCGCA‑TGGAGAAGGGCTTAAGCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGTTCACGCGAAGTGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAgcgc                                                                                                                                         <  1:939930/237‑1 (MQ=255)
                                                                                                                   aGCATGATTGCACGTTGGTTCTGGCGCGAATGGCGTTCGCCGTCGCTATTAATTGTCTGGCTGGCGCTAAGCCTGGCGGTGGCCTGCGTGCTGGCGCTGGGCAATATCAGCGATCGCA‑TGGAGAAGGGCTTAAGCCAGCAAAGCCGTGAGTTTATGGCGGGGGATCGGGGGGTGGGCAGGTTACGCGGAAGGCCGCCAGCGTGGCTGGGGGGAAGGCACAAACGCGGCCCGaaaagcg                                                                                                                            >  1:214455/1‑234 (MQ=255)
                                                                                                                      attgttGCCCGTTGGTTCGGGCGCGAATGGCGTTCGCCGTCGCTATTAATTGTCTGGCTGGCGCTAAGCCTGGCGGTGGCCTGCGTGCTGGCGCTGGGCAATATCAGCGATCGCA‑TGGAGAAGGGCTTAAGCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGTTCACGCGAAGTGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAGCGCGGCCTGAAAGTCGGCa                                                                                                                         <  2:894890/234‑1 (MQ=255)
                                                                                                                         aTTGCACGTTGGTTCTGGCGCGAATGGCGTTCGCCGTCGCTATTAATTGTCTGGCTGGCGCTAAGCCTGGCGGTGGCCTGCGTGCTGGCGCTGGGCAATATCAGCGATCGCA‑TGGGGAAGGGCTTAAGCCAGCAAAGCCGTGAGTTTATGGCGGGGGATCGGGGGGGTGGCGGGTCCCCCGAAGGGGCGCAAGCGGGGGTGGGGGAAGGGCAAAAACGCGGGCCGAAAAAAGGCAagag                                                                                                                     >  2:133670/1‑237 (MQ=255)
                                                                                                                           tGCACGTTGGTTCTGGCGCGAATGGCGTTCGCCGTCGCTATTAATTGTCTGGCTGGCGCTAAGCCTGGCGGTGGCCTGCGTGCTGGCGCTGGGCAATATCAGCGATCGCA‑TGGAGAAGGGCTTAAGCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGTTCACGCGAAGTGCCGCAAGCg                                                                                                                                                                 <  1:705632/193‑1 (MQ=255)
                                                                                                                             cACGTTGGTTCTGGCGCGAATGGCGTTCGCCGTCGCTATTAATTGTCTGGCTGGCGCTAAGCCTGGCGGTGGCCTGCGTGCTGGCGCTGGGCAATATCAGCGATCGCA‑TGGAGAAGGGCTTAAGCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGTTCACCCGAAGTGCCGCCAGCGTGGCTGGAGGAAGCGCCAAAACGCCGCCCGAAAATCGGCCagcagca                                                                                                                  >  1:571963/1‑237 (MQ=255)
                                                                                                                                         ggCGCATATTGCGTTGGCCGTCTCTATTTATTGTCTGGCTGGCGTTTAGCCTGGGGGTGGCCTGGCTTGTGGCGCTGGGCAATTTCAGCGATCGCA‑TGGAGAAGGGCTTAAGCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGTTCCCGCGAAGTGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAGCGCGGCCTGAAAGTCGGCAAGCAGCTGACTTTCG‑Ccaca                                                                                                    <  2:458471/239‑1 (MQ=255)
                                                                                                                                                          ccGTCGCTATTAATTGTCTGGCTGGCGCTAAGCCTGGCGGTGGCCTGCGTGCTGGCGCTGGGCAATATCAGCGATCGCA‑TGGAGAAGGGCTTAAGCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCCGTTCACACGAAGTGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAACGCGGCCTGAAAGTCGGCAAGCCACTGACTTTTG‑CCACACAGACCctttaccgcg                                                                                    >  1:254531/1‑231 (MQ=255)
                                                                                                                                                                  tataaggggctggcgtgggcTAAGCGTGGCGGTGGCCTGCGTGCTGGCGCTGGGCAATATCAGCGATCGCA‑TGGAGAAGGGCTTAAGCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGTTCACGCGAAGTGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAGCGCGGCCTGAAAGTCGGCAAGCAGCTGACTTTCG‑CCACAATGACCTTTGCAGGCGACACACCg                                                                            <  1:448791/220‑1 (MQ=255)
                                                                                                                                                                   ttAATTGTCTGGCTGGCGCTAAGCCTGGCGGTGGCCTGCGTGCTGGCGCTGGGCAATATCAGCGATCGCA‑TGGAGAAGGGCTTAAGCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGGTCACGCGAAGTGCCGCAAGCGTGGCTGGGGGAAGCGCAAAAACGCGGCCTGAAAGTCGGCAAGCAGCTGGCTTTTGCCCACAA‑GACCTTTTCAGGGGACACAccccc                                                                          >  1:788139/1‑236 (MQ=255)
                                                                                                                                                                           ctggctggCGCTAAGCCTGGCGGTGGCCTGCGTGCTGGCGCTGGGCAATATCAGCGATCGCA‑TGGAGAAGGGCTTAAGCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGTTCACGCGAAGTGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAGCGCGGCCTGAAAGTCGGCAAGCAGCTGACTTTCG‑CCACAATGACCTTTGCAGGCGACACACCGCAGCTGGCg                                                                   <  1:684134/238‑1 (MQ=255)
                                                                                                                                                                            tggcgggcGCTAAGGCTTGGGGGGGCCTGCGGGCCGGCGCTGGGCAATATCAGCGATCGCA‑TGGAGAAGGGCTTAAGCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGTTCACGCGAAGTGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAGCGCGGCCTGAAAGTCGGCAAGCAGCTGACTTTCG‑CCACAATGACCTTTGCAGGCGACACACCGCAGCTGGCGaa                                                                 <  2:214455/239‑1 (MQ=255)
                                                                                                                                                                                           cTGGCGTTGGCCTGCGTGCTGGCGCTGGGCAATGTCATCGATCGCA‑TGGAGAAGGGCTTAAGCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGTTCACGCGAAGTGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAGCGCGGCCTGAAAGTCGTCAAGCAGCTGACTTTCG‑CCACAATGACCTTTGCAGGCGACACACCGCAGCTGGCGAACGTCAAAGCGGTGga                                                  <  1:119671/239‑1 (MQ=255)
                                                                                                                                                                                            tGGCGGTGGCCTGCGTGCTGGCGCTGGGCAATATCAGCGATCGCA‑TGGAGAAGGGCTTAAGCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGTTCACGCGAAGTGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAGCGCGGCCTGAAAGTCGGCAAGCAGCTGACTTTCG‑CCACAATGACCTTTGCAGGCGACACACCGCAGCTGGCGAACGTCAAAGCGGTGgat                                                 <  1:145554/239‑1 (MQ=255)
                                                                                                                                                                                            gggggggggCCGGCGGGCGGGCGCGGGGAAATAACCGCGCTCGCA‑GGGGGAGGGGCGTAAGCCAGAAAAGCCGTGAGTTTATGGGGGGCGTTCGGGCGTTGCGCAGTTCACGCGAAGTGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAGCGCGGCCTGAAAGTCGGCAAGCAGCTGACTTTCG‑CCACAATGACCTTTGCAGGCGACACACCGCAGCTGGCGAACGTCAAAGCGGTGgat                                                 <  2:991178/232‑1 (MQ=255)
                                                                                                                                                                                                ggTGGCCTGCGTGCTGGCGCTGGGCAATTTCAGCGAGCGCA‑GGGAGAAGGGCTTAAGCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGTTCACGCGAAGTGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAGCGCGGCCTGAAAGTCGGCAAGCAGCTGACTTTCG‑CCACAATGACCTTTGCAGGCGACACACCGCAGCTGGCGAACGTCAAAGCGGTGgatgat                                              <  2:76837/238‑1 (MQ=255)
                                                                                                                                                                                                  tGGCCTGGGTGGTGGCGCTGGGCAATATTAGCGATCGCA‑TGGAGAAGGGCTTAAGCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGTTCACGCGAAGTGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAGCGCGGCCTGAAAGTCGGCAAGCAGCTGACTTTCG‑CCACAATGACCTTTGCAGGCGACACACCGCAGCTGGCGAACGTCAAAGCGGTGGATGatat                                            <  1:335079/238‑1 (MQ=255)
                                                                                                                                                                                                  tGGCCTGCGGGGTGGCGCTGGGCAATTTCAGCAGTCGAA‑TGGAGAAGGGCTTAAGCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGTTCACGCGAAGTGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAGCGCGGCCTGAAAGTCGGCAAGCAGCTGACTTTCG‑CCACAATGACCTTTGCAGGCGACACACCGCAGCTGGCGAACGTCAAAGCGGTGGATGatat                                            <  1:692058/238‑1 (MQ=255)
                                                                                                                                                                                                      cTGCGTGCTGGCGCTGGGCAATTTCAGCGATCGGA‑TGGAGAAGGGCTTAAGCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGTTCACGCGAAGTGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAGCGCGGCCTGAAAGTCGGCAAGCAGCTGACTTTCG‑CCACAATGACCTTTGCAGGCGACACACCGCAGCTGGCGAACGTCAAAGCGGTGGATGATATCTAc                                        <  2:453973/238‑1 (MQ=255)
                                                                                                                                                                                                       tGCGTGCTGGCGATGGGCAATATCAGCGATCGCA‑TGGAGAAGGGCTTAAGCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGTTCACGCGAAGTGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAGCGCGGCCTGAAAGTCGGCAAGCAGCTGACTTTCG‑CCACAATGACCTTTGCAGGCGACACACCGCAGCTGGCGAACGTCAAAGCGGTGGATGATATCTAcc                                       <  1:178464/238‑1 (MQ=255)
                                                                                                                                                                                                       ggCGTGCTGGCGCTGGGCAATAGCAGCGATCGCA‑TGGAGAAGGGCTTAAGCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGTTCACGCGAAGTGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAGCGCGGCCTGAAAGTCGGCAAGCAGCTGACTTTCG‑CCACAATGACCTTTGCAGGCGACACACCGCAGCTGGCGAACGTCAAAGCGGTGGATGATATCTAcc                                       <  1:302949/237‑1 (MQ=255)
                                                                                                                                                                                                        gggTGCTGGGGCTGGGCAATATCAGCGGTCGCA‑TGGAGAAGGGCTTAAGCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGTTCACGCGAAGTGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAGCGCGGCCTGAAAGTCGGCAAGCAGCTGACTTTCG‑CCACAATGACCTTTGCAGGCGACACACCGCAGCTGGCGAACGTCAAAGCGGTGGATGATATCTACCCg                                     <  2:788139/237‑1 (MQ=255)
                                                                                                                                                                                                                  gcTGGGCAATGTCAGCGATCGCA‑GGGAGAAGGGCTTAAGCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGTTCACGCGAAGTGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAGCGCGGCCTGAAAGTCGGCAAGCAGCTGACTTTCG‑CCACAATGACCTTTGCAGGCGACACACCGCAGCTGGCGAACGTCAAAGCGGTGGATGATATCTACCCGATGTATGGc                            <  1:747847/238‑1 (MQ=255)
                                                                                                                                                                                                                          ataaccgCGGTCGCA‑TGGAGAAGGGCTTAAGCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGTTCACGCGAAGTGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAGCGCGGCCTGAAAGTCGGCAAGCAGCTGACTTTCG‑CCACAATGACCTTTGCAGGCGACACACCGCAGCTGGCGAACGTCAAAGCGGTGGATGATATCTACCCGATGTATGGCGATCTGCa                    <  1:976841/232‑1 (MQ=255)
                                                                                                                                                                                                                           tatCAGCGATCGCA‑TGGAGAAGGGCTTAAGCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGTTCACGCGAAGGGCCGCAAGCGTGGCTGGAGGAAGCGCACAAGCGCGGCCTGCAAGCCGGCCAACAGCGGAGTTTCG‑CCCCAATGACCCTTTCAGGCGGCACACCGCCGCTGGGGACCGCCAAAACGGGGGGGGGTATCTTCCCCATGGATGGGGAtgcgacaa                  >  2:128659/1‑232 (MQ=255)
                                                                                                                                                                                                                               cggtgtCGCG‑TGGGGGAGGGCTTAAGCCTGCAAAGACTTGGGTTTATGGCTGGCGGTCGGGCGTTGCGCATTTCACGCGAAGTGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAGCGCGGCCTGAAAGTCGGCAAGCAGCTGACTTTCG‑CCACAATGACCTTTGCAGGCGACACACCGCAGCTGGCGAACGTCAAAGCGGTGGATGATATCTACCCGATGTATGGCGATCTGCAAACTaa              <  2:411039/234‑1 (MQ=255)
                                                                                                                                                                                                                                     cGCA‑TGGAGAAGGGCTTAAGCCAGGAAAGCCGGGAGTTTATGGGGGGCGATCGGGCGTTGCGCAGTTCACGCGAAGTGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAGCGCGGCCTGAAAGTCGGCAAGCAGCTGACTTTCG‑CCACAATGACCTTTGCAGGCGACACACCGCAGCTGGCGAACGTCAAAGCGGTGGATGATATCTACCCGATGTATGGCGATCTGCAAACTAATccccc        <  2:254531/239‑1 (MQ=255)
                                                                                                                                                                                                                                            gagaAGGGCTTAAGCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGTTCCCGCGAAGTGCCGCAAACGGGGCTGGGGGAAACGCCAAAACGCGGCCCGAAAGTCGGCAAGCAGCTTGCTTTCG‑CCCCCAAGACCTTTGCAGGCGACCAACCGCCGCCGGCGACCGCCACAACGGGGGAGGTTAACTACCCGATGAAGGGCGCGCTGccaccaacacacccccgccc  >  2:993010/1‑220 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                              
TGTTTTCCCTCAACCGTGAACATGGCACCACGTTGATTATGGTGACCCACGACCTGCAACTGGCGGCACGCTGCGACCGCTGCTTACGGCTGGTGAACGGGCAGTTGCAGGAGGAAGCATGATTGCACGTTGGTTCTGGCGCGAATGGCGTTCGCCGTCGCTATTAATTGTCTGGCTGGCGCTAAGCCTGGCGGTGGCCTGCGTGCTGGCGCTGGGCAATATCAGCGATCGCA‑TGGAGAAGGGCTTAAGCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGTTCACGCGAAGTGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAGCGCGGCCTGAAAGTCGGCAAGCAGCTGACTTTCG‑CCACAATGACCTTTGCAGGCGACACACCGCAGCTGGCGAACGTCAAAGCGGTGGATGATATCTACCCGATGTATGGCGATCTGCAAACTAATCCCCCTGGCCT  >  NC_000913/520298‑520771

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑