Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 587,076 C→G 7.3% P43R (CCG→CGG)  PP_0497 → Cro/CI family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947587,0760CG 7.3% 211.2 / 6.2 68P43R (CCG→CGG) PP_0497Cro/CI family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (29/34);  new base G (3/2);  total (32/36)
Fisher's exact test for biased strand distribution p-value = 6.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.58e-01

GCGTCACCTGGATTCCGCATATGAGCGATCACTTCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCAGCAACCCAGTGACCCGATCAGCGAGTT  >  NC_002947/586929‑587220
                                                                                                                                                   |                                                                                                                                                
gcgTCACCTGGATTCCGCATATGAGCGATCACTTCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTcc                                                                                                                                                  <  1:100798/148‑1 (MQ=255)
        tGGATTCCGCATATGAGCGATCACTTCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCtt                                                                                                                                          >  7:213036/1‑148 (MQ=255)
           aTTCCGCATATGAGCGATGACTTCGCTACAAACCTCAAACTTGACTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCGGACGGCtttgc                                                                                                                                       <  8:290602/148‑3 (MQ=255)
           aTTCCGCATATGAGCGATCACTTCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCaa                                                                                                                                      <  6:299781/149‑1 (MQ=255)
            ttcagaatatGAGCGATCACTTCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCTCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCaa                                                                                                                                      <  3:195128/142‑1 (MQ=255)
            ttCCGCATATGAGCGATCACTTCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAAc                                                                                                                                     <  4:90054/149‑1 (MQ=255)
               cGCATATGAGCGATCACTTCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTACCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCt                                                                                                                                   >  5:46403/1‑148 (MQ=255)
                 cATATGAGCGATCACTTCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTg                                                                                                                                  >  8:34734/1‑147 (MQ=255)
                  atatGAGCGATCACTTCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAg                                                                                                                               <  8:213036/149‑1 (MQ=255)
                         cGATCACTTCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATATCCGAAGTTAGCCGCCAGCAGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCg                                                                                                                        >  5:248067/1‑149 (MQ=255)
                           aTCACTTCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGc                                                                                                                      <  3:288019/149‑1 (MQ=255)
                               cTTCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGAct                                                                                                                  >  2:155443/1‑149 (MQ=255)
                                ttCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAAc                                                                                                                                     <  3:266484/129‑1 (MQ=255)
                                ttCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAAc                                                                                                                                     >  4:266484/1‑129 (MQ=255)
                                ttCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGActt                                                                                                                 >  2:7433/1‑149 (MQ=255)
                                    cTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAGCCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACttcttc                                                                                                             <  5:285720/149‑1 (MQ=255)
                                       cAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCg                                                                                                                        >  3:112924/1‑135 (MQ=255)
                                       cAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCg                                                                                                                        <  4:112924/135‑1 (MQ=255)
                                          aCCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTc                                                                                                       >  4:275931/1‑149 (MQ=255)
                                            cTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCTGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTc                                                                                                       >  1:215157/1‑147 (MQ=255)
                                            cTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTc                                                                                                       >  2:145099/1‑147 (MQ=255)
                                               aaaCTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAg                                                                                                    <  1:155443/147‑1 (MQ=255)
                                               aaaCTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAg                                                                                                    <  5:102431/147‑1 (MQ=255)
                                                 aCTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGAtt                                                                                                <  2:215157/149‑1 (MQ=255)
                                                     gCCTGCAGTCACTACCGCTCTATTACCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTAc                                                                                              <  6:46403/147‑1 (MQ=255)
                                                       cTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAAc                                                                                          <  6:25208/149‑1 (MQ=255)
                                                                 tACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAg                                                                                                    <  3:287512/129‑1 (MQ=255)
                                                                 tACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAg                                                                                                    >  4:287512/1‑129 (MQ=255)
                                                                  aCCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTAACTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGa                                                                                                                                 >  8:50989/1‑99 (MQ=255)
                                                                  aCCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTAACTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGa                                                                                                                                 <  7:50989/99‑1 (MQ=255)
                                                                   ccGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCccg                                                                              >  7:196288/1‑149 (MQ=255)
                                                                           tttcataaGTTTGCCTCCAGCTGTCGATCAACAGCGCGCAGATCAAAAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCa                                                                       <  8:136813/142‑1 (MQ=255)
                                                                             tCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGt                                                                     >  8:253282/1‑148 (MQ=255)
                                                                               cGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCg                                                                  >  2:173769/1‑149 (MQ=255)
                                                                                gAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGc                                                                 <  7:34734/149‑1 (MQ=255)
                                                                                      tGCCTCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTc                                                            <  6:124435/148‑1 (MQ=255)
                                                                                       gccgccAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCt                                                           >  4:131746/1‑148 (MQ=255)
                                                                                                tGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCg                                                      >  5:220664/1‑144 (MQ=255)
                                                                                                tGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCg                                                      <  6:220664/144‑1 (MQ=255)
                                                                                                 gTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTgcgc                                                >  7:26345/1‑149 (MQ=255)
                                                                                                  tCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTgcgc                                                <  7:253282/148‑1 (MQ=255)
                                                                                                  tCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTgcgc                                                <  3:131746/148‑1 (MQ=255)
                                                                                                    gATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAAc                                                                                     <  7:166678/109‑1 (MQ=255)
                                                                                                    gATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAAc                                                                                     >  8:166678/1‑109 (MQ=255)
                                                                                                    gATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGTCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCgg                                              <  3:275931/148‑1 (MQ=255)
                                                                                                       cAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTAt                                          >  6:143740/1‑149 (MQ=255)
                                                                                                              gcgcAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCa                                   <  8:196288/149‑1 (MQ=255)
                                                                                                                   tggcggcaaGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGcc                               <  7:39004/142‑1 (MQ=255)
                                                                                                                   gTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGccc                              <  1:7433/149‑1 (MQ=255)
                                                                                                                   gTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTCATCGGTGCGCGGGTATCGCCCCAGGccc                              <  1:173769/149‑1 (MQ=255)
                                                                                                                    ttCAACAAGTACCTGAGCTGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCGGCCGAGCAGTTGGCCCGTGTGATCGGTGCGCGGGTATCGCCCCAGGCCCa                             <  6:248067/149‑1 (MQ=255)
                                                                                                                       aacaaGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCAGc                           >  3:149544/1‑148 (MQ=255)
                                                                                                                          aaGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCAGCAAcc                       >  4:100455/1‑149 (MQ=255)
                                                                                                                                   aGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCa                             <  3:192534/134‑1 (MQ=255)
                                                                                                                                   aGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCa                             >  4:192534/1‑134 (MQ=255)
                                                                                                                                     cGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGACCAGTTCGCCCGTCTGATCGGTGCGCGGGAATCGCCCCACGCCCa                             >  4:182271/1‑132 (MQ=255)
                                                                                                                                     cGGGCAAAGCAGTCCGACGGCTTTCAACCTGAAGCGTATCGGCTACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCa                             <  3:182271/132‑1 (MQ=255)
                                                                                                                                      gttgaaaGCCGTCGGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCAGCAACCCAGTGACCCGAt           >  5:220259/5‑149 (MQ=255)
                                                                                                                                         gaaaGCCGTCGGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCAGCAACCCAGTGACCCGATCAg        >  5:199156/2‑149 (MQ=255)
                                                                                                                                         cAAAGCCGTCCGTCGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCAGCAACCCAGTGACCCGATCAg        <  1:145099/149‑1 (MQ=255)
                                                                                                                                          aaaGCCGTCGGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCAGCAACCCAGTGACCCGATCAGc       >  5:45845/1‑149 (MQ=255)
                                                                                                                                          aaaGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCa                             <  7:205933/127‑1 (MQ=255)
                                                                                                                                          aaaGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCa                             >  8:205933/1‑127 (MQ=255)
                                                                                                                                           aaGCCGTCGGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCAGCAACCCAGTGACCCGATCAGCg      <  6:220259/149‑1 (MQ=255)
                                                                                                                                            aGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCAGCAACCCAGTGAcc               <  4:282877/139‑1 (MQ=255)
                                                                                                                                            aGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCAGCAACCCAGTGAcc               >  3:282877/1‑139 (MQ=255)
                                                                                                                                               cGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCAGCAACCCAGTGACCCGATCAGCGAgtt  >  6:90106/1‑149 (MQ=255)
                                                                                                                                                gTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCa                             >  3:95306/1‑121 (MQ=255)
                                                                                                                                                gTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCa                             <  4:95306/121‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
GCGTCACCTGGATTCCGCATATGAGCGATCACTTCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCAGCAACCCAGTGACCCGATCAGCGAGTT  >  NC_002947/586929‑587220

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