Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,722,805 A→T 7.6% Y21F (TAT→TTT) ‡ PP_3288 → universal stress protein family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,722,8050AT 7.6% 234.4 / 7.8 79Y21F (TAT→TTT) ‡PP_3288universal stress protein family protein
Reads supporting (aligned to +/- strand):  ref base A (36/37);  new base T (4/2);  total (40/39)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01

TGCCTTGAGTTGATAAACATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAG  >  NC_002947/3722662‑3722944
                                                                                                                                               |                                                                                                                                           
tGCCTTGAGTTGATAAACATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGAcc                                                                                                                                        <  1:245099/149‑1 (MQ=255)
 gCCTTGAGTTGATAAACATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCg                                                                                                                                       >  8:115616/1‑149 (MQ=255)
    ttGAGTTGATAAACATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCATAGACCgggg                                                                                                                                    >  2:213379/1‑148 (MQ=255)
      gAGTTGATAAACATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGgc                                                                                                                                    <  6:187712/147‑1 (MQ=255)
      gAGTTGATAAACATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCATAGACCgggg                                                                                                                                    <  1:213379/147‑2 (MQ=255)
         ttGATAAACATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGctgctg                                                                                                                               >  3:293224/1‑149 (MQ=255)
                 cATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGCCCGGGCTGCTGCGCTggcc                                                                                                                       >  8:232865/1‑149 (MQ=255)
                      aGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTggccaaggc                                                                                                                  <  7:115616/149‑1 (MQ=255)
                       gTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTggccaaggcc                                                                                                                 <  3:325665/149‑1 (MQ=255)
                           gCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGccaaggccaagg                                                                                                             <  8:289236/149‑1 (MQ=255)
                            cAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGccaaggccaagg                                                                                                             >  4:350409/1‑148 (MQ=255)
                               cGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAAt                                                                                                         <  5:182477/149‑1 (MQ=255)
                                     gCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCt                                                                                                   >  6:123519/1‑149 (MQ=255)
                                      cTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTg                                                                                                  >  2:149739/1‑149 (MQ=255)
                                        ggCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGca                                                                                                >  6:267607/1‑149 (MQ=255)
                                         gCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGc                                                                                                 <  5:123519/147‑1 (MQ=255)
                                            gATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACAtc                                                                                            >  8:57891/1‑149 (MQ=255)
                                                ctctACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACAtcgtcg                                                                                        >  8:250115/1‑149 (MQ=255)
                                                ctctACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACAtcgtc                                                                                         <  8:325840/148‑1 (MQ=255)
                                                 tctACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGgccaaggccaag                                                                                                              >  5:159707/1‑126 (MQ=255)
                                                 tctACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGgccaaggccaag                                                                                                              <  6:159707/126‑1 (MQ=255)
                                                   tacctacTTATCGCAACCTCAGAGGTGTGCCATGGAAACTCATCAAAGACTGCTGCTGATAGCATCCCCACTGATGCGACGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCtt                                                                                     >  5:257551/1‑149 (MQ=255)
                                                      ctacGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTAGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCg                                                                                  >  6:185764/1‑149 (MQ=255)
                                                       tacGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGa                                                                                 >  5:148998/1‑149 (MQ=255)
                                                         cGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGAc                                                                                >  7:209245/1‑148 (MQ=255)
                                                            aTCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACg                                                                            >  1:192656/1‑149 (MQ=255)
                                                              cGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTg                                                                          >  4:314955/1‑149 (MQ=255)
                                                                    cTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAggg                                                                    <  3:350409/149‑1 (MQ=255)
                                                                     tCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTg                                                                          <  8:205440/142‑1 (MQ=255)
                                                                     tCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTg                                                                          >  7:205440/1‑142 (MQ=255)
                                                                     tCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCACCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGc                                                                   >  7:77827/1‑149 (MQ=255)
                                                                       agagTTGTGCAATGGAAACTCCGCAAAGTCTGCTTCTGATCGCCTCCCCCCTGATGCGCCGAACCGCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGc                                                                   <  1:355998/147‑1 (MQ=255)
                                                                              tgCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCAcgg                                                          <  4:293224/149‑1 (MQ=255)
                                                                                ccATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCAcggcg                                                        >  3:188068/1‑149 (MQ=255)
                                                                                ccATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCAcggcg                                                        >  2:355991/1‑149 (MQ=255)
                                                                                 cATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCAcggcgg                                                       <  5:267607/149‑1 (MQ=255)
                                                                                  aTGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCggg                                                      >  6:135428/1‑149 (MQ=255)
                                                                                   tGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGc                                                               <  2:124963/139‑1 (MQ=255)
                                                                                   tGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGc                                                               >  1:124963/1‑139 (MQ=255)
                                                                                     gAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCt                                                   <  1:355991/149‑1 (MQ=255)
                                                                                       aaCTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTgg                                                 >  4:343500/1‑149 (MQ=255)
                                                                                        aCTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGt                                                <  3:314955/149‑1 (MQ=255)
                                                                                          tCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTg                                               >  3:275297/1‑148 (MQ=255)
                                                                                            cGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGaa                                             <  6:148998/148‑1 (MQ=255)
                                                                                            cGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGa                                              <  7:57891/147‑1 (MQ=255)
                                                                                             gCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATg                                           <  7:250115/149‑1 (MQ=255)
                                                                                                    ctgctgctgATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTCATTCACCAGgc                                    <  8:209245/149‑1 (MQ=255)
                                                                                                     tgctgctgATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGt                                                <  2:293501/136‑1 (MQ=255)
                                                                                                     tgctgctgATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCCCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGt                                                >  1:293501/1‑136 (MQ=255)
                                                                                                       ctgctgATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCt                                 >  6:361718/1‑149 (MQ=255)
                                                                                                        tgctgATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGTAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTg                                <  4:188068/149‑1 (MQ=255)
                                                                                                        tgctgATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTg                                >  6:225776/1‑149 (MQ=255)
                                                                                                          ctgATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTgg                               <  7:232865/148‑1 (MQ=255)
                                                                                                             aTCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGc                                                               >  1:42144/1‑113 (MQ=255)
                                                                                                             aTCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGc                                                               <  2:42144/113‑1 (MQ=255)
                                                                                                             aTCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGt                           <  3:343500/149‑1 (MQ=255)
                                                                                                              tCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTc                          <  7:69449/149‑1 (MQ=255)
                                                                                                              tCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTc                          <  5:361718/149‑1 (MQ=255)
                                                                                                              tCGCCGCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTc                          <  8:270960/149‑1 (MQ=255)
                                                                                                                gCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTAGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCa                         <  5:185764/148‑1 (MQ=255)
                                                                                                                     cccccTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCGTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGcgcg                   >  5:301129/1‑149 (MQ=255)
                                                                                                                        ccTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTgg                                                 >  4:239387/1‑116 (MQ=255)
                                                                                                                        ccTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTgg                                                 <  3:239387/116‑1 (MQ=255)
                                                                                                                           gATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAAt                                            >  4:187380/1‑118 (MQ=255)
                                                                                                                               cgcCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGt         <  5:135428/149‑1 (MQ=255)
                                                                                                                                gcCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAg                 >  5:158551/1‑140 (MQ=255)
                                                                                                                                gcCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAg                 <  6:158551/140‑1 (MQ=255)
                                                                                                                                gcCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTg        <  1:149739/149‑1 (MQ=255)
                                                                                                                                 ccGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATg                                                                                                        <  7:260891/52‑1 (MQ=255)
                                                                                                                                 ccGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATg                                                                                                        >  8:260891/1‑52 (MQ=255)
                                                                                                                                 ccGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTg        >  8:69799/1‑148 (MQ=255)
                                                                                                                                     agcccGGTCATAGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTAc           >  7:289772/3‑141 (MQ=255)
                                                                                                                                     agcccGGTCATAGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTAc           <  8:289772/139‑1 (MQ=255)
                                                                                                                                     agcccGGTCATAGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTgcag     >  7:50584/3‑147 (MQ=255)
                                                                                                                                      gcccGGTCATAGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGCGGGTTACGTgcagcag  >  3:266391/2‑149 (MQ=255)
                                                                                                                                         cGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGt                           <  2:326471/121‑1 (MQ=255)
                                                                                                                                         cGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGt                           >  1:326471/1‑121 (MQ=255)
                                                                                                                                              tATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTgcagc    <  2:97817/139‑1 (MQ=255)
                                                                                                                                              tATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTgcagc    >  1:97817/1‑139 (MQ=255)
                                                                                                                                               |                                                                                                                                           
TGCCTTGAGTTGATAAACATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAG  >  NC_002947/3722662‑3722944

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: