Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,458,892:1 +A 6.6% intergenic (+17/‑45) leuS → / → PP_4795 leucine‑‑tRNA ligase/lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,458,8921.A 6.6% 238.5 / 11.2 61intergenic (+17/‑45)leuS/PP_4795leucine‑‑tRNA ligase/lipoprotein
Reads supporting (aligned to +/- strand):  ref base . (29/28);  new base A (2/2);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.75e-01

CAGCCGCGAAGAAGTCGAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCTGAGCGTGAAG  >  NC_002947/5458752‑5459035
                                                                                                                                             |                                                                                                                                               
cAGCCGCGAAGAAGTCGAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGc                                                                                                                                                  <  1:170495/141‑1 (MQ=255)
cAGCCGCGAAGAAGTCGAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGc                                                                                                                                                  >  2:170495/1‑141 (MQ=255)
  gCCGCGAAGAAGTCGAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCg                                                                                                                                       >  7:116826/1‑149 (MQ=255)
       gaagaaGTCGAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTgg                                                                                                                                  <  4:167828/149‑1 (MQ=255)
        aagaagTCGAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCAATCGATGCGGGTggg                                                                                                                                 <  6:316847/149‑1 (MQ=255)
        aagaagTCGAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTggg                                                                                                                                 >  6:319622/1‑149 (MQ=255)
             gTCGAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGaa                                                                                                                            >  2:8318/1‑149 (MQ=255)
             gTCGAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGCCAAGCTGGTCAACATCGTCGGCAACTGATGAAAACGCCCACCCGCATCGATTGCTGGTGGGCGGaa                                                                                                                            >  6:365503/1‑149 (MQ=255)
              tCGAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCAATCGATGCGGGTgg                                                                                                                                  <  8:81590/142‑1 (MQ=255)
              tCGAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCAATCGATGCGGGTgg                                                                                                                                  >  7:81590/1‑142 (MQ=255)
                gAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGATCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAg                                                                                                                         <  7:220436/149‑1 (MQ=255)
                         gCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAg                                                                                                                >  7:302783/1‑149 (MQ=255)
                            cGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTTATGAAAACGCCCACCCGCTTCGATTGCGGTTGGGCGGAAc                                                                                                                           <  3:115110/135‑1 (MQ=255)
                            cGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGa                                                                                                             <  7:242563/149‑1 (MQ=255)
                            cGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGACAACTCATGAAAACGCCCACCCGCATCGATTGCTGGTGGGCGGAAc                                                                                                                           >  4:115110/1‑135 (MQ=255)
                             gcagcaACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAg                                                                                                            <  1:297432/149‑1 (MQ=255)
                             gcagcaACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAg                                                                                                            <  5:161615/149‑1 (MQ=255)
                              cagcaACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAATACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGTAGc                                                                                                           <  5:365503/149‑1 (MQ=255)
                              cagcaACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGc                                                                                                           <  3:246620/149‑1 (MQ=255)
                                gcaACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGcaa                                                                                                         <  6:318381/149‑1 (MQ=255)
                                   aCGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGcaacaa                                                                                                      >  6:253836/1‑149 (MQ=255)
                                    cGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCaacaac                                                                                                     <  1:8318/149‑1 (MQ=255)
                                     gagaATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAacaaca                                                                                                    <  1:152103/149‑1 (MQ=255)
                                        aaTGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATg                                                                                                  >  6:360595/1‑148 (MQ=255)
                                             ccTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCaaa                                                                                            >  1:80462/1‑149 (MQ=255)
                                               tGCGCTTCATCGATGGCCTGACCATCCGCATGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGTCCCAATAGGGAGCAACAAGATGATCAAAc                                                                                           >  6:287657/1‑148 (MQ=255)
                                                       aTCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGACCAAAAGGGAGCAACAACATGATCAAACGCAATTTGc                                                                                  >  5:4269/1‑149 (MQ=255)
                                                        tCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCt                                                                                 <  8:57253/149‑1 (MQ=255)
                                                          gATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAACCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTg                                                                                <  5:319622/148‑1 (MQ=255)
                                                               ccTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATc                                                                          >  2:191881/1‑149 (MQ=255)
                                                                cTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCg                                                                         >  6:277289/1‑149 (MQ=255)
                                                                  gACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTg                                                                                >  8:112343/1‑140 (MQ=255)
                                                                  gACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTg                                                                                <  7:112343/140‑1 (MQ=255)
                                                                     cATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCCTGCGTTGCGGGTGGGCGGAACAGGGCCAAAAGGGAGGCACAACATGTAGAAACGCAATTTGGTGGTAATCGGccgg                                                                    >  2:350489/1‑147 (MQ=255)
                                                                         cGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCgg                                                                <  8:302783/149‑1 (MQ=255)
                                                                          gCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGt                                                               >  6:301293/1‑149 (MQ=255)
                                                                            aaGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCa                                                             <  5:277289/149‑1 (MQ=255)
                                                                                 cATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCt                                                         >  4:88700/1‑148 (MQ=255)
                                                                                  aTCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCa                                                       <  5:360595/149‑1 (MQ=255)
                                                                                  aTCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCa                                                       <  6:4269/149‑1 (MQ=255)
                                                                                                       gTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCt                                    <  2:40169/147‑1 (MQ=255)
                                                                                                          aaCATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGcgcg                               >  4:208647/1‑149 (MQ=255)
                                                                                                                      aaCTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGTCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCa                   <  5:287657/149‑1 (MQ=255)
                                                                                                                      aaCTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCa                   >  7:173726/1‑149 (MQ=255)
                                                                                                                          gatgaAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCAc                           <  2:274505/137‑1 (MQ=255)
                                                                                                                          gatgaAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCAc                           >  1:274505/1‑137 (MQ=255)
                                                                                                                          gatgaAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCAcc                 >  6:277204/1‑147 (MQ=255)
                                                                                                                           atgaAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAg              >  1:23264/1‑149 (MQ=255)
                                                                                                                           atgaAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAg              >  7:274710/1‑149 (MQ=255)
                                                                                                                           atgaAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAg              >  3:69390/1‑149 (MQ=255)
                                                                                                                            tgaAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGc             <  5:253836/149‑1 (MQ=255)
                                                                                                                              aaaaCGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCt            <  5:301293/148‑1 (MQ=255)
                                                                                                                               aaaCGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCTGa          <  3:208647/149‑1 (MQ=255)
                                                                                                                                  cGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCTGAGCg       <  8:274710/149‑1 (MQ=255)
                                                                                                                                   gCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCTGAGCGt      >  3:29976/1‑149 (MQ=255)
                                                                                                                                     ccACCCGCAATCGATGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGATAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCTGAGCGTg     >  7:8854/1‑148 (MQ=255)
                                                                                                                                     ccACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCTGAGCGTg     >  4:53248/1‑148 (MQ=255)
                                                                                                                                     ccACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCTGAGCGTg     >  7:135302/1‑148 (MQ=255)
                                                                                                                                      caccgccATCTATTGGGGTTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGGCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCTGAGCGTGaa   <  1:350489/143‑1 (MQ=255)
                                                                                                                                      cACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCTGAGCGTGaa   <  2:23264/149‑1 (MQ=255)
                                                                                                                                       aCCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCTGAGCGTGAAg  <  5:308912/149‑1 (MQ=255)
                                                                                                                                       aCCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCTGAGCGTGAAg  >  6:308912/1‑149 (MQ=255)
                                                                                                                                           gCATCGATTGCGGGTGGGCGGAACAGGCCCAAATGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTc                                         <  5:123829/106‑1 (MQ=255)
                                                                                                                                           gCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTc                                         >  6:123829/1‑106 (MQ=255)
                                                                                                                                             |                                                                                                                                               
CAGCCGCGAAGAAGTCGAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCTGAGCGTGAAG  >  NC_002947/5458752‑5459035

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: