Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,969,480 T→C 6.1% I255V (ATC→GTC)  fliC ← flagellin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,969,4800TC 6.1% 181.4 / 2.3 65I255V (ATC→GTC) fliCflagellin
Reads supporting (aligned to +/- strand):  ref base T (30/31);  new base C (2/2);  total (32/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.79e-01

CACCGGCCTCGGCCAGACCCAAGGCAGTGCTCAAGACGTCATCAGTGCTGGCATCGATGGTGAATGCCTGGGTGCCGTCGAGCTTCAGTGCACCAGTGGCGCTGTCAATGCTAGCGGTAACACCGGTATTGCTAGTTTCCAGGTTGATATTGTCAACCACATCCTGCAGGGTGTCGCCTGCCGCAGCAGCGATGGTGGTGGTGGTGGTACCAACGTTCAGCACCAAGTTACCGGCAGCAGCACCGATGGTCGGCTTATCAAATACGCCCGCGCCTGCAGTACCCACAACCGAAGC  >  NC_002947/4969333‑4969627
                                                                                                                                                   |                                                                                                                                                   
caccGGCCTCGGCCAGACCCAAGGCAGTGCTCAAGACGTCATCAGTGCTGGCATCGATGGTGAATGCCTGGGTGCCGTCGAGCTTCAGTGCACCAGTGGCGCTGTCAATGCTAGCGGTAACACCGGTATTGCTAGTTTCCAGGTTGata                                                                                                                                                    >  1:64012/1‑149 (MQ=255)
 accGGCCTCGGCCAGACCCAAGGCAGTGCTCAAGACGTCATCAGTGCTGGCATCGATGGTGAATGCCTGGGTGCCGTCGAGCTTCAGTGCCCCAGTGGCGCTGTCAATGCTAGCGGTAACACCGGTATTGCTAGTTTCCAGGTTGatat                                                                                                                                                   >  3:304900/1‑149 (MQ=255)
 accGGCCTCGGCCAGACCCAAGGCAGTGCTCAAGACGTCATCAGTGCTGGCATCGATGGTGAATGCCTGGGTGCCGTCGAGCTTCAGTGCACCAGTGGCGCTGTCAATGCTAGCGGTAACACCGGTATTGCTAGTTTCCAGGTTGatat                                                                                                                                                   >  1:334935/1‑149 (MQ=255)
       cTCGGCCAGACCCAAGGCAGTGCTCAAGACGTCATCAGTGCTGGCATCGATGGTGAATGCCTGGGTGCCGTCGAGCTTCAGTGCACCAGTGGCGCTGTCAATGCTAGCGGTAACACCGGTATTGCTAGTTTCCAGGTTGATATTGTCaa                                                                                                                                             >  5:363578/1‑149 (MQ=255)
        tCGGCCAGACCCAAGGCAGTGCTCAAGACGTCATCAGTGCTGGCATCGATGGTGAATGCCTGGGTGCCGTCGAGCTTCAGTGCACCAGTGGCGCTGTCAATGCTAGCGGTAACACCGGTATTGCTAGTTTCCAGGTTGATATTGTCAAc                                                                                                                                            >  2:94872/1‑149 (MQ=255)
        tCGGCCAGACCCAAGGCAGTGCTCAAGACGTCATCAGTGCTGGCATCGATGGTGAATGCCTGGGTGCCGTCGAGCTTCAGTGCACCAGTGGCGCTGTCAATGCTAGCGGTAACACCGGTATTGCTAGTTTCCAGGTTGACAATATCAAc                                                                                                                                            >  3:257451/1‑149 (MQ=255)
          ggCCAGACCCAAGGCAGTGCTCAAGACGTCATCAGTGCTGGCATCGATGGTGAATGCCTGGGTGCCGTCGAGCTTCAGTGCACCAGTGGCGCTGTCAATGCTAGCGGTAACACCGGTATTGCTAGTTTCCAGGTTGACAATATCAAcct                                                                                                                                          <  6:198119/149‑2 (MQ=255)
           gCCAGACCCAAGGCAGTGCTCAAGACGTCATCAGTGCTGGCATCGATGGTGAATGCCTGGGTGCCGTCGAGCTTCAGTGCACCAGTGGCGCTGTCAATGCTAGCGGTAACACCGGTATTGCTAGTTTCCAGGTTGATATTGTCAACca                                                                                                                                          >  7:100936/1‑148 (MQ=255)
              ataCCCAAGGCAGTGCTCAAGACGTCATCAGTGCTGGCATCGATTGTGTATGCCTGGGTGCCGTCGAGCTTCAGTGCACCAGTGGCGCTGTCAATGCTAGCGGTAACACCGGTATTGCTAGTTTCCAGGTTGATATTGTGAACCACATc                                                                                                                                      <  5:98361/147‑1 (MQ=255)
               gACCCAAGGCAGTGCTCAAGACGTCATCAGTGCTGGCATCGATGGTGAATGCCTGGGTGCCGTCGAGCTTCAGTGCACCAGTGGCGCTGTCAATGCTAGCGGTAACACCGGTATTGCTAGTTTCCAGGTTGACAATATCAAcctggaaa                                                                                                                                     <  4:257451/149‑7 (MQ=255)
                 cccAAGGCAGTGCTCAAGACGTCATCAGTGCTGGCATCGATGGTGAATGCCTGGGTGCCGTCGAGCTTCAGTGCACCAGTGGCGCTGTCAATGCTAGCGGTAACACCGGTATTGCTAGTTTCCAGGTTGATATTGTCAACCACATCCTg                                                                                                                                   >  1:134219/1‑149 (MQ=255)
                     aGGCAGTGCTCAAGACGTCATCAGTGCTGGCATCGATGGTGAATGCCTGGGTGCCGTCGAGCTTCAGTGCACCAGTGGCGCTGTCAATGCTAGCGGTAACACCGGTATTGCTAGTTTCCAGGTTGATATTGTCAACCACATCCTGCAgg                                                                                                                               <  1:410604/149‑1 (MQ=255)
                     aGGCAGTGCTCAAGACGTCATCAGTGCTGGCATCGATGGTGAATGCCTGGGTGCCGTCGAGCTTCAGTGCACCAGTGGCGCTGTCAATGCTAGCGGTAACACCGGTATTGCTAGTTTCCAGGTTGATATTGTCAACCACATCCTGCAgg                                                                                                                               <  1:417069/149‑1 (MQ=255)
                     aGGCAGTGCTCAAGACGTCATCAGTGCTGGCATCGATGGTGAATGCCTGGGTGCCGTCGAGCTTCAGTGCACCAGTGGCGCTGTCAATGCTAGCGGTAACACCGGTATTGCTAGTTTCCAGGTTGATATTGTCAACCACATCCTGCAgg                                                                                                                               <  7:204701/149‑1 (MQ=255)
                      ggCAGTGCTCAAGACGTCATCAGTGCTGGCATCGATGGTGAATGCCTGGGTGCCGTCGAGCTTCAGTGCACCAGTGGCGCTGTCAATGCTAGCGGTAACACCGGTATTGCTAGTTTCCAGGTTCATATTGTCAACCACATCCTGCAggg                                                                                                                              >  4:231010/1‑149 (MQ=255)
                       gCAGTGCTCAAGACGTCATCAGTGCTGGCATCGATGGTGAATGCCTGGGTGCCGTCGAGCTTCAGTGCACCAGTGGCGCTGTCAATGCTAGCGGTAACACCGGTATTGCTAGTTTCCAGGTTGATATTGTCAACCACATCCTGCAGGgt                                                                                                                             >  2:395982/1‑149 (MQ=255)
                       gCAGTGCTCAAGACGTCATCAGTGCTGGCATCGATGGTGAATGCCTGGGTGCCGTCGAGCTTCAGTGCACCAGTGGCGCTGTCAATGCTAGCGGTAACACCGGTATTGCTAGTTTCCAGGATGATATTGTCAAc                                                                                                                                            >  8:183004/1‑134 (MQ=255)
                             cTCAAGACGTCATCAGTGCTGGCATCGATGGTGAATGCCTGGGTGCCGTCGAGCTTCAGTGCACCAGTGGCGCTGTCAATGCTAGCGGTAACACCGGTATTGCTAGTTTCCAGGTTGATATTGTCAACCACATCCTGCAGGGTGTCGcc                                                                                                                       <  6:225085/149‑1 (MQ=255)
                               cAAGACGTCATCAGTGCTGGCATCGATGGTGAATGCCTGGGTGCCGTCGAGCTTCAGTGCACCAGTGGCGCTGTCAATGCTAGCGGTAACACCGGTATTGCTAGTTTCCAGGTTGATATTGTCAACCACATCCTGCAGGGTGTCGCCTg                                                                                                                     >  4:263707/1‑149 (MQ=255)
                                aaGACGTCATCAGTGCTGGCATCGATGGTGAATGCCTGGGTGCCGTCGAGCTTCAGTGCACCAGTGGCGCTGTCAATGATAGCGGTAACACCGGTATTGCTAGTTTCCAGGTTGATATTGTCAACCACATCCTGCAGGGTGTCGCCTGc                                                                                                                    >  6:131339/1‑149 (MQ=255)
                                 aGACGTCATCAGTGCTGGCATCGATGGTGAATGCCTGGGTGCCGTCGAGCTTCAGTGCACCAGTGGCGCTGTCAATGCTAGCGGTAACACCGGTATTGCTAGTTTCCAGGTTGATATTGTCAACCACATCCTGCAGGGTGTCGCCTGcc                                                                                                                   <  7:83399/149‑1 (MQ=255)
                                 aGACGTCATCAGTGCTGGCATCGATGGTGAATGCCTGGGTGCCGTCGAGCTTCAGTGCACCAGTGGCGCTGTCAATGCTAGCGGTAACACCGGTATTGCTAGTTTCCAGGTTGATATTGTCAACCACATCCTGCAGGGTGTCGCCTGcc                                                                                                                   <  5:260284/149‑1 (MQ=255)
                                     atcatcaGTTCAGGCATCGATGGTGAATGCCTGGGTGACGTCGAGCTTAAGTGCAACAGTGGCGATGTCAAAGCTATCGGTAACAGCGGTATTGCTAGTTTCCATGTTGATATTTTCAACCCCAGCCTGCAGGGTGAGGCCTGCGgcac                                                                                                               >  2:78913/2‑148 (MQ=255)
                                        atcaGTGCTGGCATCGATGGTGAATGCCTGGGTGCCGTCGAGCTTCAGTGCACCAGTGGCGCTGTCAATGCTAGCGGTAACACCGGTATTGCTAGTTTCCAGGTTGATATTGTCAACCACATCCTGCAGGGTGTCGCCTGCCGcagcag                                                                                                            <  4:117208/149‑1 (MQ=255)
                                          caGTGCTGGCATCGATGGTGAATGCCTGGGTGCCGTCGAGCTTCAGTGCACCAGTGGCGCTGTCAATGCTAGCGGTAACACCGGTATTGCTAGTTTCCAGGTTGATATTGTCAACCACATCCTGCAGGGTGTTGCCTGCCGCagcagc                                                                                                           >  3:347176/1‑148 (MQ=255)
                                            gTGCTGGCATCGATGGTGAATGCCTGGGTGCCGTCGAGCTTCAGTGCACCAGTGGCGCTGTCAATGCTAGCGGTAACACCGGTATTGCTAGTTTCCAGGTTGATATTGTCAACCACATCCTGCAGGGTGTCGCCTGCCGCAGCAGCGAt                                                                                                        <  5:207160/149‑1 (MQ=255)
                                                   cATCGATGGTGAATGCCTGGGTGCCGTCGAGCTTCAGTGCACCAGTGGCGCTGTCAATGCTAGCGGTAACCCCGGTATTGCTAGTTTCCAGGTTGATATTGTCAACCCAATCCTGCAGGGTGTCGCCTGCCGCAGCAGCGAtggtggtg                                                                                                 >  5:202036/1‑149 (MQ=255)
                                                    aTCGATGGTGAATGCCTGGGTGCCGTCGAGCTTCAGTGCACCAGTGGCGCTGTCAATGCTAGCGGTAACACCGGTATTGCTAGTTTCCAGGTTGATATTGTCAACCACATCCTGCAGGGTGTCGCCTGCCGCAGCAGCGAtggtggtgg                                                                                                >  6:273719/1‑149 (MQ=255)
                                                    aTCGATGGTGAATGCCTGGGTGCCGTCGAGCTTCAGTGCACCAGTGGCGCTGTCAATGCTAGCGGTAACACCGGTATTGCTAGTTTCCAGGTTGATATTGTCAACCACATCCTGCAGGGTGTCGCCTGCCGCAGCAGCGAtggtggtgg                                                                                                <  1:94872/149‑1 (MQ=255)
                                                      agATGGTGAATGCCTGGGTGCCGTAGAGCTTCAGTGCAAAAGTGGCGCTGTCAATGCTAGCGGTAAAACCGGTATTGCTAGTTTCCAGGTTGATATTGTCAACCACATACTGAATGGTGTCGACTGCCGCAGCATCGAtggtggtggtg                                                                                              >  1:304412/2‑149 (MQ=255)
                                                       gATGGTGAATGCCTGGGTGCCGTCGAGCTTCAGTGCACCAGTGGCGCTGTCAATGCTAGCGGTAACACCGGTATTGCTAGTTTCCAGGTTGATATTGTCAACCACATCCTGCAGGGTGTCGCCTGCCGCAGCAGCGAtggtggtggtgg                                                                                             <  6:435494/149‑1 (MQ=255)
                                                        aTGGTGAATGCCTGGGTGCCGTCGAGCTTCAGTGCACCAGTGGCGCTGTCAATGCTAGCGGTAACACCGGTATTGCTAGTTTCCAGGTTGATATTGTCAACCACATCCTGCAGGGTGTCGCCTGCCGCAGCAGCGAtggtggtggtggt                                                                                            <  6:363578/149‑1 (MQ=255)
                                                         tGGTGAATGCCTGGGTGCCGTCGAGCTTCAGTGCACCAGTGGCGCTGTCAATGCTAGCGGTAACACCGGTATTGCTAGTTTCCAGGTTGATATTGTCAACCACATCCTGCAGGGTGTCGCCTGCCGCAGCAGCGAtggtggtggtggtg                                                                                           >  1:115634/1‑149 (MQ=255)
                                                              aaTGCCTGGGTGCCGTCGAGCTTCAGTGCACCAGTGGCGCTGTCAATGCTAGCGGTAACACCGGTATTGCTAGTTTCCAGGTTGATATTGTCAACCACATCCTGCAGGGTGTCGCCTGCCGCAGCAGCGATGGTGGTGGTGGTGGTAc                                                                                       <  2:64012/148‑1 (MQ=255)
                                                                      ggTGCCGTCGAGCTTCAGTGCACCAGTGGCGCTGTCAATGCTAGCGGTAACACCGGTATTGCTAGTTTCCAGGTTGATATTGTCAACCACATCCTGCAGGGTGTCGCCTGCCGCAGCAGCGATGGTGGTGGTGGTGGTACCAACGTTCa                                                                              <  8:100936/149‑1 (MQ=255)
                                                                        tGCCGTCGAGCTTCAGTGCACCAGTGGCGCTGTCAATGCTAGCGGTAACACCGGTATTGCTAGTTTCCAGGTTGATATTGTCAACCACATCCTGCAGGGTGTCGCCTGCCGCAGCAGCGATGGTGGTGGTGGTGGTACCAACGTTCAGc                                                                            <  2:334935/149‑1 (MQ=255)
                                                                         gCCGTCGAGCTTCAGTGCACCAGTGGCGCTGTCAATGCTAGCGGTAACACCGGTATTGCTAGTTTCCAGGTTGATATTGTCAACCACATCCTGCAGGGTGTCGCCTGCCGCAGCAGCGATGGTGGTGGTGGTGGTACCAACGTTCAGCa                                                                           >  1:2040/1‑149 (MQ=255)
                                                                             tCGAGATACAGACCACCATTGGCGCTGTCAATGCTAGCGGTAACAACGCTATAGCTAGATTCCAGCAAGATATTGTCAACAAAAACATGAAGAGTGTCGACTGCCGCAGCAGAGATGCTTCTGGTGGTGTTACCAACGTTCAGCACCaa                                                                       >  2:408010/1‑149 (MQ=255)
                                                                                   ttCAGTGCACCAGTGGCGCTGTCAATGCTAGCGGTAACACCGGTCTTGCTAGTTTCCAGGTTGATATTGTCAACCACATCCTGCAGGGTGTCGCCTGCCGCAGCAGCGATGGTGGTGGTGGTGGTACCAAc                                                                                   <  7:228631/131‑1 (MQ=255)
                                                                                   ttCAGTGCACCAGTGGCGCTGTCAATGCTAGCGGTAACACCGGTATTGCTAGTTTCCAGGTTGATATTGTCAACCACATCCTGCAGGGTGTCGCCTGCCGCAGCAGCGATGGTGGTGGTGGTGGTACCAAc                                                                                   >  8:228631/1‑131 (MQ=255)
                                                                                   ttCAGTGCACCAGTGGCGCTGTCAATGCTAGCGGTAACACCGGTATTGCTAGTTTCCAGGTTGATATTGTCAACCACATCCTGCAGGGTGTCGCCTGCCGCAGCAGCGATGGTGGTGGTGGTGGTACCAACGtt                                                                                <  4:110119/134‑1 (MQ=255)
                                                                                   ttCAGTGCACCAGTGGCGCTGTCAATGCTAGCGGTAACACCGGTATTGCTAGTTTCCAGGTTGATATTGTCAACCACATCCTGCAGGGTGTCGCCTGCCGCAGCAGCGATGGTGGTGGTGGTGGTACCAACGtt                                                                                >  3:110119/1‑134 (MQ=255)
                                                                                     cAGTGCACCAGTGGCGCTGTCAATGCTAGCGGTAACACCGGTATTGCTAGTTTCCAGGTTGATATTGTCAACCACATCCTGCAGGGTGTCGCCTGCCGCAGCAGCGATGGTGGTGGTGGTGGTACCAACGTTCAGCACCAAGTTAcc                                                                 <  3:231010/147‑1 (MQ=255)
                                                                                       gTGCACCAGTGGCGCTGTCAATGCTAGCGGTAACACCGGTATTGCTAGTTTCCAGGTTGATATTGTCAACCACATCCTGCAGCGTGTCGCCTGCCGCAGCAGCGATGGTGGTGGTGGTGGTACCAACGTTCAGCACCAAGTTACCgg                                                               >  1:213432/1‑147 (MQ=255)
                                                                                          cACCAGTGGCGCTGTCAATGCTAGCGGTAACACCGGTATTGCTAGTTTCCAGGTTGATATTGTCAACCACATCCTGCAGGGTGTCGCCTGCCGCAGCAGCGATGGTGGTGGTGGTGGTACCAACGTTCAGCACCAAGTTACCGgcag                                                            >  4:395166/1‑147 (MQ=255)
                                                                                            ccAGTGGCGCTGTCAATGCTAGCGGTAACACCGGTATTGCTAGTTTCCAGGTTGATATTGTCAACCACATCCTGCAGGGTGTCGCCTGCCGCAGCAGCGATGGTGGTGGTGGTGGTACCAACGTTCAGCACCAAGTTACCGgcagcagc                                                        <  2:134219/149‑1 (MQ=255)
                                                                                              aGTGGCGCTGTCAATGCTAGCGGTAACTGCGGTTTTTCTAGTTTCCAGGTTGATATTGTCAACCACATCCTGCAGGGTGTCGCCTGCCGCAGCAGCGATGGTGGTGGAGGTGGTACCAACGTTCAGCACCAAGTTACCGGCAGCAGCAc                                                      <  1:408010/149‑1 (MQ=255)
                                                                                                        tCAATGCTAGCGGTAACACCGGTATTGCTAGTTTCCAGGTTGATATTGTCAACCACATCCTGCAGGGTGTCGCCTGCCGCAGCAGCGATGTTGGTGGTGGTGGTACCAACGTTCAGCACCAAGTTACCGGCAGCAGCACCGATGGTCgg                                            >  7:379284/1‑149 (MQ=255)
                                                                                                          aaTGCTAGCGGTAACACCGGTATTGCTAGTTTCCAGGTTGATATTGTCAACCACATCCTGCAGGGTGTTGCCTGCCGCAGCAGCGATGGTGGTGGTGGTGGTACCAACGTTCAGCACCAAGTTACCGGCAGCAGCACCGATGGTCGGc                                           <  4:347176/148‑1 (MQ=255)
                                                                                                          aaTGCTAGCGGTAACACCGGTATTGCTAGTTTCCAGGTTGATATTGTCAACCACATCCTGCAGGGTGTCGCCTGCCGCAGCAGCGATGGTGGTGGTGGTGGTACCAACGTTCAGCACCAAGTTACCGGCAGCAGCACCGATGGTCGGc                                           <  2:155912/148‑1 (MQ=255)
                                                                                                           aTGCTAGCGGTAACACCGGTATTGCTAGTTTCCAGGTTGATATTGTCAACCACATCCTGCAGGGTGTCGCCTGCCGCAGCAGCGATGGTGGTGGTGGTGGTACCAACGTTCAGCACCAAGTTACCGGCAGCAGCACCGATGGTCGGCtt                                         <  5:273719/149‑1 (MQ=255)
                                                                                                            tGCTAGCGGTAACACCGGTATTGCTAGTTTCCAGGTTGATATTGTCAACCACATCCTGCAGGGTGTCGCCTGCCGCAGCAGCGATGGTGGTGGTGGTGGTACCAACGTTCAGCAcc                                                                         <  6:149326/116‑1 (MQ=255)
                                                                                                            tGCTAGCGGTAACACCGGTATTGCTAGTTTCCAGGTTGATATTGTCAACCACATCCTGCAGGGTGTCGCCTGCCGCAGCAGCGATGGTGGTGGTGGTGGTACCAACGTTCAGCAcc                                                                         >  5:149326/1‑116 (MQ=255)
                                                                                                               tAGCGGTAACACCGGTATTGCTAGTTTCCAGGTTGATATTGTCAACCACATCCTTCAGGGTGTCGCCTGCCGCAGCAGCGATGGT‑‑T‑GTGGTGGTACCAACGTTCAGCACCAAGTTACCGGCAGCAGCACCGATGGTCGGCTTATCaaa                                   >  1:297970/1‑148 (MQ=255)
                                                                                                                 gCGGTAACACCGGTATTGCTAGTTTCCAGGTTGATATTGTCAACCACATCCTGCAGGGTGTCGCCTGCCGCAGCAGCGATGGTGGTGGTGGTGGTACCAACGTTCAGCACCAAGTTACCGGCAGCAGCACCGATGGTCGGCTTATCaaa                                   <  4:304900/149‑1 (MQ=255)
                                                                                                                     tAACACCGGTATTGCTAGTTTCCAGGTTGATATTGTCAACCACATCCTGCAGGGTGTCGCCTGCCGCAGCAGCGATGGTGGTGGTGGTGGTACCAACGTTCAGCACCAAGTTACCGGCAGCAGCACCGATGGTCGGCTTATCAAATACg                               <  6:202036/149‑1 (MQ=255)
                                                                                                                          ccGGTATTGCTAGTTTCCAGGTTGATATTGTCAACCACATCCTGCAGGGTGTCGCCTGCCGCAGCAGCGATGGTGGTGGTGGTGGTACCAACGTTCAGCACCAAGTTACCGGCAGCAGCACCGATGGTCGGCTTATCAAATACGCCcgc                          <  2:2040/149‑1 (MQ=255)
                                                                                                                                ttGCTAGTTTCCAGGTTGATATTGTCAACCACATCCTGCAGGGTGTCGCCTGCCGCAGCAGCGATGGTGGTGGTGGTGGTACCAACGTTCAGCACCAAGTTACCGGCAGCAGCACCGATGGTCGGCTTATCAAATACGCCCGCGCCTg                     >  1:127907/1‑148 (MQ=255)
                                                                                                                                   cTAGTTTCCATGTTGATATTGTCAACCACATACTGCAGGGTGTCTCCTGCCGCAGCAGCGATGGTGGTTGTTGTTGTACCAACGTTCAGCACCAAGTTACCGGCAGCAGCACCGATTGTCGGCTTATCAAATACGCCCGCGCCTGCat                  <  2:304412/148‑2 (MQ=255)
                                                                                                                                    tAGTTTCCAGGTTGATATTGTCAACCACATCCTGCAGGGTGTCGCCTGCCGCAGCAGCGATGGTGGTGGTGGTGGTACCAACGTTCAGCACCAAGTTACCGGCAGCAGCACCGATGGTCGGCTTATCAAATACGCCCGCGACTGCAGTa                >  7:194957/1‑149 (MQ=255)
                                                                                                                                      gTTTCCAGGTTGATATTGTCAACCACATCCTGCAGGGGGTCGCCTGCCGCAGCAGCGATGGTGGTGGTGGTGGTACCAACGTTCAGCACCAAGTTACCGGCAGCAGCACCGATGGTCGGCTTATCAAATACGCCCGCGCCTGCAGTAcc              <  6:231187/149‑1 (MQ=255)
                                                                                                                                       tttCCAGGTTGATATTGTCAACCACATCCTGCAGGGTGTCGCCTGCCGCAGCAGCGATGGTGGTGGTGGTGGTACCAACGTTCAGCACCAAGTTACCGGCAGCAGCACCGATGGTCGGCTTATCAAATACGCCCGCGCCTGCAGTAcc              >  8:61094/1‑148 (MQ=255)
                                                                                                                                        gatgtggTTGACAATATCAACCACATCCTGCAGGGTGTCGCCTGCCGCAGCAGCGATGGTGGTGGTGGTGGTACCAACGTTCAGCACCAAGTTACCGGCAGCAGCACCGATGGTCGGCTTATCAAATACGCCCGCGCCTGCAGTACCca            >  6:18507/6‑149 (MQ=255)
                                                                                                                                           cAGGTGGATCTTGTCAACCACTTCGTGGATGGTGTCGCGTGCCGCAGCAGCGATGGTGGTGTTGGTTGTACCAACGTTCAGCACCAAGTTACCGGAAGCAGCACCGATGGTCGGCTTTTCAAATAGGCACGCGCCTGCAGTAGCCACaa         <  8:194957/149‑1 (MQ=255)
                                                                                                                                            aGGTTGATATTGTCAACCACATCCTGCAGGGTGTCGCCTGCCGCAGCAGCGATGGTGGTGGTGGTGGTACCAACGTTCAGCACCAAGTTACCGGCAGCAGCACCGATGGTCGGCTTATCAAATACGCCCGCGCCTGCAGTACCCACAAc        >  7:379013/1‑149 (MQ=255)
                                                                                                                                            aGGTTGATATTGTCAACCACATCCTGCAGGGTGTCGCCTGCCGCAGCAGCGATGGTGGTGGTGGTGGTACCAACGTTCAGCACCAAGTTACCGGCAGCAGCACCGATGGTCGGCTTATCAAATACGCCCGCGCCTGCAGTACCCACAAc        >  3:345025/1‑149 (MQ=255)
                                                                                                                                              gTTGATATTGTCAACCACATCCTGCAGGGTGTCGCCTGCCGCAGCAGCGATGGTGGTGGTGGTGGTACCAACGTTCAGCACCAAGTTACCGGCAGCAGCACCGATGGTCGGCTTATCAAATACGCCCGCGCCTGCAGTACCCACAAcc       >  2:428333/1‑148 (MQ=255)
                                                                                                                                                 gATATTGTCAACCACATCCTGCAGGGTGTCGCCTGCCGCAGCAGCGATGGTGGTGGTGGTGGTACCAACGTTCAGCACCAAGTTACCGGCAGCAGCACCGATGGTCGGCTTATCAAATACGCCCGCGCCTGCAGTACCCACAACCGaa    <  8:379284/148‑1 (MQ=255)
                                                                                                                                                   tatTGTCAACCACATCCTGCAGGGTGTCGCCTGCCGCAGCAGCGATGGTGGTGGTGGTGGTACCAACGTTCAGCACCAAGTTACCGGCAGCAGCACCGATGGTCGGCTTATCAAATACGCCCGCGCCTGCAGTACCCACAACCGAAGc  <  2:213432/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CACCGGCCTCGGCCAGACCCAAGGCAGTGCTCAAGACGTCATCAGTGCTGGCATCGATGGTGAATGCCTGGGTGCCGTCGAGCTTCAGTGCACCAGTGGCGCTGTCAATGCTAGCGGTAACACCGGTATTGCTAGTTTCCAGGTTGATATTGTCAACCACATCCTGCAGGGTGTCGCCTGCCGCAGCAGCGATGGTGGTGGTGGTGGTACCAACGTTCAGCACCAAGTTACCGGCAGCAGCACCGATGGTCGGCTTATCAAATACGCCCGCGCCTGCAGTACCCACAACCGAAGC  >  NC_002947/4969333‑4969627

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑