Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,268,952 T→A 7.7% intergenic (+20/+16) 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,9520TA 7.7% 189.2 / 5.8 65intergenic (+20/+16)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.60e-01

CGCTGCGTGATGGGGACTACTTCAAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCACGATGGCACGCTTGAAGCCGTGCTTGGCCGCTTCCT  >  NC_002947/5268805‑5269098
                                                                                                                                                    |                                                                                                                                                  
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                                                                                                                    >  1:22933/1‑149 (MQ=255)
                               cAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCAt                                                                                                                    >  2:376824/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                                                             >  5:148248/1‑149 (MQ=255)
                                                                                       ggCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGcttc                                                             >  4:79201/1‑149 (MQ=255)
                                                                                           aTGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTc                                                         <  2:22933/149‑1 (MQ=255)
                                                                                           aTGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTc                                                         <  4:121620/149‑1 (MQ=255)
                                                                                            tGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCa                                                                                   <  8:4887/122‑1 (MQ=255)
                                                                                            tGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCa                                                                                   >  7:4887/1‑122 (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                                         <  7:322904/149‑1 (MQ=255)
                                                                                                           cACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTgg                                         <  6:175810/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          <  4:300527/149‑1 (MQ=255)
                                                                                                                                          ggAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCACGATGGCACGCTTGAAGCCGTGCTTGGc          >  2:125866/1‑149 (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)
                                                                                                                                                    |                                                                                                                                                  
CGCTGCGTGATGGGGACTACTTCAAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCACGATGGCACGCTTGAAGCCGTGCTTGGCCGCTTCCT  >  NC_002947/5268805‑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: