Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,651,276 G→C 9.2% intergenic (+41/+29) PP_2323 → / ← aroF‑I hypothetical protein/tyr‑sensitive phospho‑2‑dehydro‑3‑deoxyheptonate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,651,2760GC 9.2% 193.5 / 5.4 62intergenic (+41/+29)PP_2323/aroF‑Ihypothetical protein/tyr‑sensitive phospho‑2‑dehydro‑3‑deoxyheptonate aldolase
Reads supporting (aligned to +/- strand):  ref base G (31/25);  new base C (3/3);  total (34/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTCTGGAGTACGCCGAGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGC‑CTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGTCCGTGATCGACACGCCATACTGCAA  >  NC_002947/2651134‑2651421
                                                                                                                                               |                                                                                                                                                 
gCTCTGGAGTACGCCGAGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGC‑CTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCg                                                                                                                                             <  2:31141/149‑1 (MQ=255)
 ctctGGAGTACGCCGAGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGC‑CTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCgg                                                                                                                                            >  7:63764/1‑149 (MQ=255)
  tctGGAGTACGCCGAGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGC‑CTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGc                                                                                                                                           >  8:233055/1‑149 (MQ=255)
   ctGGAGTACGCCGAGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGC‑CTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCg                                                                                                                                          >  2:230717/1‑149 (MQ=255)
      gAGTACGCCGAGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGC‑CTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCACTGCCCGGCGtt                                                                                                                                        <  3:270991/148‑1 (MQ=255)
         tACGCCGAGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGC‑CTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTc                                                                                                                                    <  6:184904/149‑1 (MQ=255)
           cGCCGAGCAGATGGGGTACAGGGTGACCCCGTCTTGCTCATATGTGGAGCGTAAAATGCAGCGCAAGCAGCGGCATTCGACCAAGGC‑CTGAGCTTACGACAGGCACCACACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCa                                                                                                                                   <  4:4246/148‑1 (MQ=255)
                aGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGC‑CTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCg                                                                                                                             <  3:289361/149‑1 (MQ=255)
                                ggTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGC‑CTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCgcgc                                                                                                             >  4:3629/1‑149 (MQ=255)
                                  tGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGC‑CTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCACTGCCCGGCGttt                                                                                                                                       >  6:140653/1‑121 (MQ=255)
                                  tGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGC‑CTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCACTGCCCGGCGttt                                                                                                                                       <  5:140653/121‑1 (MQ=255)
                                   gATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGC‑CTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTa                                                                                                          >  1:200035/1‑149 (MQ=255)
                                     tCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGC‑CTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACg                                                                                                        >  3:113412/1‑149 (MQ=255)
                                       ccGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGC‑CTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGc                                                                                                       <  6:321410/148‑1 (MQ=255)
                                       ccGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGC‑CTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTACCCGGCGTTTTTCATTCCCGAAATGTTCAGCCGCGCTTACGCt                                                                                                      >  3:191819/1‑149 (MQ=255)
                                          tCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGC‑CTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCg                                                                                                   <  7:30918/149‑1 (MQ=255)
                                               cTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGC‑CTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGCGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGgcag                                                                                               <  8:63764/148‑1 (MQ=255)
                                                tCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGC‑CTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGgcagca                                                                                             >  2:165348/1‑149 (MQ=255)
                                                     tgtgGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCGCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATc                                                                                         >  1:78590/1‑149 (MQ=255)
                                                       tgGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGC‑CTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCtt                                                                                      >  6:223923/1‑149 (MQ=255)
                                                       tgGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGC‑CTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCtt                                                                                      <  3:118335/149‑1 (MQ=255)
                                                         gAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGC‑CTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCtt                                                                                      >  3:150480/1‑147 (MQ=255)
                                                          aGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGC‑CTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAg                                                                                   >  4:2190/1‑149 (MQ=255)
                                                               tACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGC‑CTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTgg                                                                              >  2:229988/1‑149 (MQ=255)
                                                                 cATGGAGCGCCAGCAGCGGCATTCGAGCAAGGC‑CTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGgc                                                                             <  4:113412/148‑1 (MQ=255)
                                                                    ggAGCGCCAGCAGCGGCATTCGAGCAAGGC‑CTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTg                                                                               >  2:282214/1‑143 (MQ=255)
                                                                    ggAGCGCCAGCAGCGGCATTCGAGCAAGGC‑CTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTg                                                                               <  1:282214/143‑1 (MQ=255)
                                                                    ggAGCGCCAGCAGCGGCATTCGAGCAAGGC‑CTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGgcatg                                                                          >  7:171096/1‑148 (MQ=255)
                                                                     gAGCGCCAGCAGCGGCATTCGAGCAAGGC‑CTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGgcatgc                                                                         <  8:171096/148‑1 (MQ=255)
                                                                       gcgcCAGCAGCGGCATTCGAGCAAGGC‑CTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGcatgcatg                                                                      >  3:189479/1‑149 (MQ=255)
                                                                         gcCAGCAGCGGCATTCGAGCAAGGC‑CTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCAtgct                                                                    >  4:43193/1‑149 (MQ=255)
                                                                           cagcagCGGCATTCGAGCAAGGC‑CTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTACCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCAtgctgc                                                                  <  4:191819/149‑1 (MQ=255)
                                                                             gcagcGGCATTCGAGCAAGGC‑CTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTgcg                                                                 >  3:188235/1‑148 (MQ=255)
                                                                              cagcGGCATTCGAGCAAGGC‑CTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTgcg                                                                 >  2:245695/1‑147 (MQ=255)
                                                                                 cGGCATTCGAGCAAGGC‑CTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAggg                                                            <  1:230717/149‑1 (MQ=255)
                                                                                 cGGCATTCGAGCAAGGC‑CTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAggg                                                            <  7:233055/149‑1 (MQ=255)
                                                                                  ggCATTCGAGCAAGGC‑CTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCATTGAGCTTGGCATGCATGCTGCGCAGGGt                                                           >  8:52850/1‑149 (MQ=255)
                                                                                     aTTCGAGCAAGGC‑CTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGgcatgca                                                                        >  5:75866/1‑133 (MQ=255)
                                                                                     aTTCGAGCAAGGC‑CTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGgcatgca                                                                        <  6:75866/133‑1 (MQ=255)
                                                                                         gAGCAAGGC‑CTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTtctc                                                     <  3:2190/148‑1 (MQ=255)
                                                                                          aGCAAGGC‑CTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCgg                                                   <  4:150480/149‑1 (MQ=255)
                                                                                           gCAAGGC‑CTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCgg                                                   >  6:279187/1‑148 (MQ=255)
                                                                                            cAAGGC‑CTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTa                                                 <  5:223923/149‑1 (MQ=255)
                                                                                             aaGGC‑CTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAg                                                >  3:124598/1‑149 (MQ=255)
                                                                                               ggC‑CTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGcc                                              <  1:229988/149‑1 (MQ=255)
                                                                                                      aGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCa                                         <  1:165348/148‑1 (MQ=255)
                                                                                                          gacgacAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCTGACCAGTCGa                                    >  2:270849/1‑149 (MQ=255)
                                                                                                             gacAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCACGCGATGc                                 >  4:342614/1‑149 (MQ=255)
                                                                                                             gacAGGCACCAGACACAACAAAAAACGCCGGGCACTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTAGAGCTTGGCATGCATGCTGCGAAGGGTTTTCTCGGAAGCCGACCACTCGATg                                  >  4:156536/1‑148 (MQ=255)
                                                                                                                aGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAgg                              >  4:366408/1‑149 (MQ=255)
                                                                                                                    aCCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGTc                          >  2:139242/1‑149 (MQ=255)
                                                                                                                       aGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCTGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGa                                    >  3:39373/1‑136 (MQ=255)
                                                                                                                       aGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGa                                    <  4:39373/136‑1 (MQ=255)
                                                                                                                        gACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGTCCGTg                      <  5:279187/149‑1 (MQ=255)
                                                                                                                             aacaaAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGGCGACCAGTCGATGCAGGCGTCCGTGATCGa                 >  4:58877/1‑149 (MQ=255)
                                                                                                                               caaAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGt                                       <  7:133710/125‑1 (MQ=255)
                                                                                                                               caaAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGt                                       >  8:133710/1‑125 (MQ=255)
                                                                                                                                   aaaCGCCGGGCACTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGTCCGTGATCGACACGc            >  2:69011/1‑148 (MQ=255)
                                                                                                                                   aaaCCCCCGGCACTGCCCGGGTTTTTTTAATACCCCAATTTTCACCCGCGCTTACGCTGGGGCATCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGTCCGTGATCGACACGcc           <  7:52850/149‑1 (MQ=255)
                                                                                                                                    aaCGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTtc                                                       >  6:130464/1‑104 (MQ=255)
                                                                                                                                    aaCGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTtc                                                       <  5:130464/104‑1 (MQ=255)
                                                                                                                                            gCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGTCCGTGATCGACACGCCATACTGCaa  <  3:3629/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                 
GCTCTGGAGTACGCCGAGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGC‑CTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGTCCGTGATCGACACGCCATACTGCAA  >  NC_002947/2651134‑2651421

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