Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,268,951 T→A 7.7% intergenic (+19/+17) PP_4643 → / ← radA xanthine/uracil permease family protein/DNA repair enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,268,9510TA 7.7% 190.9 / 4.7 65intergenic (+19/+17)PP_4643/radAxanthine/uracil permease family protein/DNA repair enzyme
Reads supporting (aligned to +/- strand):  ref base T (28/32);  new base A (3/2);  total (31/34)
Fisher's exact test for biased strand distribution p-value = 6.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

CGCGCTGCGTGATGGGGACTACTTCAAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCACGATGGCACGCTTGAAGCCGTGCTTGGCCGCTTCCT  >  NC_002947/5268803‑5269098
                                                                                                                                                     |                                                                                                                                                   
cgcgCTGCGTGATGGGGACTACTTCAAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCt                                                                                                                                                     <  3:354852/149‑1 (MQ=255)
  cgcTGCGTGATGGGGACTACTTCAAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTg                                                                                                                                                   <  4:96474/149‑1 (MQ=255)
        gTGATGGGGACTACTTCAAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTgctcgct                                                                                                                                             <  4:264506/149‑1 (MQ=255)
         tGATGGGGACTACTTCAAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTgctcgctc                                                                                                                                            <  3:56206/149‑1 (MQ=255)
         tGATGGGGACTACTTCAAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTCGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTgctcgctc                                                                                                                                            <  5:68552/149‑1 (MQ=255)
            aggggACTACTTCAAGGCTGGCCAGCTGTGAGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGGCAGGCGCACGGCCGTTGCCGAGCACTTACGGGTTGCAGGAGCTAGCTTGCTCGCTCCt                                                                                                                                          <  3:195241/147‑1 (MQ=255)
              gggACTACTTCAAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGc                                                                                                                                        >  1:96490/1‑149 (MQ=255)
                     ctTCAAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTAtt                                                                                                                                >  6:17554/1‑149 (MQ=255)
                      ttCAAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCAGTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTc                                                                                                                               >  5:374489/1‑149 (MQ=255)
                        cAAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTc                                                                                                                               <  3:143575/147‑1 (MQ=255)
                         aaGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGaa                                                                                                                            >  3:121620/1‑149 (MQ=255)
                            gATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACaa                                                                                                                         >  7:241112/1‑149 (MQ=255)
                             aTGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGCTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAAt                                                                                                                        <  7:366360/149‑1 (MQ=255)
                              tcgCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCCCGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATg                                                                                                                       <  7:172589/147‑1 (MQ=255)
                                gCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGc                                                                                                                      <  5:168500/148‑1 (MQ=255)
                                 cAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCAt                                                                                                                    >  2:376824/1‑149 (MQ=255)
                                 cAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCAt                                                                                                                    >  1:22933/1‑149 (MQ=255)
                                  aaGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATc                                                                                                                   >  8:322904/1‑149 (MQ=255)
                                        ttCGATGGCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAgct                                                                                                              <  8:365477/148‑2 (MQ=255)
                                        ttCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAgcgc                                                                                                             <  1:376824/149‑1 (MQ=255)
                                                cGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGGCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCa                                                                                                     <  4:64866/149‑1 (MQ=255)
                                                 gAGGGTAACGAAGTGCCGATGCTGTAACTGCACGACCACGGCAATGAAACGGTCCGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAg                                                                                                    <  8:241112/149‑1 (MQ=255)
                                                     gTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCgg                                                                                                <  5:17554/149‑1 (MQ=255)
                                                      tAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCggg                                                                                               >  3:300527/1‑149 (MQ=255)
                                                              tGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGc                                                                                        >  2:104300/1‑148 (MQ=255)
                                                                 ccGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGAGCCTGTTCCAGCCGGGTAACGGCAAt                                                                                     >  5:202389/1‑149 (MQ=255)
                                                                         ggAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGc                                                                              >  4:373886/1‑148 (MQ=255)
                                                                         ggAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGc                                                                              >  4:215107/1‑148 (MQ=255)
                                                                          gAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCaa                                                                            >  7:39170/1‑149 (MQ=255)
                                                                                gACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCg                                                                       <  2:96490/148‑1 (MQ=255)
                                                                                   gAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAgg                                                                   >  8:264949/1‑149 (MQ=255)
                                                                                      cACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGgcgc                                                                >  5:175810/1‑149 (MQ=255)
                                                                                        cGGCAATGTAACGGTCAGGCGCACGTCCGTTGCGGTGGACTGACGGGTTGCAGGAGCGAGCTTGCTCGGTCCGGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCTCAGGCGctt                                                              <  5:298496/149‑1 (MQ=255)
                                                                                         ggCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGcttc                                                             >  4:79201/1‑149 (MQ=255)
                                                                                         ggCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGcttc                                                             >  5:148248/1‑149 (MQ=255)
                                                                                             aTGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTc                                                         <  2:22933/149‑1 (MQ=255)
                                                                                             aTGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTc                                                         <  4:121620/149‑1 (MQ=255)
                                                                                              tGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCa                                                                                   >  7:4887/1‑122 (MQ=255)
                                                                                              tGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCa                                                                                   <  8:4887/122‑1 (MQ=255)
                                                                                                 aaCGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAAt                                                                                     <  3:191259/117‑1 (MQ=255)
                                                                                                  aCGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACgg                                                                                         >  5:326244/1‑112 (MQ=255)
                                                                                                  aCGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACgg                                                                                         <  6:326244/112‑1 (MQ=255)
                                                                                                            gcACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCg                                                        >  4:128710/1‑135 (MQ=255)
                                                                                                            gcACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCg                                                        <  3:128710/135‑1 (MQ=255)
                                                                                                            gcACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTg                                          >  1:61940/1‑149 (MQ=255)
                                                                                                            gcACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTg                                          <  2:61940/149‑1 (MQ=255)
                                                                                                             cACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTgg                                         <  6:175810/149‑1 (MQ=255)
                                                                                                             cACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTgg                                         <  7:322904/149‑1 (MQ=255)
                                                                                                                           aGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCACGATGGCACGc                           <  7:347695/149‑1 (MQ=255)
                                                                                                                             cACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCACGATGGCACGCtt                         <  7:264949/149‑1 (MQ=255)
                                                                                                                               cTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCACGATGGCACGCTTGa                       <  4:351467/149‑1 (MQ=255)
                                                                                                                                tGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCACGATGGCACGCTTg                        >  3:210870/1‑147 (MQ=255)
                                                                                                                                      cgtGCAGGAGCGAGCAAGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCACGATGGCACGCTTGa                       >  3:294368/3‑142 (MQ=255)
                                                                                                                                      cgtGCAGGAGCGAGCAAGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCACGATGGCACGCTTGa                       <  4:294368/140‑1 (MQ=255)
                                                                                                                                         gCAGGAGCGAGCAAGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCACGATGGCACGCTTGa                       <  3:50859/139‑1 (MQ=255)
                                                                                                                                         gCAGGAGCGAGCAAGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCACGATGGCACGCTTGa                       >  4:50859/1‑139 (MQ=255)
                                                                                                                                          cAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCACGATGGCACGCTTGAAGCCGTGCtt             <  2:193705/148‑1 (MQ=255)
                                                                                                                                          cAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCACGATGGCACGCTTGAAGCCGTGCTTg            >  8:81360/1‑149 (MQ=255)
                                                                                                                                           aGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCACGATGGCACGCTTGAAGCCGTGCTTgg           >  5:354170/1‑149 (MQ=255)
                                                                                                                                            ggAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAGCAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCACGATGGCACg                            >  6:53171/1‑131 (MQ=255)
                                                                                                                                            ggAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAGCAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCACGATGGCACg                            <  5:53171/131‑1 (MQ=255)
                                                                                                                                            ggAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCACGATGGCACGCTTGAAGCCGTGCTTGGc          >  2:125866/1‑149 (MQ=255)
                                                                                                                                            ggAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCACGATGGCACGCTTGAAGCCGTGCTTGGc          <  4:300527/149‑1 (MQ=255)
                                                                                                                                             gagcgagcAAGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCACGATGGCACGCTTGAAGCCGTGCTTGGcc         >  5:204588/1‑149 (MQ=255)
                                                                                                                                                cgagcTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCACGATGGCACGCTTGAAGCCGTGCTTGGCCGct      >  3:289713/1‑149 (MQ=255)
                                                                                                                                                  agcTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCACGATGGCACGCTTGAAGCCGTGCTTGGCCGcttc    <  4:325828/149‑1 (MQ=255)
                                                                                                                                                     ttGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCACGATGGCACGCTTGAAGCCGTGCTTGGCCGcttcct  >  5:320565/1‑148 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CGCGCTGCGTGATGGGGACTACTTCAAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCACGATGGCACGCTTGAAGCCGTGCTTGGCCGCTTCCT  >  NC_002947/5268803‑5269098

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: