Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,415,037 G→A 8.7% intergenic (+146/+22) PP_5644 → / ← PP_3899 membrane protein/hicB family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,415,0370GA 8.7% 192.7 / 7.2 69intergenic (+146/+22)PP_5644/PP_3899membrane protein/hicB family protein
Reads supporting (aligned to +/- strand):  ref base G (31/32);  new base A (3/3);  total (34/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.38e-01

TGGCTATGAAGCACTCCTCTTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAGGTCGTGCCCTACCCTCACATTGAAC  >  NC_002947/4414893‑4415182
                                                                                                                                                |                                                                                                                                                 
tGGCTATGAAGCACTCCTCTTTAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCACTGCAATGCg                                                                                                                                               <  5:23709/149‑1 (MQ=255)
          gCACTCCTCTTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCtt                                                                                                                                     <  1:121827/149‑1 (MQ=255)
            aCTCCTCTTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTg                                                                                                                                   <  7:270177/149‑1 (MQ=255)
             ctcctcTTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTgg                                                                                                                                  <  8:173279/149‑1 (MQ=255)
             ctcctcTTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTgg                                                                                                                                  >  7:215504/1‑149 (MQ=255)
               cctcTTGAACCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGggc                                                                                                                                <  3:232692/149‑1 (MQ=255)
                ctctTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGggc                                                                                                                                <  3:222411/148‑1 (MQ=255)
                 tctTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGggcgg                                                                                                                              >  8:54139/1‑149 (MQ=255)
                       aGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTAcgc                                                                                                                        <  5:190599/149‑1 (MQ=255)
                           gccATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGtga                                                                                                                    <  7:213566/149‑1 (MQ=255)
                            ccATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCtt                                                                                                                                     <  5:265093/131‑1 (MQ=255)
                            ccATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCtt                                                                                                                                     >  6:265093/1‑131 (MQ=255)
                            ccATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGtgat                                                                                                                   >  6:280009/1‑149 (MQ=255)
                                 aaGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGAcc                                                                                                               >  7:261147/1‑148 (MQ=255)
                                  aGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGAcc                                                                                                               >  8:188370/1‑147 (MQ=255)
                                   gTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCgg                                                                                                                                              >  4:5332/1‑115 (MQ=255)
                                   gTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCgg                                                                                                                                              <  3:5332/115‑1 (MQ=255)
                                   gTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAg                                                                                                            <  4:11682/149‑1 (MQ=255)
                                   gTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTATTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACCCGTGATGACCCAg                                                                                                            >  8:275739/1‑149 (MQ=255)
                                    tGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTATTCTTTGGGCGGCTACGCGTGATGACCCAgg                                                                                                           <  5:90400/149‑1 (MQ=255)
                                     gATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAgg                                                                                                           <  3:292159/148‑1 (MQ=255)
                                     gATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAgg                                                                                                           >  7:65744/1‑148 (MQ=255)
                                             aGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGtga                                                                                                                    <  3:284592/131‑1 (MQ=255)
                                             aGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGtga                                                                                                                    >  4:284592/1‑131 (MQ=255)
                                                  gCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAgcgc                                                                                             >  1:218705/1‑149 (MQ=255)
                                                  gCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGGGATGACCCAGGTACTCACTAAgcgc                                                                                             <  5:254739/149‑1 (MQ=255)
                                                   ccTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCt                                                                                            >  8:209163/1‑149 (MQ=255)
                                                    cTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCtt                                                                                           <  7:188370/149‑1 (MQ=255)
                                                       ggcggGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGc                                                                                        <  8:65744/149‑1 (MQ=255)
                                                            gCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTc                                                                                   <  2:218705/149‑1 (MQ=255)
                                                                  ccATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGAccc                                                                                                              >  6:64894/1‑116 (MQ=255)
                                                                  ccATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGAccc                                                                                                              <  5:64894/116‑1 (MQ=255)
                                                                   cATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCAt                                                                            <  1:252804/149‑1 (MQ=255)
                                                                       cATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTcc                                                                         <  8:215504/148‑2 (MQ=255)
                                                                       cATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTc                                                                          >  1:13140/1‑147 (MQ=255)
                                                                         ttCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAgg                                                                      <  5:280009/149‑1 (MQ=255)
                                                                                      cACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTAcgc                                                                                                                        >  1:240075/1‑86 (MQ=255)
                                                                                      cACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTAcgc                                                                                                                        <  2:240075/86‑1 (MQ=255)
                                                                                      cACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCACGTACTCACTAAGCGCTTTGCGTTTCAGATCATTCAAGAAAATACtttcc                                                           >  8:236608/1‑145 (MQ=255)
                                                                                               cccTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGc                                                >  4:233458/1‑149 (MQ=255)
                                                                                                  tGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCAc                                             >  1:50997/1‑149 (MQ=255)
                                                                                                  tGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTAGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCAc                                             >  1:107293/1‑149 (MQ=255)
                                                                                                   gCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGAt                                                                               >  5:210281/1‑114 (MQ=255)
                                                                                                   gCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGAt                                                                               <  6:210281/114‑1 (MQ=255)
                                                                                                       aGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTgg                                        <  2:107293/149‑1 (MQ=255)
                                                                                                         tGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGcc                                      >  7:178195/1‑149 (MQ=255)
                                                                                                         tGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCACTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGcc                                      >  5:230030/1‑149 (MQ=255)
                                                                                                           aaGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCg                                     <  2:248594/148‑1 (MQ=255)
                                                                                                             ggCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTa                                  <  7:275739/149‑1 (MQ=255)
                                                                                                                tCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGt                               >  3:316472/1‑149 (MQ=255)
                                                                                                                    ccGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAg                           >  7:200242/1‑149 (MQ=255)
                                                                                                                      gCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAGGt                         >  4:215591/1‑149 (MQ=255)
                                                                                                                      gCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAGGt                         >  7:14167/1‑149 (MQ=255)
                                                                                                                           ccTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATAGTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAGGTCGTg                     <  8:261147/148‑1 (MQ=255)
                                                                                                                            cTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAGGTCGTGcc                   >  2:8743/1‑149 (MQ=255)
                                                                                                                            cTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAGGTCGTGcc                   <  4:61358/149‑1 (MQ=255)
                                                                                                                            cTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAGGTCGTGc                    >  3:61358/1‑148 (MQ=255)
                                                                                                                             tACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGTGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAGGTCGTGccc                  >  1:46001/1‑149 (MQ=255)
                                                                                                                                 ttAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAGGTCGTGCCCTAcc              >  1:126628/1‑149 (MQ=255)
                                                                                                                                 ttAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAGGTCGTGCCCTAcc              <  7:54139/149‑1 (MQ=255)
                                                                                                                                  aaGCCCGCACTGCAATGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAGGTCGTGCCCTAcc              <  3:117629/147‑1 (MQ=255)
                                                                                                                                  aaGCCCGCACTGCAATGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAGGTCGTGCCCTAcc              >  4:117629/2‑148 (MQ=255)
                                                                                                                                    gCCCGCACTGCAATGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAGGTCGTGCCCTACCCt            >  4:141437/1‑148 (MQ=255)
                                                                                                                                    gCCCGCACTGCAATGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAGGTCGTGCCCTACCCt            >  3:197299/1‑148 (MQ=255)
                                                                                                                                     cccGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAGGTCGTGCCCTACCCTca          <  2:13140/149‑1 (MQ=255)
                                                                                                                                      ccGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAGGTCGTGCCCTACCCTca          >  5:110993/1‑148 (MQ=255)
                                                                                                                                       cGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGAGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAGGTCGTGCCCTACCCTcaca        <  7:236608/149‑1 (MQ=255)
                                                                                                                                           ctGCAATGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAGGTCGTGCCCTACCCTCACATTGa    <  3:141437/148‑1 (MQ=255)
                                                                                                                                              gaGTTCGGGCTTTTCTTTTGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAGGTCGTGCCCTACCCTCACATTGAAc  <  8:200242/147‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
TGGCTATGAAGCACTCCTCTTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAGGTCGTGCCCTACCCTCACATTGAAC  >  NC_002947/4414893‑4415182

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: