Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,420,652 A→T 5.3% intergenic (‑42/+67) PP_3035 ← / ← PP_5554 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,420,6520AT 5.3% 229.6 / 2.5 76intergenic (‑42/+67)PP_3035/PP_5554hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (37/35);  new base T (2/2);  total (39/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.88e-01

AGAGACAAAGGCCGACATAGAAAAAAGGCGACCCGAAGGTCGCCTGGTTTCCCCCAGTAGTACCACCACAGTGCTACCGGGCCGCGTTCAAGGATGTGGTGCACACCACGAAACTACTGTCCTTGAACCCGGGGGTCCAGGCACAAGTGCCTGGGTGCCGCCTTTCTCCACCACAGATAGGGCGGCGGGGTGCCAGTGATAGTCGGGCGACATCACTGGTAATAGCTCTTAACGAGCCATGCGGCGAGTGTAACACCACCGCTCGTCGCTGGGACTGGCAGCTTAAGTGAC  >  NC_002947/3420507‑3420797
                                                                                                                                                 |                                                                                                                                                    
agagACAAAGGCCGACATAGAAAAAAGGCGACCCGAAGGTCGCCTGGTTTCCCCCAGTAGTACCACCACAGTGCTACCGGGCCGCGTTCAAGGATGTGGTGCACACCACGAAACTACTGTCCTTGAACCCGGGGGTCCAGGCACAAGTg                                                                                                                                                  >  7:160046/1‑149 (MQ=255)
agagACAAAGGCCGACATAGAAAAAAGGCGACCCGAAGGTCGCCTGGTTTCCCCCAGTAGTACCACCACAGTGCTACCGGGCCGCGTTCAAGGATGTGGTGCACACCACGAAACTACTGTCCTTGAACCCGGGGGTCCAGGCACAAGTg                                                                                                                                                  >  7:27841/1‑149 (MQ=255)
  agaCAAAGGCCGACATAGAAAAAAGGCGACCCGAAGGTCGCCTGGTTTCCCCCAGTAGTACCACCACAGTGCTACCGGGCCGCGTTCAAGGATGTGGTGCACACCACGAAACTACTGTCCTTGAACCCGGGGGTCCAGGCACAAGTGcg                                                                                                                                                >  8:323141/1‑148 (MQ=255)
      aaaGGCCGACATAGAAAAAAGGCGACCCGAAGGTCGCCTGGTTTCCCCCAGTAGTACCACCACAGTGCTACCGGGCCGCGTTCAAGGATGTGGTGCACACCACGAAACTACTGTCCTTGAACCCGGGGGTCCAGGCACTTGTGCCTgga                                                                                                                                           <  2:222945/149‑2 (MQ=255)
      aaaGGCCGACATAGAAAAAAGGCGACCCGAAGGTCGCCTGGTTTCCCCCAGTAGTACCACCACAGTGCTACCGGGCCGCGTTCAAGGATGTGGTGCACACCACGAAACTACTGTCCTTGAACCCGGGGGTCCAGGCACAAGTGCCTggg                                                                                                                                           >  6:313674/1‑149 (MQ=255)
      aaaGGCCGACATAGAAAAAAGGCGACCCGAAGGTCGCCTGGTTTCCCCCAGTAGTACCACCACAGTGCTACCGGGCCGCGTTCAAGGATGTGGTGCACACCACGAAACTACTGTCCTTGAACCCGGGGGTCCAGGCACAAGTGCCTggg                                                                                                                                           >  2:65387/1‑149 (MQ=255)
       aaGGCCGACATAGAAAAAAGGCGACCCGAAGGTCGCCTGGTTTCCCCCAGTAGTACCACCACAGTGCTACCGGGCCGCGTTCAAGGATGTGGTGCACACCACGAAACTACTGTCCTTGAACCCGGGGGTCCAGGCACAAGTGCCTGGGt                                                                                                                                          <  5:287764/149‑1 (MQ=255)
         ggCCGACATAGAAAAAAGGCGACCCGAAGGTCGCCTGGTTTCCCCCAGTAGTACCACCACAGTGCTACCGGGCCGCGTTCAAGGATGTGGTGCACACCACGAAACTACTGTCCTTGAACCCGGGGGCCCAGGCACTTGTGCCTggaccc                                                                                                                                        <  6:90708/149‑5 (MQ=255)
            cGACATAGAAAAAAGGCGACCCGAAGGTCGCCTGGTTTCCCCCAGTAGTACCACCACAGTGCTACCGGGCCGCGTTCAAGGATGTGGTGCACACCACGAAACTACTGTCCTTGAACCCGGGGGTCCAGGCACAAGTGCCTGGGTgccgc                                                                                                                                     <  8:311876/149‑1 (MQ=255)
                aTAGAAAAAAGGCGACCCGAAGGTCGCCTGGTTTCCCCCAGTAGTACCACCACAGTGCTACCGGGCCGCGTTCAAGGATGTGGTGCACACAACGAAACTACTGTCCTTGAACCCGGGGGTCCAGGCACAAGTGCCTGGGTGCCGCCtt                                                                                                                                  >  7:143648/1‑148 (MQ=255)
                  aGAAAAAAGGCGACCCGAAGGTCGCCTGGTTTCCCCCAGTAGTACCACCACAGTGCTACCGGGCCGCGTTCAAGGATGTGGTGCACACCACGAAACTACTGTCCTTGAACCCGGGGGTCCAGGCACAAGTGCCTGGGTGCCGCCTTtct                                                                                                                               >  5:41397/1‑149 (MQ=255)
                   gAAAAAAGGCGACCCGAAGGTCGCCTGGTTTCCCCCAGTAGTACCACCACAGTGCTACCGGGCCGCGTTCAAGGATGTGGTGCACACCACGAAACTACTGTCCTTGAACCCGGGGGTCCAGGCACAAGTGCCTGGGTGCCGCCTTtctc                                                                                                                              <  2:317384/149‑1 (MQ=255)
                          ggCGACCCGAAGGTCGCCTGGTTTCCCCCAGTAGTACCACCACAGTGCTACCGGGCCGCGTTCAAGGATGTGGTGCACACCACGAAACTACTGTCCTTGAACCCGGGGGTCCAGGCACAAGTGCCTGGGTGCCGCCTTTCTCCACcaca                                                                                                                       <  1:65387/149‑1 (MQ=255)
                           gCGACCCGAAGGTCGCCTGGTTTCCCCCAGTAGTACCACCACAGTGCTACCGGGCCGCGTTCAAGGATGTGGTGCACACCACGAAACTACTGTCCTTGAACCCGGGGGTCCACGCACAAGAGCCGGGGTGCCCCCGTtatcaccca                                                                                                                         >  2:174402/1‑141 (MQ=255)
                            cGACCCGAAGGTCGCCTGGTTTCCCCCAGTAGTACCACCACAGTGCTACCGGGCCGCGTTCAAGGATGTGGTGCACACCACGAAACTACTGTCCTTGAACCCGGGGGTCCAGGCACAAGTGCCTGGGTGCCGCCTTTCTCCACCACAGa                                                                                                                     <  3:355781/149‑1 (MQ=255)
                                  gAAGGTCGCCTGGTTTCCCCCAGTAGTACCACCACAGTGCTACCGGGCCGCGTTCAAGGATGTGGTGCACACCACGAAACTACTGTCCTTGAACCCGGGGGTCCAGGCACAAGTGCCTGGGTGCCGCCTTTCTCCACCACAGATAGggc                                                                                                               >  8:310492/1‑149 (MQ=255)
                                      gTCGCCTGGTTTCCCCCAGTAGTACCACCACAGTGCTACCGGGCCGCGTTCAAGGATGTGGTGCACACCACGAAACTACTGTCCTTGAACCCGGGGGTCCAGGCACAAGTGCCTGGGTGCCGCCTTTCTCCACCACAGATAGGgcggcg                                                                                                           <  6:363765/149‑1 (MQ=255)
                                         gCCTGGTTTCCCCCAGTAGTACCACCACAGTGCTACCGGGCCGCGTTCAAGGATGTGGTGCACACCACGAAACTACTGTCCTTGAACCCGGGGGTCCAGGCACAAGTGCCTGGGTGCCGCCTTTCTCCAGCACAGATAGGGCGGCgggg                                                                                                        >  5:197784/1‑149 (MQ=255)
                                         gCCTGGTTTCCCCCAGTAGTACCACCACAGTGCTACCGGGCCGCGTTCAAGGATGTGGTGCACACCACGAAACTACTGTCCTTGAACCCGGGGGTCCAGGCACAAGGGCATGGGTGCCGCCTTTCTCCACAACAGATAGGGCGGCgggg                                                                                                        >  7:327092/1‑149 (MQ=255)
                                          ccTGGTTTCCCCCAGTAGTACCACCACAGTGCTACCGGGCCGCGTTCAAGGATGTGGTGCACACCACGAAACTACTGTCCTTGAACCCGGGGGTCCAGGCACAAGTGCCTGGGTGCCGCCTTTCTCCACCACAGATAGGGCGGCGGGGt                                                                                                       <  5:247638/149‑1 (MQ=255)
                                           cTGGATTCCCCCAGTAGTACCACCACAGTGCTACCGGGCCGCGTTCAAGGATGTGGTGCACACCACGAAACTACTGTCCTTGAACCCGGGGGTCCAGGCACAAGTGCCTGGGTGCCGCCTTTCTCCACCACAGATAGGGCGGCGGGGTg                                                                                                      >  2:103191/1‑149 (MQ=255)
                                               tttCCCCCAGTAGTACCACCACAGTGCTACCGGGCCGCGTTCAAGGATGTGGTGCACACCACGAAACTACTGTCCTTGAACCCGGGGGTCCAGGCACAAGTGCCTGGGTGCCGCCTTTCTCCACCACAGATAGGGCGGCGGGGTGCCAg                                                                                                  >  5:307278/1‑149 (MQ=255)
                                               tttCCCCCAGTAGTACCACCACAGTGCTACCGGGCCGCGTTCAAGGATGTGGTGCACACCACGAAACTACTGTCCTTGAACCCGGGGGTCCAGGCACAAGTGCCTGGGTGCCGCCTTTCTCCACCACAGATAGGGCGGCGGGGTGCCAg                                                                                                  <  6:116031/149‑1 (MQ=255)
                                                 tCCCCCAGTAGTACCACCACAGTGCTACCGGGCCGCGTTCAAGGATGTGGTGCACACCACGAAACTACTGTCCTTGAACCCGGGGGTCCAGGCACAAGTAGCCTGGGTGCCGCCTTTCTCCACCACAGATAGGGCGGCGGGGTGCCAg                                                                                                  <  1:86525/148‑1 (MQ=255)
                                                    cccAGTAGTACCACCACAGTGCTACCGGGCCGCGTTCAAGGATGTGGTGCACACCACGAAACTACTGTCCTTGAACCCGGGGGTCCAGGCACAAGTGCCTGGGTGCCGCCTTTCTCCACCACAGATAGGGCGGCGGGGTGCCAGTGa                                                                                               >  6:185780/1‑147 (MQ=255)
                                                     ccAGTAGTACCACCACAGTGCTACCGGGCCGCGTTCAAGGATGTGGTGCACACCACGAAACTACTGTCCTTGAACCCGGGGGTCCAGGCACAAGTGCCTGGGTGCCGCCTTTCTCCACCACAGATAGGGCGGCGGGGTGCCAGTGATAg                                                                                            >  3:226771/1‑149 (MQ=255)
                                                     ccAGTAGTACCACCACAGTGCTACCGGGCCGCGTTCAAGGATGTGGTGCACACCACGAAACTACTGTCCTTGAACCCGGGGGTCCAGGCACAAGTGCCTGGGTGCCGCCTTTCTCCACCACAGATAGGGCGGCGGGGTGCCAGTGATAg                                                                                            >  4:158744/1‑149 (MQ=255)
                                                         tagtaCCACCACAGTGCTACCGGGCCGCGTTCAAGGATGTGGTGCACACCACGAAACTACTGTCCTTGAACCCGGGGGTCCAGGCACAAGTGCCTGGGTGCCGCCTTTCTCCACCACAGATAGGGCGGCGGGGTGCCAGTGATAGTCgg                                                                                        <  8:160046/149‑1 (MQ=255)
                                                           gtaCCACCACAGTGCTACCGGGCCGCGTTCAAGGATGTGGTGCACACCACGAAACTACTGTCCTTGAACCCGGGGGTCCAGGCACAAGTGCCTGGGTGCCGCCTTTCTCCACCACAGATAGGGCGGCGGGGTGCCAGTGATAGTCGGGc                                                                                      <  6:307278/149‑1 (MQ=255)
                                                           gtaCCACCACAGTGCTACCGGGCCGCGTTCAAGGATGTGGTGCACACCACGAAACTACTGTCCTTGAACCCGGGGGTCCAGGCACAAGTGCCTGGGTGCCGCCTTTCTCCACCACAGATAGGGCGGCGGGGTGCCAGTGATAGTCGGGc                                                                                      <  2:67826/149‑1 (MQ=255)
                                                                  cacaGTGCTACCGGGCCGCGTTCAAGGATGTGGTGCACACCACGAAACTACTGTCCTTGAACCCGGGGGTCCAGGCACAAGTGCCTGGGTGCCGCCTTTCTCCACCACAGATAGGGCGGCGGGGTGCCAGTGATAGTCGGGCGACATCa                                                                               >  2:161791/1‑149 (MQ=255)
                                                                         cTACCGGGCCGCGTTCAAGGATGTGGTGCACACCACGAAACTACTGTCCTTGAACCCGGGGGTCCAGGCACAAGTGCCTGGGTGCCGCCTTTCTCCACCACAGATAGGGCGGCGGGGTGCCAGTGATAGTCGGGCGACATCACTGGTa                                                                         <  1:161791/148‑1 (MQ=255)
                                                                            ccGGGCCGCGTTCAAGGATGTGGTGCACACCACGAAACTACTGTCCTTGAACCCGGGGGTCCAGGCACAAGTGCCTGGGTGCCGCCTTTCTCCCCCACAGATAGGGCGGCGGGGTGCCAGTGATAGTCGGGCGACATCACTGGTAATAg                                                                    >  2:2726/1‑149 (MQ=255)
                                                                                gCCGCGTTCAAGGATGTGGTGCACACCACGAAACTACTGTCCTTGAACCCGGGGGTCCAGGCACAAGTGCCTGGGTGCCGCCTTTCTCCACCACAGATAGGGCGGCGGGGTGCCAGTGATAGTCGGGCGACATCACTGGTAATAGctct                                                                <  8:27841/149‑1 (MQ=255)
                                                                                      ttCAAGGATGTGGTGCACACCACGAAACTACTGTCCTTGAACCCGGGGGTCCAGGCACAAGTGCCTGGGTGCCGCCTTTCTCCACCACAGATAGGGCGGCGGGGTGCCAGTGATAGTCGGGCGACATCACTGGTAATAGCTCTTAACGa                                                          <  1:174402/149‑1 (MQ=255)
                                                                                         aaGGATGTGGTGCACACCACGAAACTACTGTCCTTGAACCCGGGGGTCCAGGCACAAGTGCCTGGGTGCCGCCTTTCTCCACCACAGATAGGGCGGCGGGGTGCCAGTGATAGTCGGGCGACATCACTGGTAATAGCTCTTAACGAGc                                                        <  8:233647/148‑1 (MQ=255)
                                                                                         aaGGATGTGGTGCACACCACGAAACTACTGTCCTTGAACCCGGGGGTCCAGGCACAAGTGCCTGGGTGCCGCCTTTCTCCACCACAGATAGGGCCGCCGGGTGCCAGTGATAGTCGGGCGACATCACTGGTAATAGCTCTTAACGAGc                                                        >  2:103794/1‑148 (MQ=255)
                                                                                           ggATGTGGTGCACACCACGAAACTACTGTCCTTGAACCCGGGGGTCCAGGCACAAGTGCCTGGGTGCCGCCTTTCTCCACCACAGATAGGGCGGCGGGGTGCCAGTGATAGTCGGGCGACATCACTGGTAATAGCTCTTAACGAGCCAt                                                     >  2:195039/1‑149 (MQ=255)
                                                                                                tggtgCACACCACGAAACTACTGTCCTTGAACCCGGGGGTCCAGGCACAAGTGCCTGGGTGCCGCCTTTCTCCACCACAGATAGGGCGGCGGGGTGCCAGTGATAGTCGGGCGACATCACTGGTAATAGCTCTTAACGAGCCATgcggc                                                <  3:52099/149‑1 (MQ=255)
                                                                                                tggtgCACACCACGAAACTACTGTCCTTGAACCCGGGGGTCCAGGCACAAGTGCCTGGGTGCCGCCTTTCTCCACCACAGATAGGGCGGCGGGGTGCCAGTGATAGTCGGGCGACATCACTGGTAATAGCTCTTAACGAGCCATgcggc                                                <  7:310492/149‑1 (MQ=255)
                                                                                                 ggtgCACACCACGAAACTACTGTCCTTGAACCCGGGGGTCCAGGCACAAGTGCCTGGGTGCCGCCTTTCTCCACCACAGATAGGGCGGCGGGGTGCCAGTGATAGTCGGGCGACATCACTGGTAATAGCTCTTAACGAGCCATtcggc                                                >  3:124976/1‑148 (MQ=255)
                                                                                                  gtgCACACCACGAAACTACTGTCCTTGAACCCGGGGGTCCAGGCACAAGTGCCTGGGTGCCGCCTTTCTCCACCACAGATAGGGCGGCGGGGTGCCAGTGATAGTCGGGCGACATCACTGGTAATAGCTCTTAACGAGCCATgcggc                                                >  5:283135/1‑147 (MQ=255)
                                                                                                  gtgCACACCACGAAACTACTGGCCTTGAACCCGGGGGTCCAGGCACAAGTGCCTGGGTGCCGCCTTTCTCCACCACAGATAGGGCGGCGGGGTGCCAGTGATAGTCGGGCGACATCACTGGTAATAGCTCTTAACGAGCCATGCGGCGa                                              <  4:226771/149‑1 (MQ=255)
                                                                                                   tgCACACCACGAAACTACTGTCCTTGAACCCGGGGGTCCAGGCACAAGTGCCTGGGTGCCGCCTTTCTCCACCACAGATAGGGCGGCGGGGTGCCAGTGATAGTCGGGCGACATCACTGGTACATAGCTCTTAACGAGCCATGCGGCGa                                              <  1:195039/149‑1 (MQ=255)
                                                                                                      acacCACGAAACTACTGTCCTTGAACCCGGGGGTCCAGGCACAAGTGCCTGGGTGCCGCCTTTCTCCACCACAGATAGGGCGGCGGGGTGCCAGTGATAGTCGGGCGACATCACTGGTAATAGCTCTTAACGAGCCATGCGGCGAgtgt                                          <  2:250258/149‑1 (MQ=255)
                                                                                                        accacGAAACTACTGTCCTTGAACCCGGGGGTCCAGGCACAAGTGCCTGGGTGCCGCCTTTCTCCACCACAGATAGGGCGGCGGGGTGCCAGTGATAGTCGGGCGACATCACTGGTAATAGCTCTTAACGAGCCATGCGGCGAGTGTaa                                        >  7:247491/1‑149 (MQ=255)
                                                                                                          cacGAAACTACTGTCCTTGAACCCGGGGGTCCAGGCACAAGTGCCTGGGTGCCGCCTTTCTCCACCACAGATAGGGCGGCGGGGTGCCAGTGATAGTCGGGCGACATCACTGGTaa                                                                       <  3:232476/116‑1 (MQ=255)
                                                                                                          cacGAAACTACTGTCCTTGAACCCGGGGGTCCAGGCACAAGTGCCTGGGTGCCGCCTTTCTCCACCACAGATAGGGCGGCCGGGTGCCAGTGATAGTCGGGCGACATCACTGGTaa                                                                       >  4:232476/1‑116 (MQ=255)
                                                                                                            cGAAACTACTGTCCTTGAACCCGGGGGTCCAGGCACAAGTGCCTGGGTGCCGCCTTTCTCCACCACAGATAGGGCGGCGGGGTGCCAGTGATTGTCGGGCGACATCACTGGTAATAGCTCTTAACGAGCCATGCGGCGAGTGTAAcacc                                    <  1:103191/149‑1 (MQ=255)
                                                                                                                  tactGTCCTTGAACCCGGGGGTCCAGGCACAAGTGCCTGGGTGCCGCCTTTCTCCACCACAGATAGGGCGGCGGGGTGCCAGTGATAGTCGGGCGACATCACTGGTAATAGCTCTTAACGAGCCATGCGGCGAGTGTAac                                       >  8:22707/1‑140 (MQ=255)
                                                                                                                  tactGTCCTTGAACCCGGGGGTCCAGGCACAAGTGCCTGGGTGCCGCCTTTCTCCACCACAGATAGGGCGGCGGGGTGCCAGTGATAGTCGGGCGACATCACTGGTAATAGCTCTTAACGAGCCATGCGGCGAGTGTAac                                       <  7:22707/140‑1 (MQ=255)
                                                                                                                  tactGTCCTTGAACCCGGGGGTCCAGGCACAAGTGCCTGGGTGCCGCCTTTCTCCACCACAGATAGGGCGGCGGGGTGCCAGTGATAGTCGGGCGACATCACTGGTAATAGCTCTTAACGAGCCATGCGGCGAGTGTAACACCACCGCt                              >  5:172902/1‑149 (MQ=255)
                                                                                                                     tGTCCTTGAACCCGGGGGTCCAGGCACAAGTGCCTGGGTGCCGCCTTTCTCCACCACAGATAGGGCGGCGGGGTGCCAGTGATAGTCGGGCGACATCACTGGTAATAGCTCTTAACGAGCCATGCGGCGAGTGTAACACCACCGCtcgt                           >  5:39653/1‑149 (MQ=255)
                                                                                                                     tGTCCTTGAACCCGGGGGTCCAGGCACAAGTGCCTGGGTGCCGCCTTTCTCCACCACAGATAGGGCGGCGGGGTGCCAGTGATAGTCGGGCGACATCACTGGTAATAGCTCTTAACGAGCCATGCGGCGAGTGTAACACCACCGCtcgt                           >  6:313451/1‑149 (MQ=255)
                                                                                                                       tCCTTGAACCCGGGGGTCCAGGCACAAGTGCCTGGGTGCCGCCTTTCTCCACCACAGATAGGGCGGCGGGGTGCCAGTGATAGTCGGGCGACATCACTGGTAATAGCTCTTAACGAGCCATgcggc                                                <  1:255018/126‑1 (MQ=255)
                                                                                                                       tCCTTGAACCCGGGGGTCCAGGCACAAGTGCCTGGGTGCCGCCTTTCTCCACCACAGATAGGGCGGCGGGGTGCCAGTGATAGTCGGGCGACATCACTGGTAATAGCTCTTAACGAGCCATgcggc                                                >  2:255018/1‑126 (MQ=255)
                                                                                                                         cTTGAACCCGGGGGTCCAGGCACAAGTGCCGTGGGTGCCGCCTTTCTCCACCACAGATAGGGCGGCGGGGTGCCAGTGATAGTCGGGCGACATCACTGGTAATAGCTCTTAACGAGCCATGCGGCGAGTGTAACACCACCGCTCGTCGc                        >  1:223331/1‑149 (MQ=255)
                                                                                                                         cTTGAACCCGGGGGTCCAGGCACAAGTGCCTGGGTGCCGCCTTTCTCCACCACAGATAGGGCGGCGGGGTGCCAGTGATAGTCGGGCGACATCACTGGTAATAGCTCTTAACGAGCCATGCGGCGAGTGTAACACCACCGCTCGTCGCt                       >  4:240924/1‑149 (MQ=255)
                                                                                                                            gAACCCGGGGGTCCAGGCACAAGTGCCTGGGTGTTGCCTTTCTCCTGCTCAGATAGGGCGGCGGTGTGCCAGTGATAGTCGGGCGACATCACTGGTAATATCTCTTAAGGTGCGATGCGGCGAGTGTAACACCACAGCTCGTCGCTg                      <  4:124976/147‑1 (MQ=255)
                                                                                                                            gAACCCGGGGGTCCAGGCACAAGTGCCTGGGTGCCGCCTTTCTCCACCACAGATAGGGCGGCGGGGTGCCAGTGATAGTCGGGCGACATCACTGGTAATAGCTCTTAACGAGCCATGCGGCGAGTGTAACACCACCGCTCGTCGCTgg                     <  7:223168/148‑1 (MQ=255)
                                                                                                                            gAACCCGGGGGTCCAGGCACAAGTGCCTGGGTGCCGCCTTTCTCCACCACAGATAGGGCGGCGGGGTGCCAGTGATAGTCGGGCGACATCACTGGTAATAGCTCTTAACGAGCCATGCGGCGAGTGTAACACCACCGCTCGTCGCTg                      <  6:172902/147‑1 (MQ=255)
                                                                                                                            gAACCCGGGGGTCCAGGCACAAGTGCCTGGGTGCCGCCTTTCTCCACCACAGATAGGGCGGCGGGGTGCCAGTGATAGTCGGGCGACATCACTGGTAATAGCTCTTAACGAGCCATGCGGCGAGTGTAACACCACCGCTCGTCGCTg                      <  7:205140/147‑1 (MQ=255)
                                                                                                                            gAACCCGGGGGTCCAGGCACAAGTGCCTGGGTGCCGCCTTTCTCCACCACAGATAGGGCGGCGGGGTGCCAGTGATAGTCGGGCGACATCACTGGTAATAGCTCTTAACGAGCCATGCGGCGAGTGTAACACCACCGCTCGTCGCTg                      <  6:41397/147‑1 (MQ=255)
                                                                                                                              aCCCGGGGGTCCAGGCACAAGTGCCTGGGTGCCGCCTTTCTCCACCACAGATAGGGCGGCAGGGTGCCAGTGATAGTCGGGCGACATCACTGGTAATAGCTCTTAACGAGCCATGCGGCGAGTGTAACACCACCGCTCGTCGCTGGGa                   <  8:143648/148‑1 (MQ=255)
                                                                                                                                ccGGGGGTCCAGGCACAAGTGCCTGGGTGCCGCCTTTCTCCACCACAGATAGGGCGGCGGGGTGCCAGTGATAGTCGGGCGACATCACTGGTAATAGCTCTTAACGAGCCATGCGGCGAGTGTAACACCACCGCTCGTCGCTGGGACTg                >  1:175430/1‑149 (MQ=255)
                                                                                                                                  gggggTCCAGGCACAAGTGCCTGGGTGCCGCCTTTCTCCACCACAGATAGGGCGGCGGGGTGCCAGTGATAGTCGGGCGACATCACTGGTAATAGCTCTTAACGAGCCATGCGGCGAGTGTAACACCACCGCTCGTCGCTGGGACTGGc              <  3:158744/149‑1 (MQ=255)
                                                                                                                                   ggggTCCACGCACAAGTGCCTGGGTGCCGCCTTTCTCCACCACAGATAGGGCGGCGGGGTGCCAGTGATAGTCGGGCGACATCACTGGTAATAGCTCTTAACGAGCCATGCGGCGAGTGTAACACCACCGCTCGTCGCTGGGACTGGCa             <  4:80210/149‑1 (MQ=255)
                                                                                                                                    gggTCCAGGCACAAGTGCCTGGGTGCCGCCTTTCTCCACCACAGATAGGGGGGCGGGGTGCCAGTGATAGTCGGGCGACATCACTGGAAATAGCTCTTACCGAGCCATGAGGCGAGTGTAACACCACCCCGGGTCCCTGGGCCTGGAAg            >  4:193606/1‑149 (MQ=255)
                                                                                                                                    gggTCCAGGCACAAGTGCCTGGGTGCCGCCTTTCTCCACCACAGATAGGGCGGCGGGGTGCCAGTGATAGTCGGGCGACATCACTGGTAATAGCTCTTAACGAGCCATGCGGCGAGTGTAACACCACCGCTCGTCGCTGGGACTGGCAg            <  6:197784/149‑1 (MQ=255)
                                                                                                                                      gTCCAGGCACAAGTGCCTGGGTGCCGCCTTTCTCCACCACAGATAGGGCGGCGGGGTGCCAGTGATAGTCGGGCGACATCACTGGTAATAGCTCTTAACGAGCCATGCGGCGAGTGTAACACCACCGCTCGTCGCTGGGACTGGCAGCt          >  1:15728/1‑149 (MQ=255)
                                                                                                                                      acccAGGCACTTGTGCCTGGGTGCCGCCTTTCTCCACCACAGATAGGGCGGCGGGGTGCCAGTGATAGTCGGGCGACATCACTGGTAATAGCTCTTAACGAGCCATGCGGCGAGTGTAACACCACCGCTCGTCGCTGGGACTGGCAGCt          >  8:56074/3‑149 (MQ=255)
                                                                                                                                        ccAGGCACTTGTGCCTGGGTGCCGCCTTTCTCCACCACAGATAGGGCGGCGGGGTGCCAGTGATAGTCGGGCGACATCACTGGTAATAGCTCTTAACGAGCCATGCGGCGAGTGTAACACCAACGCTCGTCGGTGGGACTGGCAGCtt         >  6:378830/1‑148 (MQ=255)
                                                                                                                                        ccAGGCACAAGTGCCTGGGTGCCGCCTTTCTCCACCACAGATAGGGCGGCGGGGTGCCAGTGATAGTCGGGCGACATCACTGGTAATAGCTCTTAACGAGCCATGCGGCGAGTGTAACACCACCGCTCGTCGCTGGGACTGGCAGCtt         >  2:185085/1‑148 (MQ=255)
                                                                                                                                          aGGCACAAGTGCCTGGGTGCCGCCTTTCTCCACCACAGATAGGGCGGCGGGGTGCCAGTGATAGTCGGGCGACATCACTGGTAATAGCTCTTAACGAGCCATGCGGCGAGTGTAACACCACCGCTCGTCGCTGGGACTGGCAGCTTAAg      >  4:119910/1‑149 (MQ=255)
                                                                                                                                           ggCACAAGTGCCTGGGTGCCGCCTTTCTCCACCACAGATAGGGCGGCGGGGTGCCAGTGATAGTCGGGCGACATCACTGGTAATAGCTCTTAACGAGCCATGCGGCGAGTGTAACACCACCGCTCGTCGCTGGGACTGGCAGCTTAAGt     <  5:313674/149‑1 (MQ=255)
                                                                                                                                              acaAGTGCCTGGGTGCCGCCTTTCTCCACCACAGATAGGGCGGCGGGGTGCCAGTGATAGTCGGGCGACATCACTGGTAATAGCTCTTAACGAGCCATGCGGCGAGTGTAACACCACCGCTCGTCGCTGGGACTGGCAGCTTAAGTGAc  >  3:52037/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
AGAGACAAAGGCCGACATAGAAAAAAGGCGACCCGAAGGTCGCCTGGTTTCCCCCAGTAGTACCACCACAGTGCTACCGGGCCGCGTTCAAGGATGTGGTGCACACCACGAAACTACTGTCCTTGAACCCGGGGGTCCAGGCACAAGTGCCTGGGTGCCGCCTTTCTCCACCACAGATAGGGCGGCGGGGTGCCAGTGATAGTCGGGCGACATCACTGGTAATAGCTCTTAACGAGCCATGCGGCGAGTGTAACACCACCGCTCGTCGCTGGGACTGGCAGCTTAAGTGAC  >  NC_002947/3420507‑3420797

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: