Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,680,919 A→G 5.0% S122P (TCC→CCC)  PP_4986 ← hemolysin III family channel protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,680,9190AG 5.0% 238.2 / 2.1 79S122P (TCC→CCC) PP_4986hemolysin III family channel protein
Reads supporting (aligned to +/- strand):  ref base A (31/44);  new base G (2/2);  total (33/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 = 5.83e-01

AAATGATGCCGACGGTGTAGAACACGCCGCCGGCAGCCAGCCAGGCGAAGCCGGCCGTGCCCAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTAC‑‑CGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAGGGCCGCGCAGGCTGACCAGGCAAAACGGCGTGTAGCTGCCGGCAATCAGCAGATAGATCGACAGGTGGTCGA  >  NC_002947/5680772‑5681065
                                                                                                                                                     |                                                                                                                                                  
aaaTGATGCCGACGGTGTAGAACACGCCGCCGGCAGCCAGCCAGGCGAAGCCGGCCGTGCCCAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTAC‑‑CGCATAGATGATGATCGACAGGATCCGCGCCTCGGAg                                                                                                                                                   >  2:251280/1‑149 (MQ=255)
   tgatgCCGACGGTGTAGAACACGCCGCCGGCAGCCAGCCAGGCGAAGCCGGCCGTGCCCAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTAC‑‑CGCATAGATGATGATCGACAGGATCCGCGCCTCGGAgcgc                                                                                                                                                >  6:386532/1‑149 (MQ=255)
      tgCCGACGGTGTAGAACACGCCGCCGGCAGCCAGCCAGGCGAAGCCGGCCGTGCCCAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTAC‑‑CGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGc                                                                                                                                             <  1:420579/149‑1 (MQ=255)
        ccGACGGTGTAGAACACGCCGCCGGCAGCCAGCCAGGCGAAGCCGGCCGTGCCCAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTAC‑‑CGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCtt                                                                                                                                           <  5:326262/149‑1 (MQ=255)
        ccGACGGTGTAGAACACGCCGCCGGCAGCCAGCCAGGCGAAGCCGGCCGTGCCCAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTAC‑‑CGCACAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCt                                                                                                                                            <  1:323325/148‑1 (MQ=255)
         cGACGGTGTAGAACACGCCGCCGGCAGCCAGCCAGGCGAAGCCGGCCGTGCCCAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTAC‑‑CGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCtt                                                                                                                                           >  5:339795/1‑148 (MQ=255)
           aCGGTGTAGAACACGCCGCCGGCAGCCAGCCAGGCGAAGCCGGCCGTGCCCAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTAC‑‑CGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGAt                                                                                                                                        <  8:204261/149‑1 (MQ=255)
              gtgtAGAACACGCCGCCGGCAGCCAGCCAGGCGAAGCCGGCCGTGCCCAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTAC‑‑CGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGAtctc                                                                                                                                     <  5:255134/149‑1 (MQ=255)
               tgtAGAACACGCCGCCGGCAGCCAGCCAGGCGAAGCCGGCCGTGCCCAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTAC‑‑CGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATctct                                                                                                                                    >  8:319064/1‑149 (MQ=255)
                   gAACACGCCGCCGGCAGCCAGCCAGGCGAAGCCGGCCGTGCCCAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTAC‑‑CGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGc                                                                                                                                 >  3:396149/1‑148 (MQ=255)
                    aaCACGCCGCCGGCAGCCAGCCAGGCGAAGCCGGCCGTGCCCAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTAC‑‑CGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCaa                                                                                                                               <  4:394099/149‑1 (MQ=255)
                    aaCACGCCGCCGGCAGCCAGCCAGGCGAAGCCGGCCGTGCCCAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTAC‑‑CGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCaa                                                                                                                               >  5:153975/1‑149 (MQ=255)
                         gccgccGGCAGCCAGCCAGGCGAAGCCGGCCGTGCCCAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTAC‑‑CGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGc                                                                                                                          <  3:123962/149‑1 (MQ=255)
                           cgccgGCAGCCAGCCAGGCGAAGCCGGCCGTGCCCAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTAC‑‑CGCATAGATGATGATCGACAGGATCCGCGCTCCGAGGCGCgga                                                                                                                                             <  5:204540/128‑2 (MQ=255)
                           cgccgGCAGCCAGCCAGGCGAAGCCGGCCGTGCCCAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTAC‑‑CGCATAGATGATGATCGACAGGATCCGCGCTCCGAGGCGCgga                                                                                                                                             >  6:204540/1‑127 (MQ=255)
                                 cagccagcCAGGCGAAGCCGGCCGTGCCCAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTAC‑‑CGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATc                                                                                                                  <  2:357814/149‑1 (MQ=255)
                                 cagccagcCAGGCGAAGCCGGCCGTGCCCAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTAC‑‑CGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATc                                                                                                                  <  5:212220/149‑1 (MQ=255)
                                   gccagccaGGCGAAGCCGGCCGTGCCCAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTAC‑‑CGCATAGATGATGATCGACAGGATCCGCGCTCCGAGGCGCgga                                                                                                                                             >  7:108607/1‑119 (MQ=255)
                                   gccagccaGGCGAAGCCGGCCGTGCCCAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTAC‑‑CGCATAGATGATGATCGACAGGATCCGCGCTCCGAGGCGCgga                                                                                                                                             <  8:108607/120‑2 (MQ=255)
                                   gccagccaGGCGAAGCCGGCCGTGCCCAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTAC‑‑CGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGc                                                                                                                >  2:323085/1‑149 (MQ=255)
                                        ccagGCGAAGCCGGCCGTGCCCAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTAC‑‑CGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGcc                                                                                                           >  7:154597/1‑149 (MQ=255)
                                         cagGCGAAGCCGGCCGTGCCCAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTAC‑‑CGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGccc                                                                                                          <  8:70819/149‑1 (MQ=255)
                                          agGCGAAGCCGGCCGTGCCCAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTAC‑‑CGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGcccc                                                                                                         <  8:154597/149‑1 (MQ=255)
                                           ggCGAAGCCGGCCGTGCCCAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTAC‑‑CGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGcccc                                                                                                         >  4:373854/1‑148 (MQ=255)
                                            gCGAAGCCGGCCGTGCCCAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTAC‑‑CGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGcccc                                                                                                         >  5:51856/1‑147 (MQ=255)
                                            gCGAAGCCGGCCGTGCCCAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTAC‑‑CGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGccac                                                                                                         >  8:85835/1‑145 (MQ=255)
                                               aaGCCGGCCGTGCCCAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTAC‑‑CGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGAc                                                                                                    <  4:259894/149‑1 (MQ=255)
                                                aGCCGGCCGTGCCCAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTAC‑‑CGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGAcc                                                                                                   <  1:251280/149‑1 (MQ=255)
                                                  ccggccgTGCCCAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTAC‑‑CGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGAtct                                                                                                                                      >  7:287581/1‑112 (MQ=255)
                                                  ccggccgTGCCCAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTAC‑‑CGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGAtct                                                                                                                                      <  8:287581/112‑1 (MQ=255)
                                                     gccgTGCCCAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTAC‑‑CGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGcc                                                                                              <  7:319064/149‑1 (MQ=255)
                                                      ccgTGCCCAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTAC‑‑CGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCg                                                                                             <  8:233768/149‑1 (MQ=255)
                                                      ccgTGCCCAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTAC‑‑CGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCg                                                                                             <  1:397461/149‑1 (MQ=255)
                                                            ccAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTAC‑‑CGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAgg                                                                                       >  2:378867/1‑149 (MQ=255)
                                                            ccAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTAC‑‑CGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAgg                                                                                       <  2:22455/149‑1 (MQ=255)
                                                             cAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTAC‑‑CGCATAGATGATGATCGACAGGATCCGCGCCTCTGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGc                                                                                      <  7:152141/149‑1 (MQ=255)
                                                              acctaGTGCAGAAGCGGCTTGACTCCCACCTGCAGCAT‑CACGCCATTACTTCACTTAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGc                                                                                      <  1:359071/145‑1 (MQ=255)
                                                               gCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTAC‑‑CGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTc                                                                                    <  2:286879/149‑1 (MQ=255)
                                                                   gTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTAC‑‑CGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATc                                                                                >  6:50302/1‑149 (MQ=255)
                                                                    ttCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTAC‑‑CGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATcc                                                                               >  2:280199/1‑149 (MQ=255)
                                                                    ttCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTAC‑‑CGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATcc                                                                               >  4:107457/1‑149 (MQ=255)
                                                                          agcGGCTTGACTGCCACCAGCACGATCCAGCCCATTAC‑‑CGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAgg                                                                         <  8:204477/149‑1 (MQ=255)
                                                                               cTTGACTGCCACCAGCACGATCCAGCCCATTAC‑‑CGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAGGGCcgc                                                                    >  4:220647/1‑149 (MQ=255)
                                                                                ttGACTGCCACCAGCACGATCCAGCCCATTAC‑‑CGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAGGGCcgcg                                                                   >  6:184076/1‑149 (MQ=255)
                                                                                  gACTGCCACCAGCACGATCCAGCCCATTAC‑‑CGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAGGGCcgcg                                                                   <  5:240137/147‑1 (MQ=255)
                                                                                   aCTGCCACCAGCACGATCCAGCCCATTAC‑‑CGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGTCCACGCCGAACAGGCTCCATCCCCAAGGGCCGCGCa                                                                 <  1:289953/148‑1 (MQ=255)
                                                                                   aCTGCCACCAGCACGATCCAGCCCATTAC‑‑CGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAGGGCCGCGCAg                                                                <  8:357215/149‑1 (MQ=255)
                                                                                    cTGCCACCAGCACGATCCAGCCCATTAC‑‑CGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAGGGCCGCGCAgg                                                               <  5:50302/149‑1 (MQ=255)
                                                                                    cTGCCACCAGCACGATCCAGCCCATTAC‑‑CGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAGGGCCGCGCAgg                                                               <  5:294840/149‑1 (MQ=255)
                                                                                    cTGCCACCAGCACGATCCAGCCCATTAC‑‑CGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAGGGCCGCGCAgg                                                               <  1:174877/149‑1 (MQ=255)
                                                                                    cTGCCACCAGCACGATCCAGCCCATTAC‑‑CGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAGGGCCGCGCAgg                                                               >  8:311429/1‑149 (MQ=255)
                                                                                       ccaccaGCACGATCCAGCCCATTAC‑‑CGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACg                                                                                                <  2:165179/113‑1 (MQ=255)
                                                                                       ccaccaGCACGATCCAGCCCATTAC‑‑AGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCCGCCCCCAGACCACGCCGAACAGCCTCCATCCCCAAGGGCCGCGCAGGCt                                                             >  3:282498/1‑148 (MQ=255)
                                                                                          ccaGCACGATCCAGCCCATTAC‑‑CGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAGGGCCGCGCAGGCTGAc                                                          >  7:176181/1‑148 (MQ=255)
                                                                                           caGCACGATCCAGCCCATTAC‑‑CGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAGGGCCGCGCAGGCTGACCa                                                        <  4:396149/149‑1 (MQ=255)
                                                                                              cACGATCCAGCCCATTAC‑‑CGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAGGGCCGCGCAGGCTGACCAgg                                                      >  8:92838/1‑148 (MQ=255)
                                                                                              cACGATCCAGCCCATTAC‑‑CGCATAGATGATGATCGACACGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAGGGCCGCGCAGGCTGACCAgg                                                      >  8:397857/1‑148 (MQ=255)
                                                                                                cGATCCAGCCCATTAC‑‑CGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAGGGCCGCGCAGGCTGACCAGGCaa                                                   <  8:176181/149‑1 (MQ=255)
                                                                                                 gATCCAGCCCATTAC‑‑CGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAGGGCCGCGCAGGCTGACCAGGCaaa                                                  <  6:153975/149‑1 (MQ=255)
                                                                                                   tCCAGCCCATTAC‑‑CGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCACCATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAGGGCCGCGCAGGCTGACCAGGCAAAAc                                                <  7:92838/149‑1 (MQ=255)
                                                                                                     cAGCCCATTAC‑‑CGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCCCGCCGAACAGGCTCCATCCCCAAGGGCCGCGCAGGCTGACCAGGCAAAACgg                                              <  2:372656/149‑1 (MQ=255)
                                                                                                        cccATTAC‑‑CGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAGGGCCGCGCAGGCTGACCAGGCGAAACGGCgt                                           <  3:220647/149‑1 (MQ=255)
                                                                                                           aTTAC‑‑CGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAGGGCCGCGCAGGCTGACCAGGCAAAACGGCGTGTa                                        <  7:311429/149‑1 (MQ=255)
                                                                                                                           gatgatCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAGGGCCGCCCAGGCTGACCAGGCAAAACGGCGTGTAGCTGCCGGCAAt                            >  5:371191/1‑147 (MQ=255)
                                                                                                                             tgatCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAGGGCCGCGCAGGCTGACCAGGCAAAACGGCGTGTAGCTGCCGGCAATcagc                        <  6:371191/149‑1 (MQ=255)
                                                                                                                               atCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAGGGCCGCGCAGGCTGACCAGGCAAAACGGCGTGTAGCTGCCGGCAATcagca                       <  5:386532/148‑1 (MQ=255)
                                                                                                                               atCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAGGGCCGCGCAGGCTGACCAGGCAAAACGGCGTGTAGCTGCCGGCAATcagca                       <  1:378867/148‑1 (MQ=255)
                                                                                                                                tCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAGGGCCGCGCAGGCTGACCAGGCAAAACGGCGTGTAGCTGCCGGCAATcag                         >  3:22388/1‑145 (MQ=255)
                                                                                                                                tCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAGGGCCGCGCAGGCTGACCAGGCAAAACGGCGTGTAGCTGCCGGCAATCAGCaga                     <  7:85835/149‑1 (MQ=255)
                                                                                                                                   aCAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAGGGCCGCGCAGGCTGACCAGGCAAAACGGCGTGTAGCTGCCGGCAATCAGCaga                     >  2:69814/1‑146 (MQ=255)
                                                                                                                                   aCAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAGGGCCGCGCAGGCTGACCAGGCAAAACGGCGTGTAGCTGCCGGCAATCAGCaga                     >  2:57870/1‑146 (MQ=255)
                                                                                                                                            gcgcCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAGGGCCGCGCAGGCTGACCAGGCAAAACGGCGTGTAGCTGCCGGCAATCAGCAGATAGATCGACAgg         >  3:113570/1‑149 (MQ=255)
                                                                                                                                            gcgcCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAGGGCCGCGCAGGCTGACCAGGCAAAACGGCGTGTAGCTGCCGGCAATCAGCAGATAGATCGACAgg         >  3:286661/1‑149 (MQ=255)
                                                                                                                                            gcgcCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAGGGCCGCGCAGGCTGACCAGGCAAAACGGCGTGTAGCTGCCGGCAATCAGCAGATAGATCGACAgg         >  4:350177/1‑149 (MQ=255)
                                                                                                                                               ccTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAGGGCCGCGCAGGCTGACCAGGCAAAACGGCGTGTAGCTGCCGGCAATCAGCAGATAGATCGACAggtgg      <  6:339795/149‑1 (MQ=255)
                                                                                                                                               ccTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAGGGCCGCGCAGGCTGACCAGGCAAAACGGCGTGTAGCTGCCGGCAATCAGCAGATAGATCGACAggtgg      <  6:51856/149‑1 (MQ=255)
                                                                                                                                               ccTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAGGGCCGCGCAGGCTGACCAGGCAAAACGGCGTGTAGCTGCCGGCAATCAGCAGATAGATCGACAggtgg      <  1:323085/149‑1 (MQ=255)
                                                                                                                                                cTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAGGGCCGCGCAGGCTGACCAGGCAAAACGGCGTGTAGCTGCCGGCAATCAGCAGATAGATCGACAggtggt     >  2:278693/1‑149 (MQ=255)
                                                                                                                                                   ggAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAGGGCCGCGCAGGCTGACCAGGCAAAACGGCGTGTAGCTGCCGGCAATCAGCAGATAGATCGACAGGTGGTCGa  <  5:23209/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
AAATGATGCCGACGGTGTAGAACACGCCGCCGGCAGCCAGCCAGGCGAAGCCGGCCGTGCCCAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTAC‑‑CGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAGGGCCGCGCAGGCTGACCAGGCAAAACGGCGTGTAGCTGCCGGCAATCAGCAGATAGATCGACAGGTGGTCGA  >  NC_002947/5680772‑5681065

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