Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 111,820 G→A 6.5% D81N (GAC→AAC)  PP_0105 → cytochrome c oxidase assembly protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947111,8200GA 6.5% 179.0 / 2.6 61D81N (GAC→AAC) PP_0105cytochrome c oxidase assembly protein
Reads supporting (aligned to +/- strand):  ref base G (26/31);  new base A (2/2);  total (28/33)
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.37e-01

ACGTCATGTGCAAAGCCTTCGGTATCAACGGCAAGACCGGCGGCCAATACGAAGGTAGCCAGGTCAGTGACCCGACGCGTGCTGTGCGGGTGCAGTTCATGTCGACCAATGCCAGCGACATGGTCTGGGACTTCTACTCCACTGCCGACCAGCTGGAGGTCAACCCGGGGGCGGTGAACCAGATGATCTTCATTGCCCACAACCCGACCGACCGGCCGATGAGTGCGCAGGCAATACCCAGCA‑TTACCCCGGCCGAGGCCGCGGCCTACTTCCACAAGACCGAGTGCTTC  >  NC_002947/111674‑111963
                                                                                                                                                  |                                                                                                                                                
acgTCATGTGCAAAGCCTTCGGTATCAACGGCAAGACCGGCGGCCAATACGAAGGTAGCCAGGTCAGTGACCCGACGCGTGCTGTGCGGGTGCAGTTCATGTCGACCAATGCCAGCGACATGGTCTGGGACTTCTACTCCACTGCCGAc                                                                                                                                                <  7:121590/149‑1 (MQ=255)
   tCATGTGCAAAGCCTTCGGTATCAACGGCAAGACCGGCGGCCAATACGAAGGTAGCCAGGTCAGTGACCCGACGCGTGCTGTGCGGGTGCAGTTCATGTCGACCAATGCCAGCGACATGGTCTGGGACTTCTACTCCACTGCCGACCAg                                                                                                                                             >  1:356066/1‑149 (MQ=255)
        tgCAAAGCCTTCGGTATCAACGGCAAGACCGGCGGCCAATACGAAGGTAGCCAGGTCAGTGACCCGACGCGTGCTGTGCGGGTGCAGTTCATGTCGACCAATGCCAGCGACATGGTCTGGGACTTCTACTCCACTGCCGACCAGCTGGa                                                                                                                                        <  2:113683/149‑1 (MQ=255)
           aaaGCCTTCGGTATCAACGGCAAGACCGGCGGCCAATACGAAGGTAGCCAGGTCAGTGACCCGACGCGTGCTGTGCGGGTGCAGTTCATGTCGACCAATGCCAGCGACATGGTCTGGGACTTCTACTCCACTGCCGACCAGCTGGAGGt                                                                                                                                     <  4:58916/149‑1 (MQ=255)
                attCGGTATCAACGGCAAGACCGGCGGCCAATACGAAGGTAGCCATGTCAGTGACCCGACGCGTGCTGTGCGGGTGCAGTTCATGTCGACCAATGCCAGCGACATGCTCTGGGACTTCTACTCCACTGCCGACCAGCTGGAGGTCAAcc                                                                                                                                >  6:154776/2‑149 (MQ=255)
                 ttCGGTATCAACGGCAAGACCGGCGGCCAATACGAAGGTAGCCAGGTCAGTGACCCGACGCGTGCTGTGCGGGTGCAGTTCATGTCGACCAATGCCAGCGACATGGTCTGGGACTTCTACTCCACTGCCGACCAGCTGGAGGTCAAccc                                                                                                                               >  3:183805/1‑149 (MQ=255)
                     gTATCAACGGCAAGACCGGCGGCCAATACGAAGGTAGCCAGGTCAGTGACCCGACGCGTGCTGTGCGGGTGCAGTTCATGTCGACCAATGCCAGCGACATGGTCTGGGACTTCTACTCCACTGCCGACCAGCTGGAGGTCAAcccttgg                                                                                                                           <  6:347620/149‑5 (MQ=255)
                     gTATCAACGGCAAGACCGGCGGCCAATACGAAGGTAGCCAGGTCAGTGACCCGACGCGTGCTGTGCGGGTGCAGTTCATGTCGACCAATGCCAGCGACATGGTCTGGGACTTCTACTCCACTGCCGACCAGCTGGAGGTCAACCCgggg                                                                                                                           <  6:26780/149‑1 (MQ=255)
                      tATCAACGGCTAGACCGGAGGCCAATACGAAGGTAGCCAGGTCAGTGACCGGACTCGTGCTGTGCGGGTGCAGTTGATGTCGACCTATGCGAGCGACATGGTCTGGGACTTCTACTCCACTGTCGACCAGCTGGAGGt                                                                                                                                     <  4:63715/138‑1 (MQ=255)
                      tATCAACGGCAAGACCGGCGGCCAATACGAAGCTAGCCAGGTCAGTGACCCGACGCGTGCTGTGCGGGTGCAGTTCATGTCGACCAATGACAGCGACATGGTCTGGGACTTCTACTCCACTGTCGACCACCTGGAGGt                                                                                                                                     >  3:63715/1‑138 (MQ=255)
                        tCAACGGCAAGACCGGCGGCCAATACGAAGGTAGCCAGGTCAGTGACCCGACGCGTGCTGTGCGGGTGCAGTTCATGTCGACCAATGCCAGCGACATGGTCTGGGACTTCTACTCCACTGCCGACCAGCTGGAGGTc                                                                                                                                    <  1:239304/137‑1 (MQ=255)
                        tCAACGGCAAGACCGGCGGCCAATACGAAGGTAGCCAGGTCAGTGACCCGACGCGTGCTGTGCGGGTGCAGTTCATGTCGACCAATGCCAGCGACATGGTCTGGGACTTCTACTCCACTGCCGACCAGCTGGAGGTc                                                                                                                                    >  2:239304/1‑137 (MQ=255)
                        tCAACGGCAAGACCGGCGGCCAATACGAAGGTAGCCAGGTCAGTGACCCGACGCGTGCTGTGCGGGTGCAGTTCATGTCGACCAATGCCAGCGACATGGTCTGGGACTTCTACTCCACTGCCGACCAGCTGGAGGTCAACCCGGGGGCg                                                                                                                        >  8:245984/1‑149 (MQ=255)
                         cAACGGCAAGACCGGCGGCCAATACGAAGGTAGCCAGGTCAGTGACCCGACGCGTGCTGTGCGGGTGCAGTTCATGTCGACCAATGCCAGCGACATGGTCTGGGACTTCTACTCCACTGCCGACCAGCTGGAGGTCAACCCgctgggg                                                                                                                        >  1:131509/1‑142 (MQ=255)
                          aaCGGAAAGACCGGCGGCCAGTACGAAGGTAGCCAGGTCAGTGACCCGACGCGTGCTGTGCGGGTGCAGTTCATGTCGACCAATGCCAGCGACATGGTCTGGGACTTCTACTCCACTGCCGACCAGCTGGAGGTCAACCCGGGGGCGGt                                                                                                                      >  4:226069/1‑149 (MQ=255)
                           aCGGCAAGACCGGCGGCCAATACGAAGGTAGCCAGGTCAGTGACCCGACGCGTGCTGTGCGGGTGCAGTTCATGTCGACCAATGCCAGCGACATGGTCTGGGACTTCTACTCCACTGCCGACCAGCTGGAGGTCAACCCGGGGGCGGt                                                                                                                      >  5:165901/1‑148 (MQ=255)
                           aCGGCAAGACCGGCGGCCAATACGAAGGTAGCCAGGTCAGTGACCCGACGCGTGCTGTGCGGGTGCAGTTCATGTCGACCAATGCCAGCGACATGGTCTGGGACTTCTACTCCACTGCCGACCAGCTGGAGGTCAACCCGGGGGCGGt                                                                                                                      >  4:349782/1‑148 (MQ=255)
                           aCGGCAAGACCGGCGGCCAATACGAAGGTAGCCAGGTCAGTGACCCGACGCGTGCTGTGCGGGTGCAGTTCATGTCGACCAATGCCAGCGACATGGTCTGGGACTTCTACTCCACTGCCGACCAGCTGGAGGTCAACCCGGGGGCGGt                                                                                                                      >  4:151516/1‑148 (MQ=255)
                           aCGGCAAGACCGGCGGCCAATACGAAGGTAGCCAGGTCAGTGACCCGACGCGTGCTGTGCGGGTGCAGTTCATGTCGACCAATGCCAGCGACATGGTCTGGGACTTCTACTCCACTGCCGACCAGCTGGAGGTCAACCCGGGGGCGGt                                                                                                                      >  3:97194/1‑148 (MQ=255)
                             ggCAAGACCGGCGGCCAATACGAAGGTAGCCAGGTCAGTGACCCGACGCGTGCTGTGCGGGTGCAGTTCATGTCGACCAATGCCAGCGACATGGTCTGGGACTTCTACTCCACTGCCGACCAGCTGGAGGTCAACCCGGGGGCGGTGaa                                                                                                                   <  1:23564/149‑1 (MQ=255)
                              gCAAGACCGGCGGCCAGTACGAAGGTAGCCAGGTCAGTGACCCGACGCGTGCTGTGCGGGTGCAGTTCATGTCGACCAATGCCAGCGACATGGTCTGGGACTTCTACTCCACTGCCGACCAGCTGGAGGTCAACCCGGGGGCGGTGAAc                                                                                                                  <  3:226069/149‑1 (MQ=255)
                              gCAAGACCGGCGGCCAATACGAAGGTAGCCAGGTCAGTGACCCGACGCGTGCTGTGCGGGTGCAGTTCATGTCGACCAATGCCAGCGACATGGTCTGGGACTTCTACTCCACTGCCAACCAGCTGGAGGTCAACCCGGGGGCGGTGAAc                                                                                                                  >  2:357311/1‑149 (MQ=255)
                               cAAGACCGGCGGCCAATACGAAGGTAGCCAGGTCAGTGACCCGACGCGTGCTGTGCGGGTGCAGTTCATGTCGACCAATGCCAGCGACATGGTCTGGGACTTCTACTCCACTGCCGACCAGCTGGAGGTCAACCCGGGGGCGGTGAAc                                                                                                                  >  1:107949/1‑148 (MQ=255)
                                  gACCGGCGGCCAATACGAAGGTAGCCAGGTCAGTGACCCGACGCGTGCTGTGCGGGTGCAGTTCATGTCGACCAATGCCAGCGACATGGTCTGGGACTTCTACTCCACTGCCGACCAGCTGGAGGTCAACCCGGGGGCGGTGAAc                                                                                                                  >  7:161493/1‑145 (MQ=255)
                                  gACCGGCGGCCAATACGAAGGTAGCCAGGTCAGTGACCCGACGCGTGCTGTGCGGGTGCAGTTCATGTCGACCAATGCCAGCGACATGGTCTGGGACTTCTACTCCACTGCCGACCAGCTGGAGGTCAACCCGGGGGCGGTGAAc                                                                                                                  >  2:199529/1‑145 (MQ=255)
                                  gACCGGCGGCCAATACGAAGGTAGCCAGGTCAGTGACCCGACGCGTGCTGTGCGGGTGCAGTTCATGTCGACCAATGCCAGCGACATGGTCTGGGACTTCTACTCCACTGCCGACCAGCTGGAGGTCAACCCGGGGGCGGTGAACCAga                                                                                                              <  7:245984/149‑1 (MQ=255)
                                  gACCGGCGGCCAATACGAAGGTAGCCAGGTCAGTGACCCGACGCGTGCTGTGCGGGTGCAGTTCATGTCGACCAATGCCAGCGACATGGTCTGGGACTTCTACTCCACTGCCGACCAGCTGGAGGTCAACCCGGGGGCGGTGAACCAga                                                                                                              <  2:107949/149‑1 (MQ=255)
                                     cggcggCCAATACGAAGGTAGCCAGGTCAGTGACCCGACGCGTGCTGTGCGGGTGCAGTTCATGTCGACCAATGCCAGCGACATGGTCTGGGACTTCTACTCCACTGCCGCCCAGCTGGAGGTCAACCCGGGGGCGGTGAACCAgatga                                                                                                           >  4:326172/1‑149 (MQ=255)
                                           ccAATACGAAGGTAGCCAGGTCAGTGACCCGACGCGTGCTGTGCGGGTGCAGTTCATGTCGACCAATGCCAGCGACATGGTCTGGGACTTCTACTCCACTGCCGACCAGCTGGAGGTCAACCCGGGGGCGGTGAACCAGATGATCTTCa                                                                                                     <  2:131509/149‑1 (MQ=255)
                                             aaTACGAAGGTAGCCAGGTCAGTGACCCGACGCGTGCTGTGCGGGTGCAGTTCATGTCGACCAATGCCAGCGACATGGTCTGGGACTTCTACTCCACTGCCGACCAGCTGGAGGTCAACCCGGGGGCGGTGAACCAGATGATCTTCAtt                                                                                                   <  4:183805/149‑1 (MQ=255)
                                                 cGAAGGTAGCCAGGTCAGTGACCCGACGCGTGCTGTGCGGGTGCAGTTCATGTCGACCAATGCCAGCGACATGGTCTGGGACTTCTACTCCACTGCCGACCAGCTGGAGGTCAACCCGGGGGCGGTGAACCAGATGATCTTCATTGcc                                                                                                <  8:161493/148‑1 (MQ=255)
                                                  gAAGGTAGCCAGGTCAGTGACCCGACGCGTGCTGTGCGGGTGCAGTTCATGTCGACCAATGCCAGCGACATGGTCTGGGACTTCTACTCCACTGCCGACCAGCTGGAGGTCAACCCGGGGGCGGTGAACCAGATGATCTTCATTGCCca                                                                                              <  6:165901/149‑1 (MQ=255)
                                                                cAGTGACCCGACGCGTGCTGTGCGGGTGCAGTTCATGTCGACCAATGCCAGCGACATGGTCTGGGACTTCTACTCCACTGCCGACCAGCTGGAGGTCAACCCGGGGGCGGTGAACCAGATGATCTTCATTGCCCACAACCCgaccgacc                                                                                <  3:157325/149‑1 (MQ=255)
                                                                 aGTGACCCGACGCGTGCTGTGCGGGTGCAGTTCATGTCGACCAATGCCAGCGACATGGTCTGGGACTTCTACTCCACTGCCGACCAGCTGGAGGTCAACCCGGGGGCGGTGAACCAGATGATCTTCATTGCCCACAACCCgaccgacc                                                                                <  3:163055/148‑1 (MQ=255)
                                                                       ccGACGCGTGCTGTGCGGGTGCAGTTCATGTCGACCAATGCCAGCGACATGGTCTGGGACTTCTACTCCACTGCCGACCAGCTGGAGGTCAACCCGGGTGCGGTGAACCAGATGATCTTCATTGCCCACAACCCGACCGACCGGCCGAt                                                                         >  8:134396/1‑149 (MQ=255)
                                                                       ccGACGCGTGCTGTGCGGGTGCAGTTCATGTCGACCAATGCCAGCGACATGGTCTGGGACTTCTACTCCACTGCCGACCAGCTGGAGGTCAACCCGGGGGCGGTGAACCAGATGATCTTCATTGCCCACAACCCGACCGACCGGCCGAt                                                                         <  3:349782/149‑1 (MQ=255)
                                                                       ccGACGCGTGCTGTGCGGGTGCAGTTCATGTCGACCAATGCCAGCGACATGGTCTGGGACTTCTACTCCACTGCCGACCAGCTGGAGGTCAACCCGGGGGCGGTGAACCAGATGATCTTCATTGCCCACAACCCGACCGACCGGCCGAt                                                                         <  4:97194/149‑1 (MQ=255)
                                                                       ccGACGCGTGCTGTGCGGGTGCAGTTCATGTCGACCAATGCCAGCGACATGGTCTGGGACTTCTACTCCACTGCCGACCAGCTGGAGGTCAACCCGGGGGCGGTGAACCAGATGATCTTCATTGCCCACAACCCGACCGACCGGCCGAt                                                                         >  6:380822/1‑149 (MQ=255)
                                                                       ccGACGCGTGCTGTGCGGGTGCAGTTCATGTCGACCAATGCCAGCGACATGGTCTGGGACTTCTACTCCACTGCCAACCAGCTGGAGGTCAACCCGGGGGCGGTGAACCAGATGATCTTCATTGCCCACAACCCGACCGACCGGCCGAt                                                                         <  1:357311/149‑1 (MQ=255)
                                                                                  tgtgCGGGTGCAGTTCATGTCGACCAATGCCAGCGACATGGTCTGGGACTTCTACTCCACTGCCGACCAGCTGGAGGTCAACCCGGGGGCGGCGAACCAGATGATCTTCATTGCCCACAACCCGACCGACCGGCCGATGAGTGCGCAgg                                                              >  7:205868/1‑149 (MQ=255)
                                                                                      cGGGTGCAGTTCATGTCGACCAATGCCAGCGACATGGTCTGGGACTTCTACTCCACTGCCGACCAGCTGGAGGTCAACCCGGGGGCGGTGAACCAGATGATCTTCATTGCCCACAACCCGACCGACCGGCCGATGAGTGCGCAGGCAAt                                                          <  1:199529/149‑1 (MQ=255)
                                                                                        ggTGCAGTTCATGTCGACCAATGCCAGCGACATGGTCTGGGACTTCTACTCCACTGCCGACCAGCTGGAGGTCAACCCGGGGGCGGTGAACCAGATGATCTTCATTGCCCACAACCCGACCGACCGGCCGATGAGTGCGCAGGCAATAc                                                        >  2:258187/1‑149 (MQ=255)
                                                                                          tGCAGTTCATGTCGACCAATGCCAGCGACATGGTCTGGGACTTCTACTCCACTGCCGACCAGCTGGAGGTCAACCCGGGGGCGGTGAACCAGATGATCTTCATTGCCCACAACCCGACCGACCGGCCGAt                                                                         >  7:292778/1‑130 (MQ=255)
                                                                                          tGCAGTTCATGTCGACCAATGCCAGCGACATGGTCTGGGACTTCTACTCCACTGCCGACCAGCTGGAGGTCAACCCGGGGGCGGTGAACCAGATGATCTTCATTGCCCACAACCCGACCGACCGGCCGAt                                                                         <  8:292778/130‑1 (MQ=255)
                                                                                              gTTCATGTCGACCAATGCCAGCGACATGGTCTGGGACTTCTACTCCACTGCCGACCAGCTGGAGGTCAACCCGGGGGCGGTGAACCAGATGATCTTCATTGCCCACAACCCGACCGACCGGCCGATGAGTGCGCAGGCAATACCCAGCa                                                  <  1:42422/149‑1 (MQ=255)
                                                                                                   tGTCGACCAATGCCAGCGACATGGTCTGGGACTTCTACTCCACTGCCGACCAGCTGGAGGTCAACCCGGGGGCGGTGAACCAGATGATCTTCATTGCCCACAACCCGACCGACCGGCCGATGAGTGCGCAGGCAATACCCAGCA‑TTAcc                                            <  5:380822/149‑1 (MQ=255)
                                                                                                            aTGCCAGCGACATGGTCTGGGACTTCTACTCCACTGCCGACCAGCTGGAGGTCAACCCGGGGGCGGCGAACCAGATGATCTTCATTGCCCACAACCCGACCGACCGGCCGATGAGTGCGCAGGCAATACCCAGCA‑TTACCCCGGCCGAg                                   <  8:205868/149‑1 (MQ=255)
                                                                                                                          gTCTGTGACTTCTACTCCACTGCCAACCAGCTGGAGGTCAACCCGGGGGCGGTGAACCAGATGATCTTCATTGCCCACAACCCGACCGACCGGCCGATGAGTGCGCAGGCAATACCCAGCA‑TTACCCCGGCCGAGGCCGCGGc                           <  5:97068/143‑1 (MQ=255)
                                                                                                                          gTCTGGGACTTCTACTCCACTGCCAACCAGCTGGAGGTCAACCCGGGGGCGGTGAACCAGATGATCTTCATTGCCCACAACCCGACCGACCGGCCGATGAGTGCGCAGGCAATACCCAGCA‑TTACCCCGGCCGAGGCCGCGGc                           >  6:97068/1‑143 (MQ=255)
                                                                                                                           tCTGGGACTTCTACTCCACTGCCGACCAGCTGGAGGTCAACCCGGGGGCGGTGAACCAGATGATCTTCATTGCCCACAACCCGACCGACCGGCCGATGAGTGCGCAGGCAATACCCAGCA‑TTACCCCGGCCGAGGCCGCGGCCTACTTc                    <  3:326172/149‑1 (MQ=255)
                                                                                                                               ggACTTCTACTCCACTGCCGACCAGCTGGAGGTCAACCCGGGGGCGGTGAACCAGATGATCTTCATTGCCCACAACCCGACCGACCGGCCGATGAGTGCGCAGGCAATACCCAGCA‑TTACCCCGGCCGAGGCCGCGGCCTACTTCcac                 <  7:74120/148‑1 (MQ=255)
                                                                                                                                gACTTCTACTCCACTGCCGACCAGCTGGAGGTCAACCCGGGGGCGGTGAACCAGATGATCTTCATTGCCCACAACCCGACCGACCGGCCGATGAGTGCGCAGGCAATACCCAGCA‑TTACCCCGGCCTAGGCCGCGGCCTACTTCCACaa               <  5:154776/149‑1 (MQ=255)
                                                                                                                                gACTTCTACTCCACTGCCGACCAGCTGGAGGTCAACCCGGGGGCGGTGAACCAGATGATCTTCATTGCCCACAACCCGACCGACCGGCCGATGAGTGCGCAGGCAATACCCAGCA‑TTACCCCGGCCGAGGCCGCGGCCTACTTCcac                 >  2:316406/1‑147 (MQ=255)
                                                                                                                                   ttCTACTCCACTGCCGACCAGCTGGAGGTCAACCCGGGGGCGGTGAACCAGATGATCTTCATTGCCCACAACCCGACCGACCGGCCGATGAGTGCGCAg                                                               >  2:100804/1‑99 (MQ=255)
                                                                                                                                   ttCTACTCCACTGCCGACCAGCTGGAGGTCAACCCGGGGGCGGTGAACCAGATGATCTTCATTGCCCACAACCCGACCGACCGGCCGATGAGTGCGCAg                                                               <  1:100804/99‑1 (MQ=255)
                                                                                                                                      tACTCCACTGCCGACCAGCTGGAGGTCAACCCGGGGGCGGTGAACCAGATGATCTTCATTGCCCACAACCCGACCGACCGGCCGATGAGTGCGCAGGCAATACCCAGCA‑TTACCCCGGCCGAGGCCGCGGCCTACTTCCACAAGAcc           >  8:13825/1‑147 (MQ=255)
                                                                                                                                        cTCCACTGCCGACCAGCTGGAGGTCAACCCGGGGGCGGTGAACCAGATGATCTTCATTGCCCACAACCCGACCGACCGGCCGATGAGTGCGCAGGCAATACCCAGCA‑TTACCCCGGCCGAGGCCGCGGCCTACTTCCACAAGACCGa         <  3:151516/147‑1 (MQ=255)
                                                                                                                                         tCCACTGCCGACCAGCTGGAGGTCAACCCGGGGGCGGTGAACCAGATGATCTTCATTGCCCACAACCCGACCGACCGGCCGATGAGTGCGCAGGCAATACCCAGCACTTACCCCGGCCGAGGCCGCGGCCTACTTCCACAAGACCGAGt       <  1:258187/149‑1 (MQ=255)
                                                                                                                                           cACTGCCGACCAGCTGGAGGTCAACCCGGGGGCGGTGAACCAGATGATCTTCATTGCCCACAACCCGACCGACCGGCCGATGAGTGCGCAGGCAATACCCAGCA‑TTACCCCGGCCGAGGCCGCGGCCTACTTCCACAAGACCGAGTGCt    >  3:368064/1‑149 (MQ=255)
                                                                                                                                           cACTGCCGACCAGCTGGAGGTCAACCCGGGGGCGGTGAACCAGATGATCTTCATTGCCCACAACCCGACCGACCGGCCGATGAGTGCGCAGGCAATACCCAGCA‑TTACCCCGGCCGAGGCCGCGGCCTACTTCCACAAGACCGAGTGCt    >  6:105456/1‑149 (MQ=255)
                                                                                                                                               gCCGACCAGCTGGAGGTCAACCCGGGGGCGGTGAACCAGATGATCTTCATTGCCCACAACCCGACCGACCGGCCGATGAGTGCGCAGGCAATACCCAGCA‑TTACCCCGGCCGAGGCCGCGGCCTACTTCCACAAGACCGAGTGCTTc  <  1:316406/147‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
ACGTCATGTGCAAAGCCTTCGGTATCAACGGCAAGACCGGCGGCCAATACGAAGGTAGCCAGGTCAGTGACCCGACGCGTGCTGTGCGGGTGCAGTTCATGTCGACCAATGCCAGCGACATGGTCTGGGACTTCTACTCCACTGCCGACCAGCTGGAGGTCAACCCGGGGGCGGTGAACCAGATGATCTTCATTGCCCACAACCCGACCGACCGGCCGATGAGTGCGCAGGCAATACCCAGCA‑TTACCCCGGCCGAGGCCGCGGCCTACTTCCACAAGACCGAGTGCTTC  >  NC_002947/111674‑111963

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