Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,244,358 G→C 7.6% S221R (AGC→AGG)  aruH ← arginine‑‑pyruvate transaminase AruH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,244,3580GC 7.6% 205.8 / 6.4 66S221R (AGC→AGG) aruHarginine‑‑pyruvate transaminase AruH
Reads supporting (aligned to +/- strand):  ref base G (27/34);  new base C (2/3);  total (29/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.99e-01

GCCGTAGAGCATGCACAGGGCGAGGTTTTCCAGATGCGTGGCCAGTGCTGTCGAGCCGACCACCCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAGCGCCTCCCAGGTGGCGCGCGG‑CAGGCTCGCCCCTGAAGGGTTGTGCGGG  >  NC_002947/4244214‑4244502
                                                                                                                                                |                                                                                                                                                 
gCCGTAGAGCATGCACAGGGCGAGGTTTTCCAGATGCGTGGCCAGTGCTGTCGAGCCGACCACCCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGc                                                                                                                                               >  1:121595/1‑149 (MQ=255)
   gTAGAGCATGCACAGGGCGAGGTTTTCCAGATGCGTGGCCAGTGCTGTCGAGCCGACCACCCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCgg                                                                                                                                            >  3:102764/1‑149 (MQ=255)
     agagCATGCACAGGGCGAGGTTTTCCAGATGCGTGGCCAGTGCTGTCGAGCCGACCACCCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCgggg                                                                                                                                          >  5:14943/1‑149 (MQ=255)
         cATGCACAGGGCGAGGTTTTCCAGATGCGTGGCCAGTGCTGTCGAGCCGACCACCCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGa                                                                                                                                      <  5:2940/149‑1 (MQ=255)
         cATGCACAGGGCGAGGTTTTCCAGATGCGTGGCCAGTGCTGTCGAGCCGACCACCCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGCCTGCCgggcatgg                                                                                                                                      <  5:80022/149‑6 (MQ=255)
          aTGCACAGGGCGAGGTTTTCCAGATGCGTGGCCAGTGCTGTCGAGCCGACCACCCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGAc                                                                                                                                     <  2:89361/149‑1 (MQ=255)
             cacaGGGCGAGGTTTTCCAGATGCGTGGCCAGTGCTGTCGAGCCGACCACCCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCTGGCAGGCTGCCGGGGCTGACATg                                                                                                                                  <  7:208358/149‑1 (MQ=255)
              acaGGGCGAGGTTTTCCAGATGCGTGGCCAGTGCTGTCGAGCCGACCACCCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGCCTGCCgggc                                                                                                                                          >  4:53976/1‑139 (MQ=255)
              acaGGGCGAGGTTTTCCAGATGCGTGGCCAGTGCTGTCGAGCCGACCACCCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGCCTGCCgggc                                                                                                                                          <  3:53976/140‑2 (MQ=255)
              acaGGGCGAGGTTTTCCAGATGCGTGGCCAGTGCTGTCGAGCCGACCACCCACCCCATGCGCCAGCCGGTCATGCCATCCGATTTGGACAGGCTGTTCAACGTTGCAGTGAGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATg                                                                                                                                  >  2:309218/1‑148 (MQ=255)
                aGGGCGAGGTTTTCCAGATGCGTGGCCAGTGCTGTCGAGCCGACCACCCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTc                                                                                                                               >  7:12544/1‑149 (MQ=255)
                aGGGCGAGGTTTTCCAGATGCGTGGCCAGTGCTGTCGAGCCGACCACCCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTc                                                                                                                               <  6:14943/149‑1 (MQ=255)
                  ggCGAGGTTTTCCAGATGCGTGGCCAGTNCTGTCGAGCCGACCACCCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGc                                                                                                                             >  1:183891/1‑149 (MQ=255)
                  ggCGAGGTTTTCCAGATGCGTGGCCAGTGCTGTCGAGCCGACCACCCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGCCTGCCgggcat                                                                                                                                        >  6:184069/1‑135 (MQ=255)
                  ggCGAGGTTTTCCAGATGCGTGGCCAGTGCTGTCGAGCCGACCACCCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGCCTGCCgggcat                                                                                                                                        <  5:184069/138‑4 (MQ=255)
                             ccAGATGCGTGGCCAGTGCTGTCGAGCCGACCACCCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATAcag                                                                                                                   >  5:38807/1‑148 (MQ=255)
                              cAGATGCGTGGCCAGTGCTGTCGAGCCGACCACCCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATAcagca                                                                                                                 <  6:279681/149‑1 (MQ=255)
                                     gTGGCCAGTGCTGTCGAGCCGACCACCCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCt                                                                                                          >  1:193382/1‑149 (MQ=255)
                                         ccAGTGCTGTCGAGCCGACCACCCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGa                                                                                                                                      <  7:120976/117‑1 (MQ=255)
                                         ccAGTGCTGTCGAGCCGACCACCCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGa                                                                                                                                      >  8:120976/1‑117 (MQ=255)
                                                tGTCGAGCCGACCACCCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGt                                                                                               <  8:223115/149‑1 (MQ=255)
                                                 gTCGAGCCGACCACCCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGt                                                                                               >  1:135258/1‑148 (MQ=255)
                                                  tCGAGCCGACCACCCACCCCATGCGCCAGCCGGTCATTGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTc                                                                                              <  6:142622/148‑1 (MQ=255)
                                                  tCGAGCCGACCACCCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCa                                                                                             <  6:258116/149‑1 (MQ=255)
                                                  acGAGCCGACCACCCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCa                                                                                             <  6:225698/148‑1 (MQ=255)
                                                    gAGCCGACCACCCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCaga                                                                                           <  4:262290/149‑1 (MQ=255)
                                                       ccGACCACCCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGAt                                                                                        <  2:183891/149‑1 (MQ=255)
                                                         gACCACCCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGAt                                                                                        >  6:180282/1‑147 (MQ=255)
                                                          accaccCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAg                                                                                     >  1:137741/1‑149 (MQ=255)
                                                            cacccaccCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGc                                                                                    >  5:260695/1‑148 (MQ=255)
                                                             acccacccCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCca                                                                                  <  2:137741/149‑1 (MQ=255)
                                                               ccacccCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTc                                                                                                                               >  5:237203/1‑102 (MQ=255)
                                                               ccacccCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTc                                                                                                                               <  6:237203/102‑1 (MQ=255)
                                                               ccacccCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCca                                                                                  >  6:289520/1‑147 (MQ=255)
                                                               ccacccCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCca                                                                                  >  8:128102/1‑147 (MQ=255)
                                                                 acccCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGa                                                                              <  2:113007/149‑1 (MQ=255)
                                                                   cccATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTAATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGAt                                                                             >  6:260123/1‑148 (MQ=255)
                                                                    ccATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCa                                                                           >  7:182300/1‑149 (MQ=255)
                                                                     cATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATAca                                                                                                                    >  8:26838/1‑107 (MQ=255)
                                                                     cATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATAca                                                                                                                    <  7:26838/107‑1 (MQ=255)
                                                                      aTGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATg                                                                         <  2:193382/149‑1 (MQ=255)
                                                                       tGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATgg                                                                        <  2:135258/149‑1 (MQ=255)
                                                                            cAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATGGGCAAc                                                                   <  8:12544/149‑1 (MQ=255)
                                                                            cAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATGGGCAAc                                                                   <  4:102764/149‑1 (MQ=255)
                                                                                 ggTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATgg                                                                        <  1:288890/139‑1 (MQ=255)
                                                                                 ggTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATgg                                                                        >  2:288890/1‑139 (MQ=255)
                                                                                  gTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATGGGCAACGCACAg                                                             >  2:68251/1‑149 (MQ=255)
                                                                                           tGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGAGCTCGTCAGAGATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAg                                                    >  2:119379/1‑149 (MQ=255)
                                                                                             cGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGATATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAgc                                                   <  1:309218/148‑1 (MQ=255)
                                                                                             cGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAgcg                                                  <  6:38807/149‑1 (MQ=255)
                                                                                               aTTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAGCGcc                                                >  8:4321/1‑149 (MQ=255)
                                                                                                    gACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAGCGCCTcc                                             >  8:130270/1‑147 (MQ=255)
                                                                                                    gACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAGCGCCTCCCa                                           <  1:68251/149‑1 (MQ=255)
                                                                                                                aaCGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAGCGCCTCCCAGGTGGcgcgcg                                <  5:289520/148‑1 (MQ=255)
                                                                                                                aaCGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAGCGCCACCTGGGAGgcgctgg                                <  5:180282/148‑4 (MQ=255)
                                                                                                                aaCGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGACCATGGGCAACGCACAGTTCGGCCAGCGCCTCCCAGGTGGCGCGCgg                               >  6:194399/1‑149 (MQ=255)
                                                                                                                 tcGTTGCAGTGCGCTGGGCCATGCCCGGCAGGGTGCCGGGGCTGACATGCTCGCCTTCATTCAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGACACAGATCATGGGCAACGCACAGTTCGGCCAGCGCCTCCCAGGTg                                       <  2:71681/139‑1 (MQ=255)
                                                                                                                 aCGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAGCGCCTCCCAGGTGGCGCGCGG‑c                             <  7:4321/149‑1 (MQ=255)
                                                                                                                 aCGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCACCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAGCGCCTCCCAGGTggagatagg                               >  1:71681/1‑141 (MQ=255)
                                                                                                                   gTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAGCGCCTCCCAGGTGGCGCGCGG‑CAg                           <  7:130270/149‑1 (MQ=255)
                                                                                                                       cAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGACCATGGGCAACGCACAGTTCGGCCAGCGCCTCCCAGGTGGCGCGCGG‑CAGGCt                        <  5:194399/148‑1 (MQ=255)
                                                                                                                                 ggCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAGCGCCTCCCAGGTGGCGCGCGG‑CAGGCTCGCCCCTGaa              <  2:121595/148‑1 (MQ=255)
                                                                                                                                  gCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAGCGCCTCCCAGGTGGCGCGCGG‑CAGGCTCGCCCCTGAAgg            <  7:128102/149‑1 (MQ=255)
                                                                                                                                         ccGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAGCGCCTCCCAGGTGGCGCGCGG‑CAGGCTCGCCCCTGAAGGGTTGTGc     >  2:18885/1‑149 (MQ=255)
                                                                                                                                            gCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAGCGCCTCCCAGGTGGCGCGCGG‑CAGGCTCGCCCCTGAAGGGTTGTGCggg  <  6:260695/149‑1 (MQ=255)
                                                                                                                                             cAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAGCGCCTCCCAGGTGGCGCGCGGTCAGGCTCGCCCCTGAAGGGTTGTGCggg  <  1:18885/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
GCCGTAGAGCATGCACAGGGCGAGGTTTTCCAGATGCGTGGCCAGTGCTGTCGAGCCGACCACCCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAGCGCCTCCCAGGTGGCGCGCGG‑CAGGCTCGCCCCTGAAGGGTTGTGCGGG  >  NC_002947/4244214‑4244502

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