Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,329,709 C→G 5.5% G30A (GGT→GCT) ‡ yfaQ ← tandem DUF2300 domain protein, putative host defense protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,329,7090CG 5.5% 322.0 / 8.2 91G30A (GGT→GCT) ‡yfaQtandem DUF2300 domain protein, putative host defense protein
Reads supporting (aligned to +/- strand):  ref base C (43/43);  new base G (2/3);  total (45/46)
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

GGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTTGCTGCGCAGGTG  >  NC_000913/2329562‑2329856
                                                                                                                                                   |                                                                                                                                                   
ggCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCAcc                                                                                                                                                    >  1:435281/1‑149 (MQ=255)
      tCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCgcag                                                                                                                                              >  1:529485/1‑149 (MQ=255)
          gTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCGGTGCGCAGCAcg                                                                                                                                          <  2:201280/149‑2 (MQ=255)
          gTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCAcc                                                                                                                                          >  2:527995/1‑149 (MQ=255)
           ttGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCAcc                                                                                                                                          <  4:62892/148‑1 (MQ=255)
             gATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGt                                                                                                                                       <  1:957377/149‑1 (MQ=255)
              aTGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGtt                                                                                                                                      <  2:435281/149‑1 (MQ=255)
              aTGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGtt                                                                                                                                      <  1:322391/149‑1 (MQ=255)
              aTGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGt                                                                                                                                       <  4:16655/148‑1 (MQ=255)
                 ggTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGt                                                                                                                                    >  4:224965/1‑148 (MQ=255)
                 ggTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGt                                                                                                                                    >  4:275213/1‑148 (MQ=255)
                  gTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTaa                                                                                                                                  >  1:27525/1‑149 (MQ=255)
                  gTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTaa                                                                                                                                  <  1:527995/149‑1 (MQ=255)
                     tcctccAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACgg                                                                                                                               >  1:434983/1‑149 (MQ=255)
                      cctccAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACgg                                                                                                                               >  3:38260/1‑148 (MQ=255)
                       ctccAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCg                                                                                                                             >  1:713228/1‑149 (MQ=255)
                        tccAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCGGTGCgcagca                                                                                                                                            >  4:412850/1‑133 (MQ=255)
                        tccAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCGGTGCgcagca                                                                                                                                            <  3:412850/133‑1 (MQ=255)
                        tccAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCg                                                                                                                             >  4:351777/1‑148 (MQ=255)
                        tccAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGt                                                                                                                            <  1:337171/149‑1 (MQ=255)
                        tccAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTAAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGt                                                                                                                            <  1:402279/149‑1 (MQ=255)
                         ccagccagGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGt                                                                                                                            >  4:354582/1‑148 (MQ=255)
                         ccagccagGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGt                                                                                                                            >  4:597032/1‑148 (MQ=255)
                         ccagccagGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCCAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGtc                                                                                                                           <  2:699728/149‑1 (MQ=255)
                              cagGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTctcct                                                                                                                       >  3:541909/1‑148 (MQ=255)
                               agGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCtcctcc                                                                                                                     <  2:27525/149‑1 (MQ=255)
                                   gTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCaa                                                                                                                 <  1:1047153/149‑1 (MQ=255)
                                    tAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCaa                                                                                                                 <  3:535009/148‑1 (MQ=255)
                                     aaaCATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGTTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAg                                                                                                                <  4:598512/148‑1 (MQ=255)
                                       aCATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAgtgt                                                                                                             <  2:782548/149‑1 (MQ=255)
                                        cATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGtt                                                                                                            <  2:684977/149‑1 (MQ=255)
                                          tACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGGGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTgg                                                                                                          >  2:869923/1‑149 (MQ=255)
                                          tACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGtt                                                                                                            >  1:776397/1‑147 (MQ=255)
                                          tACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGtt                                                                                                            <  2:776397/147‑1 (MQ=255)
                                          tACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTg                                                                                                           <  3:354582/148‑1 (MQ=255)
                                              gTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTctcctc                                                                                                                      >  3:58007/1‑133 (MQ=255)
                                              gTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTctcctc                                                                                                                      <  4:58007/133‑1 (MQ=255)
                                              gTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCag                                                                                                       <  3:577420/148‑1 (MQ=255)
                                              gTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCag                                                                                                       <  4:429108/148‑1 (MQ=255)
                                               tttCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGTTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTATCGGCGTCTCCTCCGCAAGTGTTGGCaga                                                                                                      <  3:675270/148‑1 (MQ=255)
                                               tttCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCaga                                                                                                      <  3:275213/148‑1 (MQ=255)
                                                ttCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCagag                                                                                                     <  4:541909/148‑1 (MQ=255)
                                                   cacaAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAAc                                                                                                 >  1:374216/1‑149 (MQ=255)
                                                       aaCTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCaa                                                                                              >  3:625839/1‑148 (MQ=255)
                                                       aaCTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAAc                                                                                             <  2:713228/149‑1 (MQ=255)
                                                             cTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCaa                                                                                        <  4:789221/148‑1 (MQ=255)
                                                              tAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATa                                                                                      >  2:606677/1‑149 (MQ=255)
                                                                aCGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAg                                                                                     >  4:925049/1‑148 (MQ=255)
                                                                  ggCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGcc                                                                                   <  3:1003990/148‑1 (MQ=255)
                                                                    cGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAg                                                                                 <  4:38260/148‑1 (MQ=255)
                                                                    cGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGa                                                                                <  2:434983/149‑1 (MQ=255)
                                                                     gtggtgAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGa                                                                                <  3:597032/148‑1 (MQ=255)
                                                                       ggtgAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACg                                                                              >  3:926501/1‑148 (MQ=255)
                                                                       ggtgAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGa                                                                             >  1:819512/1‑149 (MQ=255)
                                                                         tgAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCGGTGCGCAGCAcg                                                                                                                                          <  3:844539/86‑2 (MQ=255)
                                                                         tgAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCGGTGCGCAGCAcg                                                                                                                                          >  4:844539/1‑85 (MQ=255)
                                                                                  cAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCa                                                                   >  3:703254/1‑148 (MQ=255)
                                                                                       aaTGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAAc                                                                                                                                 >  4:641013/1‑81 (MQ=255)
                                                                                       aaTGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAAc                                                                                                                                 <  3:641013/81‑1 (MQ=255)
                                                                                       aaTGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGcc                                                                                   >  1:1010826/1‑127 (MQ=255)
                                                                                       aaTGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGcc                                                                                   <  2:1010826/127‑1 (MQ=255)
                                                                                           cTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTa                                                          <  3:224965/148‑1 (MQ=255)
                                                                                           cTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGGTGATGGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGAGTCTCCTCTGCAAGTGTTGGCAGACTAACCAACGCCAATTGCCATACGCTGCTGCGCCAGTCCACGTAt                                                         >  2:443671/1‑149 (MQ=255)
                                                                                            tGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTaa                                                                                                                                  <  1:707848/75‑1 (MQ=255)
                                                                                            tGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTaa                                                                                                                                  >  2:707848/1‑75 (MQ=255)
                                                                                            tGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCa                                                                                                        >  1:229030/1‑101 (MQ=255)
                                                                                            tGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCa                                                                                                        <  2:229030/101‑1 (MQ=255)
                                                                                            tGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAAc                                                                                             <  2:979549/112‑1 (MQ=255)
                                                                                            tGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAAc                                                                                             >  1:979549/1‑112 (MQ=255)
                                                                                             gagaCTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTAtt                                                        >  4:595091/1‑148 (MQ=255)
                                                                                              agaCTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTAttt                                                       >  3:531372/1‑148 (MQ=255)
                                                                                                aCTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCAc                                                    >  2:757643/1‑149 (MQ=255)
                                                                                                   ttGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCt                                                  >  3:776623/1‑148 (MQ=255)
                                                                                                           gCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCgg                                         >  2:488415/1‑149 (MQ=255)
                                                                                                           gCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCgg                                         <  2:484246/149‑1 (MQ=255)
                                                                                                            ccTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCgg                                         >  4:838603/1‑148 (MQ=255)
                                                                                                              tgatgaCGACCGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTc                                       >  3:617449/1‑148 (MQ=255)
                                                                                                                  gacgacAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCAtt                                   >  4:133503/1‑148 (MQ=255)
                                                                                                                      acAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACt                               >  4:366931/1‑148 (MQ=255)
                                                                                                                       cAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTcc                             <  2:529485/149‑1 (MQ=255)
                                                                                                                       cAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTc                              >  3:84076/1‑148 (MQ=255)
                                                                                                                         gCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTccc                            <  3:838603/148‑1 (MQ=255)
                                                                                                                             aTACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGcc                        <  4:172193/148‑1 (MQ=255)
                                                                                                                                 aGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCaa                    <  3:351777/148‑1 (MQ=255)
                                                                                                                                   ttGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAaca                  >  3:644464/1‑148 (MQ=255)
                                                                                                                                   ttGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAaca                  >  4:357979/1‑148 (MQ=255)
                                                                                                                                   ttGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCACCCCGGTCCATTCACTCCCTGCCGCAaca                  <  4:1041812/148‑1 (MQ=255)
                                                                                                                                      aTCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACttt              >  1:494777/1‑149 (MQ=255)
                                                                                                                                      aTCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACttt              >  1:861051/1‑149 (MQ=255)
                                                                                                                                              gcgcACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTTGCTGCGCa      <  1:869923/149‑1 (MQ=255)
                                                                                                                                                 cACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTTGCTGCGCAGGt   <  1:606677/149‑1 (MQ=255)
                                                                                                                                                  aCCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTTGCTGCGCAGGTg  >  2:701744/1‑149 (MQ=255)
                                                                                                                                                   ccGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTc                                  >  3:22695/1‑116 (MQ=255)
                                                                                                                                                   ccGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTc                                  <  4:22695/116‑1 (MQ=255)
                                                                                                                                                   ccGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTTGCTGCGCAGGTg  >  4:714678/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTTGCTGCGCAGGTG  >  NC_000913/2329562‑2329856

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