Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 272,007 A→T 11.4% intergenic (+12/+17) PP_0218 → / ← metP sensory box protein/L,D‑methionine D‑methionine ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947272,0070AT11.4% 161.0 / 9.8 61intergenic (+12/+17)PP_0218/metPsensory box protein/L,D‑methionine D‑methionine ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (24/30);  new base T (3/4);  total (27/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

GGATGCAGGTGCATGCCGAGGGCATCGAGCAGGCGGAACAGGCCAGCTTCCTGCTGCAGGAGGAGTGCCAGTTGGGGCAGGGGTACTGGTTCGGGCGGCCAGTGCCGGCTGAGGATCTTCGCTGGGATTGAGACCGGGCTGCAAAGCAGCCCTCAGATTTCAACGCTTGTTCAAACCCTTCGCCAGGCGGTCCCCGCCCAGCTGGATCAACGCCACTAACGCTACCAGCATGGCGATCACGGTCAGCATGATCTGGCTGTCGAAACGCTGGTAGCCATACCGATAGGCGA  >  NC_002947/271865‑272154
                                                                                                                                                |                                                                                                                                                   
ggATGCAGGTGCATGCCGAGGGCATCGAGCAGGCGGAACAGGCCAGCTTCCTGCTGCAGGAGGAGTGCCAGTTGGGGCAGGGGTACTGGTTCGGGCGGCCAGTGCCGGCTGAGGATCTTCGCTGGGATTGAGACCGGGCTGCAAagcag                                                                                                                                               <  1:51365/149‑1 (MQ=255)
 gATGCAGGTGCATGCCGAGGGCATCGAGCAGGCGGAACAGGCCAGCTTCCTGCTGCAGGAGGAGTGCCAGTTGGGGCAGGGGTACTGGTTCGGGCGGCCAGTGCCGGCTGAGGATCTTCGCTGGGATTGAGACCGGGCTGCAAagcagc                                                                                                                                              >  1:160362/1‑149 (MQ=255)
   tGCAGGTGCATGCCGAGGGCATCGAGCAGGCGGAACAGGCCAGCTTCCTGCTGCAGGAGGAGTGCCAGTTGGGGCAGGGGTACTGGTTCGGGCGGCCAGTGCCGGCTGAGGATCTTCGCTGGGATTGAGACCGGGCTGCAAAGCAGccc                                                                                                                                            <  1:27832/149‑1 (MQ=255)
     cAGGTGCATGCCGAGGGCATCGAGCAGGCGGAACAGGCCAGCTTCCTGCTGCAGGAGGAGTGCCAGTTGGGGCAGGGGTACTGGTTCGGGCGGCCAGTGCCGGCTGAGGATCTTCGCTGGGATTGAGACCGGGCTGCAAAGCAGCCCt                                                                                                                                           >  2:219289/1‑148 (MQ=255)
      aGGTGCATGCCGAGGGCATCGAGCAGGCGGAACAGGCCAGCTTCCTGCTGCAGGAGGAGTGCCAGTTGGGGCAGGGGTACTGGTTCGGGCGGCCAGTGCCGGCTGAGGATCTTCGCTGGGATTGAGACCGGGCTGCTTTGCAGcccggt                                                                                                                                         <  2:11782/149‑4 (MQ=255)
       ggTGCATGCCGAGGGCATCGAGCAGGCGGAACAGGCCAGCTTCCTGCTGCAGGAGGAGTGCCAGTTGGGGCAGGGGTACTGGTTCGGGCGGCCAGTGCCGGCTGAGGATCTTCGCTGGGATTGAGACCGGGCTGCTTTGCAGcccggtc                                                                                                                                        <  1:219124/149‑5 (MQ=255)
             tGCCGAGGGCATCGAGCAGGCGGAACAGGCCAGCTTCCTGCTGCAGGAGGAGTGCCAGTTGGGGCAGGGGTACTGGTTCGGGCGGCCAGTGCCGGCTGAGGATCTTCGCTGGGATTGAGACCGGGCTGCAAAGCAGCCCTCAGATTTc                                                                                                                                   >  7:236067/1‑148 (MQ=255)
             tGCCGAGGGCATCGAGCAGGCGGAACAGGCCAGCTTCCTGCTGCAGGAGGAGTGCCAGTTGGGGCAGGGGTACTGGTTCGGGCGGCCAGTGCCGGCTGAGGATCTTCGCTGGGATTGAGACCGGGCTGCAAAGCAGCCCTCAGATTTCa                                                                                                                                  <  7:86777/149‑1 (MQ=255)
              gCCGAGGGCATCGAGCAGGCGGAACAGGCCAGCTTCCTGCTGCAGGAGGAGTGCCAGTTGGGGCAGGGGTACTGGTTCGGGCGGCCAGTGCCGGCTGAGGATCTTCGCTGGGATTGAGACCGGGCTGCAAAGCAGCCCTCAGATTTc                                                                                                                                   >  8:62494/1‑147 (MQ=255)
              gCCGAGGGCATCGAGCAGGCGGAACAGGCCAGCTTCCTGCTGCAGGAGGAGTGCCAGTTGGGGCAGGGGTACTGGTTCGGGCGGCCAGTGCCGGCTGAGGATCTTCGCTGGGATTGAGACCGGGCTGCAAAGCAGCCCTCAGATTTc                                                                                                                                   <  7:62494/147‑1 (MQ=255)
               ccGAGGGCATCGAGCGGGCGGAACAGGCCAGCTTCCTGCTGCAGGAGGAGTGCCAGTTGGGGCAGGGGTACTGGTTCGGGCGGCCAGTGCCGGCTGAGGATCTTCGCTGGGATTGAGACCGGGCTGCAAAGCAGCCCTCAGATTTCaa                                                                                                                                 <  8:294657/148‑1 (MQ=255)
               ccGAGGGCATCGAGCAGGCGGAACAGGCCAGCTTCCTGCTGCAGGAGGAGTGCCAGTTGGGGCAGGGGTACTGGTTCGGGCGGCCAGTGCCGGCTGAGGATCTTCGCTGGGATTGAGACCGGGCTGCAAAGCAGCCCTCAGATTTCaa                                                                                                                                 <  7:251848/148‑1 (MQ=255)
                   gggCATCGAGCAGGCGGAACAGGCCAGCATTCCTGCTGCAGGAGGAGTGCCAGTTGGGGCAGGGGTACTGGTTCGGGCGGCCAGTGCCGGCTGAGGATCTTCGCTGGGATTGAGACCGGGCTGCAAAGCAGCCCTCAGATTTCAACGCt                                                                                                                             >  1:25217/1‑149 (MQ=255)
                      cATCGAGCAGGCGGAACAGGCCAGCTTCCTGCTGCAGGAGGAGTGCCAGTTGGGGCAGGGGTACTGGTTCGGGCGGCCAGTGCCGGCTGAGGATCTTCGCTGGGATTGAGACCGGGCTGCAAAGCAGCCCTCAGATTTCAACGCTTGt                                                                                                                          >  3:174110/1‑148 (MQ=255)
                       aTCGAGCAGGCGGAACAGGCCAGCTTCCTGCTGCAGGAGGAGTGCCCGTTGGGGCAGGGGTACTGGTTCGGGCGGCCAGTGCCGGCTGAGGATCTTCGCTGGGATTGAGACCGGGCTGCAAAGCAGCCCTAAGATTTAAACGCTTGtt                                                                                                                         >  5:400427/1‑148 (MQ=255)
                                gCGGAACAGGCCAGCTTCCTGCTGCAGGAGGAGTGCCAGTTGGGGCAGGGGTACTGGTTCGGGCGGCCAGTGCCGGCTGAGGATCTTCGCTGGGATTGAGACCGGGGTGCAAAGCAGCCCTCa                                                                                                                                         >  1:135663/1‑123 (MQ=255)
                                gCGGAACAGGCCAGCTTCCTGCTGCAGGAGGAGTGCCAGTTGGGGCAGGGGTACTGGTTCGGGCGGCCAGTGCCGGCTGAGGATCTTCGCTGGGATTGAGACCGGGCTGCAAAGCAGCCCTCa                                                                                                                                         <  2:135663/123‑1 (MQ=255)
                                      cAGGCCAGCTTCCTGCTGCAGGAGGAGTGCCAGTTGGGGCAGGGGTACTGGTTCGGGCGGCCAGTGCCGGCTGAGGATCTTCGCTGGGATTGAGACCGGGCTGCAAAGCAGCCCTCAGATTTCAACGCTTGTTCAAACCCTTCGCCAgg                                                                                                         <  5:331494/149‑1 (MQ=255)
                                       aGGCCAGCTTCCTGCTGCAGGAGGAGTGCCAGTTGGGGCAGGGGTACTGGTTCGGGCGGCCAGTGCCGGCTGAGGATCTTCGCTGGGATTGAGACCGGGCTGCAAAGCAGCCCTCAGATTTCAACGCTTGTTCAAACCCTTCGCCAGGc                                                                                                        <  4:169060/149‑1 (MQ=255)
                                        ggCCAGCTTCCTGCTGCAGGAGGAGTGCCAGTTGGGGCAGGGGTACTGGTTCGGGCGGCCAGTGCCGGCTGAGGATCTTCGCTGGGATTGAGACCGGGCTGCATAGCAGCCCTCAGATTTCAACGCTTGTTCAAACCCTTCGCCAGGCg                                                                                                       <  1:288331/149‑1 (MQ=255)
                                        ggCCAGCTTCCTGCTGCAGGAGGAGTGCCAGTTGGGGCAGGGGTACTGGTTCGGGCGGCCAGTGCCGGCGGAGGATCTTCGCTGGGATTGAGACCGGGCTGCAAAGCAGCCCTCAGATTTCAACGGTTGTTCAAACCCTTCGCCAGGCg                                                                                                       >  3:305113/1‑149 (MQ=255)
                                         gCCAGCTTCCTGCTGCAGGAGGAGTGCCAGTTGGGGCAGGGGTACTGGTTCGGGCGGCCAGTGCCGGCTGAGGATCTTCGCTGGGATTGAGACCGGGCTGCAAAGCAGCCCTCAGATTTCAACGCTTGTTCAAACCCTTCGCCAGGCgg                                                                                                      >  1:22380/1‑149 (MQ=255)
                                             gCTTCCTGCTGCAGGAGGAGTGCCAGTTGGGGCAGGGGTACTGGTTCGGGCGGCCAGTGCCGGCTGAGGATCTTCGCTGGGATTGAGACCGGGCTGCAAAGCAGCCCTCAGATTTCAACGCTTGTTCAAACCCTTCGCCAGGCGGTccc                                                                                                  <  6:83541/149‑1 (MQ=255)
                                                 tCCTGCTGCAGGAGGAGTGCCAGTTGGGGCAGGGGTACTGGTTCGGGCGGCCAGTGCCGGCTGAGGATCTTCGCTGGGATTGAGACCGGGCTGCAAAGCAGCCCTCAGATTTCAACGCTTGTTCAAACCCTTCGCCAGGCGGTCCCCGc                                                                                               >  8:147387/1‑149 (MQ=255)
                                                    tgctgcAGGAGGAGTGCCAGTTGGGGCAGGGGTACTGGTTCGGGCGGCCAGTGCCGGCTGAGGATCTTCGCTGGGATTGAGACCGGGCTGCaaa                                                                                                                                                   >  2:384046/1‑94 (MQ=255)
                                                    tgctgcAGGAGGAGTGCCAGTTGGGGCAGGGGTACTGGTTCGGGCGGCCAGTGCCGGCTGAGGATCTTCGCTGGGATTGAGACCGGGCTGCaaa                                                                                                                                                   <  1:384046/94‑1 (MQ=255)
                                                        gcAGGAGGAGTGCCAGTTGGGGCAGGGGTACTGGTTCGGGCGGCCAGTGCCGGCTGAGGATCTTCGCTGGGATTGAGACCGGGCTGCAAAGCAGCCCTCAGATTTCAACGCTTGTTCAAACCCTTCGCCAGGCGGTCCCCGCCCAGCTg                                                                                        >  6:319892/1‑149 (MQ=255)
                                                         cAGGAGGAGTGCCAGTTGGGGCAGGGGTACTGGTTCGGGCGGCCAGTGCCGGCTGAGGATCTTCGCTGGGATTGAGACCGGGCTGCAAAGCAGCCCTCAGATTTCAACGCTTGTTCAAACCCTTCGCCAGGCGGTCCCCGCCCAGCTgg                                                                                       <  2:22380/149‑1 (MQ=255)
                                                           ggaggaGTGCCAGTTGGGGCAGGGGTACTGGTTCGGGCGGCCAGTGCCGGCTGAGGATCTTCGCTGGGATTGAGACCGGGCTGCAAAGCAGCCCTCAGATTTCAACGCTTGTTCAAACCCTTCGCCAGGCGGTCCCCGCCCAGCTGGAt                                                                                     <  7:147387/149‑1 (MQ=255)
                                                                 gTGCCAGTTGGGGCAGGGGTACTGGTTCGGGCGGCCAGTGCCGGCTGAGGATCTTCGCTGGGATTGAGACCGGGCTGCAAAGCAGCCCTCAGATTTCAACGCTTGTTCAAACCCTTCGCCAGGCGGTCCCCGCCCAGCTGGATCAACGc                                                                               >  5:405055/1‑149 (MQ=255)
                                                                  tGCCAGTTGGGGCAGGGGTACTGGTTCGGGCGGCCAGTGCCGGCTGAGGATCTTCGCTGGGATTGAGACCGGGCTGCAAAGCAGCCCTCAGATTTCAACGCTTGTTCAAAACCTTCGCCAGGCGGTCCCCGCCGAGc                                                                                          >  6:191597/1‑137 (MQ=255)
                                                                    ccAGTTGGGGCAGGGGTACTGGTTCGGGCGGCCAGTGCCGGCTGAGGATCTTCGCTGGGATTGAGACCGGGCTGCAAAGCAGCCCTCAGATTTCAACGCTTGTTCAAACCCTTCGCCAGGCGGTCCCCGCCCAGCTGGATCAACGCCa                                                                             <  5:319892/148‑1 (MQ=255)
                                                                           gggCAGTGGTACTGGTTCGGGCGGCCAGTGCCGGCTGAGGATCTTCGCTGGGATTGAGACCGGGCTGCAAAGCAGCCCTCAGATTTCAACGCTTGTTCAAACCCTTCGCCAGGCGGTCCCCGCCCAGCTGGATCAACGCCACTAACGc                                                                      <  4:271916/148‑1 (MQ=255)
                                                                           gggCAGGGGTACTGGTTCGGGCGGCCAGTGCCGGCTGAGGATCTTCGCTGGGATTGAGACCGGGCTGCAAAGCAGCCCTCAGATTTCAACGCTTGTTCAAACCCTTCGCCAGGCGGTCCCCGCCCAGCTGGATCAACGCCACTAACGCt                                                                     <  4:174110/149‑1 (MQ=255)
                                                                           gggCAGGGGTACTGGTTCGGGCGGCCAGTGCCGGCTGAGGATCTTCGCTGGGATTGAGACCGGGCTGCAAAGCAGCCCTCAGATTTCAACGCTTGTTCAAACCCTTCGCCAGGCGGTCCCCCCCCAGCTGGATCAACGCCACTAACGc                                                                      <  1:219289/148‑1 (MQ=255)
                                                                                  ggTACTGGTTCGGGCGGCCAGTGCCGGCTGAGGATCTTCGCTGGGATTGAGACCGGGCTGCAAAGCAGCCCTCAGATTTCAACGCTTGTTCAAACCCTTCGCCAGGCGGTCCCCGCCCAGCTGGATCAACGCCACTAACGCTACCAGCa                                                              >  6:397358/1‑149 (MQ=255)
                                                                                   gTACTGGTTCGGGCGGCCAGTGCCGGCTGAGGATCTTCGCTGGGATTGAGACCGGGCTGCAAAGCAGCCCTCAGATTTCAACGCTTGTTCAAACCCTTCg                                                                                                              >  1:285426/1‑100 (MQ=255)
                                                                                   gTACTGGTTCGGGCGGCCAGTGCCGGCTGAGGATCTTCGCTGGGATTGAGACCGGGCTGCAAAGCAGCCCTCAGATTTCAACGCTTGTTCAAACCCTTCg                                                                                                              <  2:285426/100‑1 (MQ=255)
                                                                                      cTGGTTCGGGCGGCCAGTGCCGGCTGAGGATCTTCGCTGGGATTGAGACCGGGCTGCAAAGCAGCCCTCAGATTTCAACGCTTGTTCAAACCCTTCGCCAGGCGGTCCCCGCCCAGCTGGATCAACGCCACTAACGCTACCAGCATgg                                                           <  2:155948/148‑1 (MQ=255)
                                                                                       tGGTTCGGGCGGCCAGTGCCGGCTGAGGATCTTCGCTGGGATTGAGGACCGGGCTGCAAAGCAGCCCTCAGATTTCAACGCTTGTTCAAACCCTTCGCCAGGCGGTCCCCGCCCAGCTGGATCAACGCCACTAACGCTACCAGCATGGc                                                          >  1:297195/1‑149 (MQ=255)
                                                                                       tGGTTCGGGCGGCCAGTGCCGGCTGAGGATCTTCGCTGGGATTGAGACCGGGCTGCAAAGCAGCCCTCAGATTTCAACGCTTGTTCAAACCCTTCGCCAGGCGGTCCCCGCCCAGCTGGATCAACGCCACTAACGCTa                                                                    >  6:331451/1‑138 (MQ=255)
                                                                                       tGGTTCGGGCGGCCAGTGCCGGCTGAGGATCTTCGCTGGGATTGAGACCGGGCTGCAAAGCAGCCCTCAGATTTCAACGCTTGTTCAAACCCTTCGCCAGGCGGTCCCCGCCCAGCTGGATCAACGCCACTAACGCTa                                                                    <  5:331451/138‑1 (MQ=255)
                                                                                       tGGTTCGGGCGGCCAGTGCCGGCTGAGGATCTTCGCTGGGATTGAGACCGGGCTGCAAAGCAGCCCTCAGATTTCAACGCTTGTTCAAACCCTTCGCCAGGCGGTCCCCGCCCAGCTGGATCAACGCCACTAACGCTACCAGCATGGc                                                          >  4:348637/1‑148 (MQ=255)
                                                                                             gggCGGCCAGTGCCGGCTGAGGATCTTCGCTGGGATTGAGACCGGGCTGCAAAGCAGCCCTCAGATTTCAACGCTTGTTCAAACCCTTCGCCAGGCGGTCCCCGCCCAGCTGGATCAACGCCACTAACGCTACCAGCATGGCGATCACg                                                   >  1:347631/1‑149 (MQ=255)
                                                                                               gcggcCAGTGCCGGCTGAGGATCTTCGCTGGGATTGAGACCGGGCTGCAAAGCAGCCCTCAGATTTCAACGCTTGTTCAAACCCTTCGCCAGGCGGTCCCCGCCCAGCTGGATCAACGCCACTAACGCTACCAGCATgg                                                           <  1:258408/139‑1 (MQ=255)
                                                                                               gcggcCAGTGCCGGCTGAGGATCTTCGCTGGGATTGAGACCGGGCTGCAAAGCAGCCCTCAGATTTCAACGCTTGTTCAAACCCTTCGCCAGGCGGTCCCCGCCCAGCTGGATCAACGCCACTAACGCTACCAGCATgg                                                           >  2:258408/1‑139 (MQ=255)
                                                                                                   ccAGTGCCGGCTGAGGATCTTCGCTGGGATTGAGACCGGGCTGCAAAGCAGCCCTCAGATTTCAACGCTTGTTCAAACCCTTCGCCAGGCGGTCCCCGCCCAGCTGGATCAACGCCACTAACGCTACCAGCATGGCGATCACGGTCAGc                                             <  5:397358/149‑1 (MQ=255)
                                                                                                        gCCGGCTGAGGATCTTCGCTGGGATTGAGACCGGGCTGCAAAGCAGCCCTCAGATTTCAACGCTTGTTCAAACCCTTCGCCAGGCGGTCCCCGCCCAGCTGGATCAACGCCACTAACGCTACCAGCATGGCGATCACGGTCAGCATGa                                         >  4:87788/1‑148 (MQ=255)
                                                                                                                  gATCTTCGCTGGGATTGAGACCGGGCTGCAAAGCAGCCCTCAGATTTCAACGCTTGTTCAAACCCTTCGCGAGGCGGTCCCCGCCCAGCTGGATCAACGCCACTAACGCTACCAGCATGGCGATCACGGTCAGCATGATCTGGCTGTc                               <  4:305113/148‑1 (MQ=255)
                                                                                                                        cGCTGGGATTGAGACCGGGCTGCAAAGCAGCCCTCAGATTTCAACGCTTGTTCAAACCCTTCGCCAGGCGGTCCCCGCCCAGCTGGATCAACGCCACTAACg                                                                       >  6:139572/1‑102 (MQ=255)
                                                                                                                        cGCTGGGATTGAGACCGGGCTGCAAAGCAGCCCTCAGATTTCAACGCTTGTTCAAACCCTTCGCCAGGCGGTCCCCGCCCAGCTGGATCAACGCCACTAACg                                                                       <  5:139572/102‑1 (MQ=255)
                                                                                                                        cGCTGGGATTGAGACCGGGCTGCAAAGCAGCCCTCAGATTTCAACGCTTGTTCAAACCCTTCGCCAGGCGGTCCCCGCCCAGCTGGATCAACGCCACTAACGCTACCAGCATGGCGATCACGGTCAGCATGATCTGGCTGTCGAAAc                          <  2:347631/147‑1 (MQ=255)
                                                                                                                             ggATTGAGACCGGGCTGCAAAGCAGCCCTCAGATTTCAACGCTTGTTCAAACCCTTCGCCAGGCGGTCCCCGCCCAGCTGGATCAACGCCACTAACGCTACCAGCATGGCGATCACGGTCAGCATGATCTGGCTGTCGAAACGCTGGTa                   <  3:348637/149‑1 (MQ=255)
                                                                                                                               aTTGAGACCGGGCTGCAAAGCAGCCCTCAGATTTCAACGCTTGTTCAAACCCTTCGCCAGGCGGTCCCCGCCCAGCTGGATCAACGCCACTAACGCTACCAGCATGGCGATCACGGTCAGCATGATCTGGCTGTTGAAACGCTGGTAGc                 >  8:258433/1‑149 (MQ=255)
                                                                                                                               aTTGAGACCGGGCTGCAAAGCAGCCCTCAGATTTCAACGCTTGTTCAAACCCTTCGCCAGGCGGTCCCCGCCCAGCTGGATCAACGCCACTAACGCTACCAGCATGGCGATCACGGTCAGCATGATCTGGCTGTCGAAACGCTGGTAGc                 <  3:87788/149‑1 (MQ=255)
                                                                                                                                  gagaCCGGGCTGCAAAGCAGCCCTCAGATTTCAACGCTTGTTCAAACCCTTCGCCAGGCGGTCCCCGCCCAGCTGGATCAACGCCACTAACGCTACCAGCATGGCGATCACGGTCAGCATGATCTGGCTGTCGAAACGCTGGTAGCCa               <  2:25217/148‑1 (MQ=255)
                                                                                                                                   agaCCGGGCTGCAAAGCAGCCCTCAGATTTCAACGCTTGTTCAAACCCTTCGCCAGGCGGTCCCCGCCCAGCTGGATCAACGCCACTAACGCTACCAGCATGGCGATCACGGTCAGCATGATCTGGCTGTCGAAACGCTGGTAGCCAt              <  6:405055/148‑1 (MQ=255)
                                                                                                                                      tgagggCTGCTTTGCAGCCCTCAGATTTCAACGCTTGTTCAAACCCTTCGCCAGGCGGTCCCCGCCCAGCTGGATCAACGCCACTAACGCTACCAGCATGGCGATCACGGTCa                                               <  3:25414/110‑1 (MQ=255)
                                                                                                                                      tgagggCTGCTTTGCAGCCCTCAGATTTCAACGCTTGTTCAAACCCTTCGCCAGGCGGTCCCCGCCCAGCTGGATCAACGCCACTAACGCTACCAGCATGGCGATCACGGTCa                                               >  4:25414/4‑113 (MQ=255)
                                                                                                                                      tgagggCTGCTTTGCAGCCCTCAGATTTCAACGCTTGTTCAAACCCTTCGCCAGGCGGTCCCCGCCCAGCTGGATCAACGCCACTAACGCTACCAGCATGGCGATCACGGTCAGCATGATCTGGCTGTCGAAACGCTGGTAGCCATAcc           >  3:26389/4‑149 (MQ=255)
                                                                                                                                        agggCTGCTTTGCAGCCCTCAGATTTCAACGCTTGTTCAAACCCTTCGCCAGGCGGTCCCCGCCCAGCTGGATCAACGCCACTAACGCTACCAGCATGGCGATCACGGTCAGCATGATCTGGCTGTCg                              <  5:234245/127‑1 (MQ=255)
                                                                                                                                        agggCTGCTTTGCAGCCCTCAGATTTCAACGCTTGTTCAAACCCTTCGCCAGGCGGTCCCCGCCCAGCTGGATCAACGCCACTAACGCTACCAGCATGGCGATCACGGTCAGCATGATCTGGCTGTCg                              >  6:234245/2‑128 (MQ=255)
                                                                                                                                            cTGCAAAGCAGCCCTCAGATTTCAACGCTTGTTCAAACCCTTCGCCAGGCGGTCCCCGCCCAGCTGGATCAACGCCACTAACGCTACCAGCATGGCGATCACGGTCAGCATGATCTGGCTGTTGAACCGCTGGTAGCCATACCGATAgg     <  7:258433/149‑1 (MQ=255)
                                                                                                                                                aaaGCAGCCCTCAGATTTCAACGCTTGTTCAAACCCTTCGCCAGGCGGTCCCCGCCCAGCTGGATCAACGCCACTAACGCTACCAGCATGGCGATCACGGTCAGCATGATCTGGCTGTCGAAACGCTGGTAGCCATACCGATAGGCGa  >  1:393827/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
GGATGCAGGTGCATGCCGAGGGCATCGAGCAGGCGGAACAGGCCAGCTTCCTGCTGCAGGAGGAGTGCCAGTTGGGGCAGGGGTACTGGTTCGGGCGGCCAGTGCCGGCTGAGGATCTTCGCTGGGATTGAGACCGGGCTGCAAAGCAGCCCTCAGATTTCAACGCTTGTTCAAACCCTTCGCCAGGCGGTCCCCGCCCAGCTGGATCAACGCCACTAACGCTACCAGCATGGCGATCACGGTCAGCATGATCTGGCTGTCGAAACGCTGGTAGCCATACCGATAGGCGA  >  NC_002947/271865‑272154

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: