Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 587,076 C→G 7.9% 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.9% 234.6 / 8.8 76P43R (CCG→CGG) PP_0497Cro/CI family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (33/37);  new base G (4/2);  total (37/39)
Fisher's exact test for biased strand distribution p-value = 4.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.73e-01

CACCTGGATTCCGCATATGAGCGATCACTTCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCAGCAACCCAGTGACCCGATCAGCGAGTTGTTC  >  NC_002947/586933‑587224
                                                                                                                                               |                                                                                                                                                    
cACCTGGATTCCGCATATGAGCGATCACTTCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGACGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACgg                                                                                                                                                 <  4:18324/149‑1 (MQ=255)
  ccTGGATTCCGCATATGAGCGATCACTTCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCt                                                                                                                                               >  1:303919/1‑149 (MQ=255)
   cTGGATTCCGCATATGAGCGATCACTTCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCt                                                                                                                                               >  5:61623/1‑148 (MQ=255)
   cTGGATTCCGCATATGAGCGATCACTTCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCt                                                                                                                                               <  6:61623/148‑1 (MQ=255)
      gATTCCGCATATGAGCGATCACTTCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTc                                                                                                                                            <  4:258235/148‑1 (MQ=255)
        ttCCGCATATGAGCGATCACTTCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCGGACGGCtttgccc                                                                                                                                         <  7:345080/149‑5 (MQ=255)
        ttCCGCATATGAGCGATCACTTCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAAc                                                                                                                                         <  1:131148/149‑1 (MQ=255)
            gCATATGAGCGATCACTTCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGa                                                                                                                                     <  6:321418/149‑1 (MQ=255)
            gCATATGAGCGATCACTTCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGa                                                                                                                                     <  5:39970/149‑1 (MQ=255)
             cATATGAGCGATCACTTCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGaa                                                                                                                                    <  4:74575/149‑1 (MQ=255)
               tatGAGCGATCACTTCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGc                                                                                                                                  <  1:273624/149‑1 (MQ=255)
               tatGAGCGATCACTTCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGc                                                                                                                                  <  2:303919/149‑1 (MQ=255)
                atGAGCGATCACTTCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCg                                                                                                                                 >  1:332059/1‑149 (MQ=255)
                 tGAGCGATCACTTCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGt                                                                                                                                <  3:7876/149‑1 (MQ=255)
                    gCGATCACTTCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGTCGGCTTTCAACCTGAAGCGTATc                                                                                                                             <  8:74675/149‑1 (MQ=255)
                      gATCACTTCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCgg                                                                                                                           >  7:27353/1‑149 (MQ=255)
                            ttCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGAct                                                                                                                      >  7:309613/1‑148 (MQ=255)
                               gCTACAAACCTCAAACTTTCCTTCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAAACTGAAGCGTATCGGCGActtctt                                                                                                                  <  7:229350/149‑1 (MQ=255)
                                 tACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCg                                                                                                                >  6:218955/1‑149 (MQ=255)
                                   cAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCg                                                                                                                <  2:332059/147‑1 (MQ=255)
                                        cTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACttcttc                                                                                                                 <  7:94894/141‑1 (MQ=255)
                                        cTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACttcttc                                                                                                                 >  8:94894/1‑141 (MQ=255)
                                        cTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTc                                                                                                           >  8:68150/1‑147 (MQ=255)
                                        cTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTc                                                                                                           >  1:154221/1‑147 (MQ=255)
                                         tCAAACTTGCCTGGAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCg                                                                                                                <  6:53303/141‑1 (MQ=255)
                                         tCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCggagaacg                                                                                                         >  5:53303/1‑142 (MQ=255)
                                           aaaCTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGa                                                                                                      <  5:157288/149‑1 (MQ=255)
                                               ttGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCg                                                                                                                <  1:267381/135‑1 (MQ=255)
                                               ttGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCg                                                                                                                >  2:267381/1‑135 (MQ=255)
                                                  ccTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGaa                                                                                               <  6:306383/149‑1 (MQ=255)
                                                  ccTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGaa                                                                                               <  5:218955/149‑1 (MQ=255)
                                                    tGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACt                                                                                             >  4:343114/1‑149 (MQ=255)
                                                           aCTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTg                                                                                      >  1:1142/1‑149 (MQ=255)
                                                           aCTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTg                                                                                      >  6:137247/1‑149 (MQ=255)
                                                           aCTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCCAGGATTACGAACTGAACCTg                                                                                      >  3:269817/1‑149 (MQ=255)
                                                                     tATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGa                                                                              <  8:309613/147‑1 (MQ=255)
                                                                     tATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAg                                                                             >  7:4618/1‑148 (MQ=255)
                                                                          ccGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGtt                                                                        >  8:244423/1‑148 (MQ=255)
                                                                          ccGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGtt                                                                        <  7:244423/148‑1 (MQ=255)
                                                                           cGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCg                                                                      >  4:12486/1‑149 (MQ=255)
                                                                                       ccAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTc                                                                       <  7:121505/136‑1 (MQ=255)
                                                                                       ccAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTc                                                                       >  8:121505/1‑136 (MQ=255)
                                                                                          gCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTg                                                       <  8:4618/149‑1 (MQ=255)
                                                                                            tGTCGATAAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGTAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTgcg                                                     >  4:89583/1‑149 (MQ=255)
                                                                                               cGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCgg                                                  <  3:12486/149‑1 (MQ=255)
                                                                                                gATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGa                                                                                                      >  7:255316/1‑96 (MQ=255)
                                                                                                gATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGa                                                                                                      <  8:255316/96‑1 (MQ=255)
                                                                                                 ctCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCTACGGCTTTCAACCTGAAGCGTATCGGCGACTGCTTCGGTCGCGAGGATTACGAACTGAACCTGCGCGCCGAGCAGTTCGCCCGTCTGATCGGCGCGCGGGt                                                <  1:104105/148‑1 (MQ=255)
                                                                                                  tCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTa                                               >  6:50977/1‑149 (MQ=255)
                                                                                                               gTTCAACAAGTACCTGAGCGGGTAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGccc                                  <  3:89583/149‑1 (MQ=255)
                                                                                                                ttCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGcc                                   >  3:46444/1‑147 (MQ=255)
                                                                                                                ttCAAAAAGTACCTGAGCGGGCAAAGCCGTACGACGGCTTTCAACCTGAAGCGTATCGGCTACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTGCGCCCGTCTGATCGTTGCGCGTGTATCGACCCAGGccaa                                 >  4:272236/1‑147 (MQ=255)
                                                                                                                   aacaaGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGActtct                                                                                                                   >  8:125840/1‑64 (MQ=255)
                                                                                                                   aacaaGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGActtct                                                                                                                   <  7:125840/64‑1 (MQ=255)
                                                                                                                   aacaaGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCAGc                               <  7:68150/148‑1 (MQ=255)
                                                                                                                         tACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAg                                                                             >  8:129552/1‑96 (MQ=255)
                                                                                                                         tACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAg                                                                             <  7:129552/96‑1 (MQ=255)
                                                                                                                            cTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCAGCAACCCAGTGa                     >  8:278655/1‑149 (MQ=255)
                                                                                                                              gAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCAGCAACCCAGTGAcc                   >  7:320076/1‑149 (MQ=255)
                                                                                                                               aGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCAGCAACCCAGTGAcc                   <  8:27353/148‑1 (MQ=255)
                                                                                                                                 ggttgaaaGCCGTCGGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCAGCAACCCAGTGACCCGa                >  1:349128/6‑149 (MQ=255)
                                                                                                                                  gggCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCACCCCAGGCCCAGCAACCCAGTGACCCGAt               <  2:1142/149‑1 (MQ=255)
                                                                                                                                   ttgaaaGCCGTCGGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCa                                 >  5:153324/4‑130 (MQ=255)
                                                                                                                                   ttgaaaGCCGTCGGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCa                                 <  6:153324/127‑1 (MQ=255)
                                                                                                                                   ttgaaaGCCGTCGGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCAGCAACCCAGTGACCCGAt               >  7:51205/4‑148 (MQ=255)
                                                                                                                                    gCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCAGCAACCCAGTGACCCGATCa             <  5:137247/149‑1 (MQ=255)
                                                                                                                                      aaaGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCAGCAACCCAGTGAcc                   <  4:234679/141‑1 (MQ=255)
                                                                                                                                      aaaGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCAGCAACCCAGTGAcc                   >  3:234679/1‑141 (MQ=255)
                                                                                                                                      aaaGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCAGCAACCCAGTGACCCGATCAGc           >  7:304911/1‑149 (MQ=255)
                                                                                                                                        aGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCAGCAACCCAGTGACCCGATCAGCGa         <  8:320076/149‑1 (MQ=255)
                                                                                                                                        aGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCAAGCAACCCAGTGACCCGATCAGc           >  2:196827/1‑147 (MQ=255)
                                                                                                                                         gCCGTCGGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTAAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCAGCAACCCAGTGACCCGATCAGCGAg        >  4:162632/1‑149 (MQ=255)
                                                                                                                                          ccGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATc                                                           <  2:159704/97‑1 (MQ=255)
                                                                                                                                          ccGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATc                                                           <  2:165220/97‑1 (MQ=255)
                                                                                                                                          ccGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATc                                                           >  1:165220/1‑97 (MQ=255)
                                                                                                                                          ccGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATc                                                           >  1:159704/1‑97 (MQ=255)
                                                                                                                                               cGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCAGCAACCCAGTGACCCGATCAGCGAGTTGTTc  >  5:92577/1‑149 (MQ=255)
                                                                                                                                               cGACGGAGTTCAAACTGAAGCGTATCGGCGACTTCTTCTTTGTCGATGATTACGAACTTCACCTGCCCGCCGAGCAGTTCGCCCGTCTTATCGGTGCGAGGGTATCGCCCCAGGCCCAGCAACCCAGTGACCCGATCAGCGAgttgtt   <  3:272236/148‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
CACCTGGATTCCGCATATGAGCGATCACTTCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCAGCAACCCAGTGACCCGATCAGCGAGTTGTTC  >  NC_002947/586933‑587224

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: