Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 358,085 A→G 13.7% intergenic (+53/‑77) tdcG‑I → / → PP_5441 L‑serine dehydratase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947358,0850AG13.7% 130.7 / 12.5 59intergenic (+53/‑77)tdcG‑I/PP_5441L‑serine dehydratase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (24/26);  new base G (5/3);  total (30/29)
Fisher's exact test for biased strand distribution p-value = 7.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.16e-01

AATCCGCACCATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAGTGTGGGC  >  NC_002947/357939‑358229
                                                                                                                                                   |                                                                                                                                                
aaTCCGCACCATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGaa                                                                                                                                                 >  1:92100/1‑149 (MQ=255)
 aTCCGCACCATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACt                                                                                                                                               <  6:379602/149‑1 (MQ=255)
    cGCACCATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACtgtg                                                                                                                                            <  5:136097/149‑1 (MQ=255)
             cgcgATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACCGGAACTGTGCAATCCCTg                                                                                                                                   <  7:343611/149‑1 (MQ=255)
                gATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGGATTGCACAGTTCCTGTAg                                                                                                                                <  7:94828/149‑1 (MQ=37)
                  tACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGGATTGCACAGTTCCTGTAGGa                                                                                                                              <  8:308840/149‑1 (MQ=37)
                      ggCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCgg                                                                                                                          <  2:344267/149‑1 (MQ=255)
                      ggCGCCGATATCCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCgg                                                                                                                          <  4:38695/149‑1 (MQ=255)
                       gcgcCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGCACTGTGCAATCCCTGTAGGAGCcgg                                                                                                                         >  7:282938/1‑149 (MQ=255)
                       gcgcCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCggg                                                                                                                         <  1:194774/149‑1 (MQ=255)
                            gATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGc                                                                                                                    >  6:162058/1‑149 (MQ=255)
                                    cGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAag                                                                                                            >  8:260185/1‑149 (MQ=255)
                                      aCAAGTCCACAGAGCCCGCGCGGGGCGGCGTGGCGGTCTGCTCTAGCGCGGGCGGATTTTCTGCGCTGCCCTGCTGGGCCGCTTCGCGGGCAAGCCCGCTCCTCCCGGACCGGGGCACTCCGTGTCGGCGGGGGCTTGCCCGCGAagag                                                                                                          <  2:102704/149‑1 (MQ=255)
                                             cAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCAATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCGTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGTAGCGGGCTTGCCCGCGAAGAGGCCCtt                                                                                                    <  4:221955/148‑1 (MQ=255)
                                               aaGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAg                                                                                                 <  7:260185/149‑1 (MQ=255)
                                               aaGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACAATCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATACCTGGAGGCGCggg                                                                                                                         >  3:294492/1‑125 (MQ=38)
                                                agagACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCGGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAgg                                                                                                >  1:199913/1‑149 (MQ=255)
                                                agagACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCCGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAgg                                                                                                >  5:109929/1‑149 (MQ=255)
                                                 gagaCCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGc                                                                                               <  5:162058/149‑1 (MQ=255)
                                                      ccTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGAt                                                                                           >  8:172307/1‑149 (MQ=255)
                                                       cTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCtt                                                                                                                      <  4:204833/121‑1 (MQ=255)
                                                       cTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCtt                                                                                                                      >  3:204833/1‑121 (MQ=255)
                                                       cTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATg                                                                                          >  3:260484/1‑149 (MQ=255)
                                                       cTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCGACAGGAACTGTGGCATGCCGGTAGGAGCGGGCTTGGCCGCGAAGAGGGCCTTGAGGGTGATg                                                                                          >  1:274568/1‑149 (MQ=38)
                                                       cTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCCCTCCGACAGGACCTGTGGCATGCCGGGAGGGGCGGGGGTGGCCGGGAAGAGGCCCGTGAGGGTGATg                                                                                          >  8:155118/1‑149 (MQ=255)
                                                        tCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTTTGAGCTGACCTTCTGGGCCTCTTCGCGGGCCAGCCCCCTCCGACAGGAACGGTGCAATCCCTGGAGGAGCGGGCGTGCCCGCGAAGAGGCCCTTGAGGCTGGTg                                                                                          >  6:122443/1‑148 (MQ=37)
                                                          gcgcgGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGa                                                                                                              <  6:278018/126‑1 (MQ=255)
                                                          gcgcgGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGa                                                                                                              >  5:278018/1‑126 (MQ=255)
                                                                       gCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCaa                                                                                                                                         <  8:324733/86‑1 (MQ=255)
                                                                       gCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCaa                                                                                                                                         >  7:324733/1‑86 (MQ=255)
                                                                           tCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAACCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCt                                                                      <  2:92100/149‑1 (MQ=255)
                                                                                gcTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCGACAGGAACAGTGCAATACCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGGCACACCACCGCCTATggc                                                                 >  6:84266/1‑147 (MQ=37)
                                                                                    tCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCtt                                                             <  4:121294/149‑1 (MQ=255)
                                                                                      gAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTa                                                           <  2:199913/149‑1 (MQ=255)
                                                                                       aGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCATGAGGCTGATGCAGGTGCAACAAAACCGCCTATGCctgatt                                                            >  4:290535/1‑143 (MQ=255)
                                                                                        gTGCTCATTTACTCCGCTGCCCTTCTGGCCCTCTTCGCGGGCACGCCCGCGCCTACCGGAACTGTGCAACCCCTGTCGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATg                                                         <  7:155118/149‑1 (MQ=37)
                                                                                           cTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTg                                                      >  3:115629/1‑149 (MQ=255)
                                                                                            tGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGc                                                     >  8:346502/1‑149 (MQ=255)
                                                                                                      aGCTGACCTTCTGGGCCTCTTCGCGGGCAAACCCCCAACAAAAGGGAAAGTGGAAAACCAGTTGGAGCGGGCATGCCCCCGAAGAGGCCCATGAGGGGGGTTCAGGTGGAACACAACCCCCAAAGGCGCATTATTGTTCGGTgcttggg                                           >  4:73515/1‑143 (MQ=255)
                                                                                                          gACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGc                                                                                               >  5:117459/1‑93 (MQ=255)
                                                                                                                     ccTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGAt                            >  1:108744/1‑149 (MQ=255)
                                                                                                                       tctTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGGAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGAtcg                          >  1:251970/1‑149 (MQ=255)
                                                                                                                         ttCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGAtcgtc                        >  5:328506/1‑149 (MQ=255)
                                                                                                                           ccggggCAAGCCCGCTCCTACAGTAACTGTGCAATCCCTGTAGTAGGGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGAtcgtcg                       <  8:227650/145‑1 (MQ=255)
                                                                                                                            gcgGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGc                                                                        >  8:7953/1‑98 (MQ=255)
                                                                                                                            gcgGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGc                                                                        <  7:7953/98‑1 (MQ=255)
                                                                                                                            gcgGCCAAGCCCGCTCCTACAGCAACGGTGCAACCCCTGTAGCAGCGGGCTTCCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGa                     <  5:122443/149‑1 (MQ=255)
                                                                                                                             cgGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGAt                    <  4:115629/149‑1 (MQ=255)
                                                                                                                                    gCCCGCTCCTACAGGGATTGCACAGTTCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCCTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGa                     <  7:259521/141‑1 (MQ=25)
                                                                                                                                    gCCCGCTCCTACAGGGATTGCACAGTTCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCCTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGa                     >  8:259521/1‑141 (MQ=25)
                                                                                                                                     cccGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGAt            <  7:172307/149‑1 (MQ=255)
                                                                                                                                      ccGCTCCTACAGGGATTGCACAGTTCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGAt            >  7:236478/1‑148 (MQ=37)
                                                                                                                                       ctcTCCTACTGCAACTGGGCATTCCCTGTTGTAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCa          <  2:328306/147‑1 (MQ=255)
                                                                                                                                          tCCTACAGGGATTGCACAGTTCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAgtg       >  5:27424/1‑149 (MQ=37)
                                                                                                                                            cTACAGGGATTGCACAGTTCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGAAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATCTGCGCCCCGCATACCGATCGTCGGATGTCGCGATCAGtgtg     >  8:117048/1‑149 (MQ=25)
                                                                                                                                            cTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAgtgt      <  6:109929/148‑1 (MQ=255)
                                                                                                                                             tACCGGAACTGTGCAAGCCGTGTCGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAGTGTgg    <  3:73515/149‑1 (MQ=255)
                                                                                                                                              aCAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAGTGTggg   >  7:251144/1‑149 (MQ=255)
                                                                                                                                              aCAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGAGCAGTGTggg   >  7:236200/1‑149 (MQ=255)
                                                                                                                                               cAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAGTGTGGGc  <  8:251144/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
AATCCGCACCATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAGTGTGGGC  >  NC_002947/357939‑358229

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