Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,555,065 A→T 11.1% E379D (GAA→GAT) ‡ sbcB → exodeoxyribonuclease I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,555,0650AT11.1% 159.7 / 11.2 63E379D (GAA→GAT) ‡sbcBexodeoxyribonuclease I
Reads supporting (aligned to +/- strand):  ref base A (26/30);  new base T (3/4);  total (29/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.23e-01

CCAATCAACAGGCGCAATGGCAAGAAAAGCTGGAGCACATCTATGGCAAGGACGAGTTCGCCCCGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTGGTATGGCTTCTGCCAACAGCGCCTGA  >  NC_002947/1554917‑1555210
                                                                                                                                                    |                                                                                                                                                 
ccAATCAACAGGCGCAATGGCAAGAAAAGCTGGAGCACATCTATGGCAAGGACGAGTTCGCCCCGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGaa                                                                                                                                                   <  6:207532/149‑1 (MQ=255)
  aaTCAACAGGCGCAATGGCAAGAAAAGCTGGAGCACATCTATGGCAAGGACGAGTTGGCCCCGAGCGAGGACCCGGAACAGCAGTTGTATGACGGTTTTATAGGAGACCGCGACCGGCGCTTAAGCGAGCAAGTTCGTGCGCTGGAAcc                                                                                                                                                 >  1:255970/1‑149 (MQ=255)
      aaaaGGAGCAATGGCAAGAAAAGCTGGAGAACATCTATGGCAAGGACGAGTTAGCCCATAGCGAGGACCCGGAACAGCAGTTTTATGAGGGAGTTATAGGTGACCGAGACCGTCTCTAATTCGATCAAATTCGTGCGCTAGAACCGgct                                                                                                                                             >  8:69621/4‑148 (MQ=255)
       aCAGGCGCAATGGCAAGAAAAGCTGGAGCACATCTATGGCAAGGACGAGTTCGCCCCGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGgcg                                                                                                                                             >  5:287076/1‑148 (MQ=255)
        cAGGCGCAATGGCAAGAAAAGCTGGAGCACATCTATGGCAAGGACGAGTTCGCCCCGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCCGGTTCCAGCGCa                                                                                                                                           <  7:286359/149‑1 (MQ=255)
        cAGGCGCAATGGCAAGAAAAGCTGGAGCACATCTATGGCAAGGACGAGTTCGCCCCGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCCGGTTCCAGCGCa                                                                                                                                           <  6:35522/149‑1 (MQ=255)
         aGGCGCAATGGCAAGAAAAGCTGGAGCACATCTATGGCAAGGACGAGTTCGCCCCGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCATGCGCCGGTTCCAGCGCac                                                                                                                                          <  1:213867/149‑2 (MQ=255)
           gcgcAATGGCAAGAAAAGCTGGAGCACATCTATGGCAAGGACGAGTTCGCCCCGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGtt                                                                                                                                        <  5:112207/149‑1 (MQ=255)
             gcAATGGCAAGAAAAGCTGGAGCACATCTATGGCAAGGACGAGTTCGCCCCGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGa                                                                                                                                      <  5:210264/149‑1 (MQ=255)
               aaTGGCAAGAAAAGCTGGAGCACATCTATGGCAAGGACGAGTTCGCCCCGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGc                                                                                                                                    <  7:50695/149‑1 (MQ=255)
                   gCAAGAAAAGCTGGAGCACATCTATGGCAAGGACGAGTTCGCCCCGAGCGAGGACCCGGAACCGCAGTTGTATGAGGGTTTTATAGGCGACCGCGCCCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATg                                                                                                                                <  4:21074/149‑1 (MQ=255)
                    cAAGAAAAGCTGGAGCACATCTATGGCAAGGACGAGTTCGCCCCGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATgg                                                                                                                               >  7:173899/1‑149 (MQ=255)
                     aaGAAAAGCTGGAGCACATCTATGGCAAGGACGAGTTCGCCCCGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGc                                                                                                                              >  3:118730/1‑149 (MQ=255)
                     aaGAAAAGCTGGAGCACATCTATGGCAAGGACGAGTTCGCCCCGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGc                                                                                                                              <  4:38581/149‑1 (MQ=255)
                            gCTGGAGCACATCTATGGCAAGGACGAGTTCGCCCCGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGa                                                                                                                       <  7:304217/149‑1 (MQ=255)
                                      aTCTATGGCAAGGACGAGTTCGCCCCGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCgatga                                                                                                             >  5:306486/1‑149 (MQ=255)
                                        cTATGGCAAGGACGAGTTCGCCACGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACAGTCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGCCCAGAGGATGtt                                                                                                                   >  2:236521/1‑141 (MQ=255)
                                        cTATGGCAAGGACGAGATCGCCCCGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCATGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGtt                                                                                                                   <  1:236521/141‑1 (MQ=255)
                                            ggCAAGGACGAGTTCGCCCCGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCcgcg                                                                                                       >  3:187998/1‑149 (MQ=255)
                                               aaGGACGAGTTCGCCCCGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGttt                                                                                                    <  6:287076/149‑1 (MQ=255)
                                                 ggACGAGTTCGCCCCGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGc                                                                                                  >  3:244738/1‑149 (MQ=255)
                                                 ggACGAGTTCGCCCCGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGc                                                                                                  >  1:34092/1‑149 (MQ=255)
                                                   aCGAGTTCGCCCCGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCg                                                                                                <  6:306486/149‑1 (MQ=255)
                                                      aGTTCGCCCCGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGa                                                                                                                                      <  6:125710/108‑1 (MQ=255)
                                                      aGTTCGCCCCGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGa                                                                                                                                      >  5:125710/1‑108 (MQ=255)
                                                           gCCCCGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCtgtt                                                                                        >  4:301517/1‑149 (MQ=255)
                                                              ccGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTGTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTgttgtt                                                                                     <  2:255970/149‑1 (MQ=255)
                                                               cgagcgagGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTc                                                                                    <  4:138027/149‑1 (MQ=255)
                                                                          ccGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCggg                                                                          >  1:264870/1‑148 (MQ=255)
                                                                              aaCAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAAc                                                                                                                                                  <  4:153520/72‑1 (MQ=255)
                                                                              aaCAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAAc                                                                                                                                                  >  3:153520/1‑72 (MQ=255)
                                                                                 agcagTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGATGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTacc                                                                  >  8:71932/1‑147 (MQ=255)
                                                                                     gTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTcc                                                              <  4:89648/149‑1 (MQ=255)
                                                                                        gTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCt                                                            >  3:18680/1‑148 (MQ=255)
                                                                                         tATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTg                                                           >  7:27238/1‑148 (MQ=255)
                                                                                         tATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGa                                                          <  2:150269/149‑1 (MQ=255)
                                                                                          aTGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGCGCCATGGCCAGTGGATGTTCGATGATCCGCGTTGGCCGGAGATGTTGTTCCGCTACCGGGCCGGTAGCTTccaggaa                                                         >  8:54898/1‑144 (MQ=255)
                                                                                               ggTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCtt                                                    >  5:159660/1‑149 (MQ=255)
                                                                                                 ttttATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGa                                                  <  4:24574/149‑1 (MQ=255)
                                                                                                      taGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCgg                                              <  3:301517/148‑1 (MQ=255)
                                                                                                       aGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGc                                             >  8:19038/1‑148 (MQ=255)
                                                                                                          tGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCg                                            >  8:189435/1‑146 (MQ=255)
                                                                                                              cgcgACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCg                                     >  8:299418/1‑149 (MQ=255)
                                                                                                                cgACCGGCGCTTATGCGAGCAAGGTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGc                                   <  8:173899/149‑1 (MQ=255)
                                                                                                                    cGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCAAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACg                               >  1:86139/1‑149 (MQ=255)
                                                                                                                     ggCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAAc                                                       >  4:170195/1‑124 (MQ=255)
                                                                                                                     ggCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAAc                                                       <  3:170195/124‑1 (MQ=255)
                                                                                                                       cgcTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCTAAGAGCGGCAACGTTg                            <  4:244738/149‑1 (MQ=255)
                                                                                                                       cgcTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTg                            <  4:118730/149‑1 (MQ=255)
                                                                                                                       cgcTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTg                            <  2:34092/149‑1 (MQ=255)
                                                                                                                       cgcTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTg                            <  7:189435/149‑1 (MQ=255)
                                                                                                                       cgcTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTg                            <  5:246594/149‑1 (MQ=255)
                                                                                                                       cgcTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGAGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTg                            <  7:299418/149‑1 (MQ=255)
                                                                                                                        gcTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTgg                           <  4:187998/149‑1 (MQ=255)
                                                                                                                         cTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTGGt                          <  4:18680/149‑1 (MQ=255)
                                                                                                                                      tcaactGCGCCGGTTCCAGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCg                                     <  2:192472/120‑1 (MQ=255)
                                                                                                                                      tcaactGCGCCGGTTCCAGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTGGTATGGCTTCTGcc              >  2:191317/6‑148 (MQ=255)
                                                                                                                                      tcaactGCGCCGGTTCCAGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTGGTATGGCTTCTGcc              >  6:135481/6‑148 (MQ=255)
                                                                                                                                      gTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTcc                                                              >  1:279495/1‑100 (MQ=255)
                                                                                                                                      gTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTcc                                                              <  2:279495/100‑1 (MQ=255)
                                                                                                                                       caactGCGCCGGTTCCAGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTGGTATGGCTTCTGCCaa            >  7:230324/5‑149 (MQ=255)
                                                                                                                                           tGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCt                                                            <  8:284523/97‑1 (MQ=255)
                                                                                                                                           tGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCt                                                            >  7:284523/1‑97 (MQ=255)
                                                                                                                                                  gAACCTGCTCAGTTTAGCCATGGCCAGTGGATGTTCGATTATCCGCGTTTGCCGGAGCTGTTGTTCCGTTACGGGGGACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTGGTATGGCTTCTGCCAACAGCGCCtta  <  7:69621/148‑3 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
CCAATCAACAGGCGCAATGGCAAGAAAAGCTGGAGCACATCTATGGCAAGGACGAGTTCGCCCCGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTGGTATGGCTTCTGCCAACAGCGCCTGA  >  NC_002947/1554917‑1555210

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