Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,984,040 G→T 8.3% intergenic (+34/+156) ilvE → / ← PP_3512 branched‑chain‑amino acid aminotransferase/membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,984,0400GT 8.3% 182.0 / 5.9 60intergenic (+34/+156)ilvE/PP_3512branched‑chain‑amino acid aminotransferase/membrane protein
Reads supporting (aligned to +/- strand):  ref base G (32/23);  new base T (2/3);  total (34/26)
Fisher's exact test for biased strand distribution p-value = 6.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.46e-01

TTCCACGACCTGGAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGG  >  NC_002947/3983893‑3984186
                                                                                                                                                   |                                                                                                                                                  
ttCCACGACCTGGAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAg                                                                                                                                                    >  2:249968/1‑148 (MQ=255)
  ccACGACCTGGAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCaa                                                                                                                                                 >  5:234405/1‑149 (MQ=255)
    aCGACCTGGAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATg                                                                                                                                               <  1:249968/149‑1 (MQ=255)
    aCGACCTGGAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATg                                                                                                                                               <  4:276351/149‑1 (MQ=255)
      gACCTGGAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATgg                                                                                                                                              >  5:407783/1‑148 (MQ=255)
         cTGGAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATg                                                                                                                                          >  7:22307/1‑149 (MQ=255)
          tGGAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGc                                                                                                                                         >  1:416319/1‑149 (MQ=255)
           ggAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGcc                                                                                                                                        >  2:226602/1‑149 (MQ=255)
            ttccaGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGc                                                                                                                                         <  3:241243/143‑1 (MQ=255)
            gAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCATTGCTGGCATGcc                                                                                                                                        <  1:412531/148‑1 (MQ=255)
             aaaaGGTCGGCCCGGTTACCCATAAGCTCTACAACGATCTGACCGGTATCCAGAGAGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCATTGCTGGCATGcccc                                                                                                                                      <  4:258445/149‑1 (MQ=255)
                 ggTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGttt                                                                                                                                  <  4:60840/149‑1 (MQ=255)
                 ggTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCAGACCGCTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGACTGAAACGCATCAGATGAATGAAAACGGGGCATGCCAGCAAAGGCATGCACCGttt                                                                                                                                  >  4:90198/1‑149 (MQ=255)
                     ggCCCGGTTAGCCAGAAGGTCTACAACGAGGTGTCGTGTTTGCCGAGCGGCGATGTAGAGGCACCGGCAGGCTGGAGCGTCTAGGTCGCGTGATACGCATCAGCTGTATGAAAACGGGGCATGCCAGCTATGGCATGCCCCGtttttt                                                                                                                               <  3:90198/148‑1 (MQ=255)
                       cccGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGc                                                                                                                                             >  8:238493/1‑132 (MQ=255)
                       cccGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGc                                                                                                                                             <  7:238493/132‑1 (MQ=255)
                         cGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGGGTTGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGGATCAGCTGTATGAAAACGGGGAATGCCAGCAATGGCATGCCCCGTTTTTTCAtgt                                                                                                                          >  5:349721/1‑149 (MQ=255)
                               cccAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTac                                                                                                                    <  7:28961/149‑1 (MQ=255)
                                ccAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTaca                                                                                                                   <  1:178119/149‑1 (MQ=255)
                                   gAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACgg                                                                                                                >  6:34908/1‑149 (MQ=255)
                                           aCAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCg                                                                                                        >  7:348158/1‑149 (MQ=255)
                                            cAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTGATGTGTTTACACGGCCCCCTCGt                                                                                                       <  5:95985/149‑1 (MQ=255)
                                              aCGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAg                                                                                                     <  2:325193/149‑1 (MQ=255)
                                               cGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCACCGTTTTTTCATGTG‑TTACACGGCCCCCTCGTAg                                                                                                     >  4:397695/1‑147 (MQ=255)
                                                 aGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAg                                                                                                  <  2:416319/149‑1 (MQ=255)
                                                 aGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAg                                                                                                  >  5:3640/1‑149 (MQ=255)
                                                  gCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGc                                                                                                 <  8:22307/149‑1 (MQ=255)
                                                        cGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGtt                                                                                            >  7:235915/1‑148 (MQ=255)
                                                          gTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTc                                                                                                         >  4:148264/1‑133 (MQ=255)
                                                          gTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTc                                                                                                         <  3:148264/133‑1 (MQ=255)
                                                          gTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTAc                                                                                         >  6:306380/1‑149 (MQ=255)
                                                                agagCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGa                                                                                   <  6:407783/149‑1 (MQ=255)
                                                                      gcgATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGc                                                                             <  4:410682/149‑1 (MQ=255)
                                                                           gTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCacac                                                                        >  1:149308/1‑149 (MQ=255)
                                                                                 gCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCgc                                                                    >  7:45999/1‑147 (MQ=255)
                                                                                  cACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCgccag                                                                 >  5:384644/1‑149 (MQ=255)
                                                                                        cAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCaata                                                           <  6:234405/149‑1 (MQ=255)
                                                                                          ggCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCaataat                                                         <  1:226602/149‑1 (MQ=255)
                                                                                                   gTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGga                                                 <  5:306380/148‑1 (MQ=255)
                                                                                                      aaGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTG‑TTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCa                                            <  3:397695/149‑1 (MQ=255)
                                                                                                            gCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGtt                                       >  8:209624/1‑149 (MQ=255)
                                                                                                              cTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGga                                                 <  4:236628/137‑1 (MQ=255)
                                                                                                              cTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGga                                                 >  3:236628/1‑137 (MQ=255)
                                                                                                              cTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACACGCGAATAGGACACACCCGCCAGCCAATAATGCCTGGGAGATAACAGTTg                                      >  8:96627/1‑148 (MQ=255)
                                                                                                                 aaaCGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAg                                  <  6:3640/149‑1 (MQ=255)
                                                                                                                               tATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAg                    >  1:255891/1‑149 (MQ=255)
                                                                                                                                aTGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGc                   >  4:152046/1‑149 (MQ=255)
                                                                                                                                     aaCGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGc               >  5:316308/1‑148 (MQ=255)
                                                                                                                                      aCGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGc               >  5:255679/1‑147 (MQ=255)
                                                                                                                                        ggggCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCTAATAGGCCACACCCGCCCGCCAATAATGCCTGGGAGATCACCGTTGTCTGTGCAGCCGCATTTGCGGGCAAGc           <  2:149308/149‑1 (MQ=255)
                                                                                                                                         gggCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGcc          >  4:344585/1‑149 (MQ=255)
                                                                                                                                          ggCATGCCATTTCTGGCATGACCCGTTTTTTCATGTGTTTACACGGCCCCCTTGTTGGAGTGTGTTTACTCGCGAATAGGCTACACTCGCCAGCTAATAATGTCTGGGAGATCATTGTTGTTAGTGCa                              <  8:64955/128‑1 (MQ=255)
                                                                                                                                          ggCATGCCATTGCTGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGccc         >  7:5000/1‑149 (MQ=255)
                                                                                                                                          ggCATGCCATTGCTGGCATGCCCCGTTTTTTAATGTGTTTACAAGGCCCCCTCGTAGGAGCGGGTTTACCAGAGAATAGGCCACACCCGCAAGAAAAAAATGCCTGGGAGAACACCGTTGTCAGTgaa                              >  7:64955/1‑126 (MQ=255)
                                                                                                                                          ggCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGccc         >  4:394060/1‑149 (MQ=255)
                                                                                                                                           gCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCcg        >  5:177984/1‑149 (MQ=255)
                                                                                                                                                ccAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCg   >  7:55491/1‑149 (MQ=255)
                                                                                                                                                 cAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGTCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTACAGCCGCATTAGCGGGCAAGCCCGCGCCgg  <  8:235915/149‑1 (MQ=255)
                                                                                                                                                 cAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCgg  >  3:320511/1‑149 (MQ=255)
                                                                                                                                                 cAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCg   >  1:157713/1‑148 (MQ=255)
                                                                                                                                                 cAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCATGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGTCAAGCCCGCGCCgg  >  8:180872/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TTCCACGACCTGGAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGG  >  NC_002947/3983893‑3984186

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: ‑