Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,860,817 A→C 6.2% M78R (ATG→AGG)  ligA ← DNA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,860,8170AC 6.2% 202.1 / 4.7 65M78R (ATG→AGG) ligADNA ligase
Reads supporting (aligned to +/- strand):  ref base A (30/31);  new base C (2/2);  total (32/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 = 5.97e-01

GGTACAGCAGGCTCACGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAGGCGGTCGTATTCGG  >  NC_002947/4860677‑4860955
                                                                                                                                            |                                                                                                                                           
ggTACAGCAGGCTCACGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCATGCTCAGCCTGGGCaac                                                                                                                                     <  6:228589/149‑3 (MQ=255)
 gTACAGCAGGCTCACGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACtt                                                                                                                                    >  5:60961/1‑149 (MQ=255)
    cagcagGCTCACGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGt                                                                                                                                 <  1:14730/149‑1 (MQ=255)
      gcagGCTCACGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTgg                                                                                                                               >  2:4652/1‑149 (MQ=255)
            tCACGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTg                                                                                                                                >  1:154287/1‑142 (MQ=255)
            tCACGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTg                                                                                                                                <  2:154287/142‑1 (MQ=255)
            tCACGGCCAGGCCATCGAGCTTGGGTTCGCAACTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTAcc                                                                                                                         >  1:220862/1‑149 (MQ=255)
              aCGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCt                                                                                                                        <  2:100239/148‑1 (MQ=255)
               cGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTgg                                                                                                                      <  4:184584/149‑1 (MQ=255)
                ggccaggccaTCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGc                                                                                                                     <  6:295583/149‑1 (MQ=255)
                      gccaTCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAAcg                                                                                                               >  6:192701/1‑149 (MQ=255)
                         aTCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAAcgccg                                                                                                            <  4:133627/149‑1 (MQ=255)
                             aGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACgccgcc                                                                                                          >  7:338410/1‑147 (MQ=255)
                              gCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAg                                                                                                       >  5:233030/1‑149 (MQ=255)
                               cTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAgg                                                                                                      <  6:60961/149‑1 (MQ=255)
                                  gggTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGc                                                                                                     <  2:19015/147‑1 (MQ=255)
                                          aGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCa                                                                                                                                       <  5:160960/105‑1 (MQ=255)
                                          aGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCa                                                                                                                                       >  6:160960/1‑105 (MQ=255)
                                               tAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACg                                                                                     >  8:139160/1‑149 (MQ=255)
                                                  tCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCtg                                                                                  >  1:133274/1‑149 (MQ=255)
                                                  tCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCtg                                                                                  >  1:138564/1‑149 (MQ=255)
                                                    gACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGc                                                                                                     >  2:356415/1‑129 (MQ=255)
                                                    gACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGc                                                                                                     <  1:356415/129‑1 (MQ=255)
                                                    gACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCtgg                                                                                 <  5:192701/148‑1 (MQ=255)
                                                      ccGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCtgggtg                                                                              <  7:139160/149‑1 (MQ=255)
                                                           cccGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGa                                                                         <  1:358284/149‑1 (MQ=255)
                                                              ggCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTc                                                                      <  6:114079/149‑1 (MQ=255)
                                                                  gATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCgggg                                                                  >  2:282316/1‑149 (MQ=255)
                                                                   aTCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCggggg                                                                 >  4:82676/1‑149 (MQ=255)
                                                                   aTCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCgggg                                                                  <  4:118478/148‑1 (MQ=255)
                                                                      gAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCtgggtggg                                                                            <  2:309663/135‑1 (MQ=255)
                                                                      gAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACATCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCtgggtggg                                                                            >  1:309663/1‑135 (MQ=255)
                                                                      gAGGCCTTCTACCACGCGACGGCCCAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTa                                                              >  6:155411/1‑149 (MQ=255)
                                                                       aGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTAc                                                             >  3:287764/1‑149 (MQ=255)
                                                                           cTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAgg                                                         >  4:142028/1‑149 (MQ=255)
                                                                           cTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAgg                                                         >  2:175746/1‑149 (MQ=255)
                                                                                  cACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTccac                                                                                            <  6:344111/107‑1 (MQ=255)
                                                                                  cACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTccac                                                                                            >  5:344111/1‑107 (MQ=255)
                                                                                      cgACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTc                                              <  4:287764/149‑1 (MQ=255)
                                                                                       gACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCgg                                                    >  6:325185/1‑142 (MQ=255)
                                                                                       gACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCgg                                                    <  5:325185/142‑1 (MQ=255)
                                                                                          ggCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTg                                                                                                                   <  8:160573/77‑1 (MQ=255)
                                                                                          ggCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTg                                                                                                                   >  7:160573/1‑77 (MQ=255)
                                                                                              gAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTc                                        >  2:272669/1‑147 (MQ=255)
                                                                                              gAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCa                                      <  7:56951/149‑1 (MQ=255)
                                                                                               aaCTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAg                                     >  2:350242/1‑149 (MQ=255)
                                                                                               aaCTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAg                                     <  6:233030/149‑1 (MQ=255)
                                                                                                      gCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCtt                               >  1:169742/1‑148 (MQ=255)
                                                                                                             ggCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGt                       <  2:169742/149‑1 (MQ=255)
                                                                                                                  cttcGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGaa                   >  4:260022/1‑148 (MQ=255)
                                                                                                                  cttcGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGaa                   >  4:95041/1‑148 (MQ=255)
                                                                                                                  cttcGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGa                    <  1:350242/147‑1 (MQ=255)
                                                                                                                  cttcGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGa                    <  3:82676/147‑1 (MQ=255)
                                                                                                                   ttcGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGtt                      >  1:7403/1‑144 (MQ=255)
                                                                                                                   ttcGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGtt                      <  2:7403/144‑1 (MQ=255)
                                                                                                                     cGAAGGCGTTGCCCAGGCTGAGCATGGGCACTGCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCATCGCCTTGAGCTCGTTGAACAg                <  2:220862/148‑1 (MQ=255)
                                                                                                                      gAAGGCGTTGCCCAGGCTGAGCATGGGCGCTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAGGc              >  6:117713/1‑149 (MQ=255)
                                                                                                                       cacgaagtGCCCATGCTCAGCCTGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAGGCg             >  5:97540/8‑149 (MQ=255)
                                                                                                                       aaGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAGGCg             <  1:282316/149‑1 (MQ=255)
                                                                                                                          gCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGTGCTCGTTGAACAGGCGGTc          <  5:155411/149‑1 (MQ=255)
                                                                                                                             gtGCCCATGCTCAGCCTGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGt                       <  6:162539/132‑1 (MQ=255)
                                                                                                                             gtGCCCATGCTCAGCCTGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTCAGCTCGt                       >  5:162539/2‑133 (MQ=255)
                                                                                                                                 ccAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGACTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAGGCGGTCGTATTc    >  1:274366/1‑149 (MQ=255)
                                                                                                                                 ccAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAGGCGGTCGTATTCg   <  1:76487/149‑1 (MQ=255)
                                                                                                                                  cAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAGGCGGTCGTATTCgg  <  3:95041/149‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                           
GGTACAGCAGGCTCACGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAGGCGGTCGTATTCGG  >  NC_002947/4860677‑4860955

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: