Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,078,877 G→C 7.9% intergenic (‑122/+275) mgtE ← / ← PP_t56 magnesium transporter/tRNA‑Arg

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,078,8770GC 7.9% 190.6 / 5.6 63intergenic (‑122/+275)mgtE/PP_t56magnesium transporter/tRNA‑Arg
Reads supporting (aligned to +/- strand):  ref base G (32/26);  new base C (2/3);  total (34/29)
Fisher's exact test for biased strand distribution p-value = 6.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

TTTTAACTTCTACTTCGGTCATAGGCGAACTCCACTCCCAGCAGCGGAGCACGCCGAGGGAATCAATCGGTCAGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAATTTTTGGCCAGAACAATCGTTTGCGCGATAAAGCTGGGACAACGCGTGAAGCTATCAGTGCGACGATTGTTGTGGATGGGAAAGCACATTTGATTTT  >  NC_002947/5078734‑5079025
                                                                                                                                               |                                                                                                                                                    
ttttAACTTCTACTTCGGTCATAGGCGAACTCCACTCCCAGCAGCGGAGCACGCCGAGGGAATCAATCGGTCAGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCggg                                                                                                                                                  <  5:139637/148‑1 (MQ=255)
     aCTTCTACTTCGGTCATAGGCGAACTCCACTCCCAGCAGCGGAGCACGCCGAGGGAATCAATCGGTCAGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGCCGGGGCttgt                                                                                                                                            <  2:9207/149‑4 (MQ=255)
       ttCTACTTCGGTCATAGGCGAACTCCACTCCCAGCAGCGGAGCACGCCGAGGGAATCAATCGGTCAGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCcggcgg                                                                                                                                                   >  6:25679/1‑140 (MQ=255)
       ttCTACTTCGGTCATAGGCGAACTCCACTCCCAGCAGCGGAGCACGCCGAGGGAATCAATCGGTCAGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCcggcgg                                                                                                                                                   <  5:25679/140‑1 (MQ=255)
         cTACTTCGGTCATAGGCGAACTCCACTCCCAGCAGCGGAGCACGCCGAGGGAATCAATCGGTCAGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGgc                                                                                                                                        >  7:246345/1‑149 (MQ=255)
            cTTCGGTCATAGGCGAACTCCACTCCCAGCAGCGGAGCACGCCGAGGGAATCAATCGGTCAGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGgcg                                                                                                                                       <  8:246345/147‑1 (MQ=255)
             ttCGGTCATAGGCGAACTCCACTCCCAGCAGCGGAGCACGCCGAGGGAATCAATCGGTCAGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCaa                                                                                                                                    >  4:50240/1‑149 (MQ=255)
                    aTAGGCGAACTCCACTCCCAGCAGCGGAGCACGCCGAGGGAATCAATCGGTCAGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATa                                                                                                                               <  4:345588/147‑1 (MQ=255)
                    aTAGGCGAACTCCACTCCCAGCAGCGGAGCACGCCGAGGGAATCAATCGGTCAGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATa                                                                                                                               >  3:345588/1‑147 (MQ=255)
                      aGGCGAACTCCACTCCCAGCAGCGGAGCACGCCGAGGGAATCAATCGGTCAGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGCCGGGGCttgtg                                                                                                                                           >  1:241818/1‑129 (MQ=255)
                      aGGCGAACTCCACTCCCAGCAGCGGAGCACGCCGAGGGAATCAATCGGTCAGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGCCGGGGCttgtg                                                                                                                                           <  4:295225/133‑5 (MQ=255)
                      aGGCGAACTCCACTCCCAGCAGCGGAGCACGCCGAGGGAATCAATCGGTCAGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGCCGGGGCttgtg                                                                                                                                           >  3:295225/1‑129 (MQ=255)
                      aGGCGAACTCCACTCCCAGCAGCGGAGCACGCCGAGGGAATCAATCGGTCAGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGCCGGGGCttgtg                                                                                                                                           <  2:241818/133‑5 (MQ=255)
                              tccactccCAGCAGCGGAGCACGCCGAGGGAATCAATCGGTCAGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATacac                                                                                                                            <  1:255962/140‑1 (MQ=255)
                              tccactccCAGCAGCGGAGCACGCCGAGGGAATCAATCGGTCAGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATacac                                                                                                                            >  2:255962/1‑140 (MQ=255)
                                 actccCAGCAGCGGAGCACGCCGAGGGAATCAATCGGTCAGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGt                                                                                                                >  8:98623/1‑149 (MQ=255)
                                   tccCAGCAGCGGAGCACGCCGAGGGAATCAATCGGTCAGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTc                                                                                                               >  6:205316/1‑148 (MQ=255)
                                     cctgcagcGGAGGACGCCGAGGGAATCATTCTGTCAGATCTTTCTGTATGATTCGTGTTAGCGATTAACTTCTGGGTAAGTCCATGGTGGTTTTCCACAAGCGCCTGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCa                                                                                                              <  4:14639/144‑1 (MQ=255)
                                     ccAGCAGCGGAGCACGCCGAGGGAATCAATCGGTCAGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTc                                                                                                               <  5:205316/146‑1 (MQ=255)
                                     ccAGCAGCGGAGCACGCCGAGGGAATCAATCGGTCAGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCACGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGt                                                                                                                >  3:14639/1‑145 (MQ=255)
                                      cagcagCGGAGCACGCCGAGGGAATCAATCGGTCAGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGt                                                                                                                >  4:94854/1‑144 (MQ=255)
                                                    gCCGAGGGAATCAATCGGTCAGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCgggg                                                                                                                                                 >  1:102649/1‑97 (MQ=255)
                                                    gCCGAGGGAATCAATCGGTCAGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCgggg                                                                                                                                                 <  2:102649/97‑1 (MQ=255)
                                                             atcaatcGGTCAGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAAttttt                                                                                              <  5:217930/139‑1 (MQ=255)
                                                             atcaatcGGTCAGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAAttttt                                                                                              >  6:217930/1‑139 (MQ=255)
                                                                aatcGGTCAGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAATTTTTGGCCAGAACAATc                                                                                 >  4:168713/1‑149 (MQ=255)
                                                                aatcGGTCAGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAATTTTTGGCCAGAACAATc                                                                                 >  4:279945/1‑149 (MQ=255)
                                                                    ggTCAGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTTGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAATTTTTGGCCAGAACAATCGttt                                                                             <  7:217539/149‑1 (MQ=255)
                                                                     gTCAGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAATTTTTGGCCAGAACAATCGTTTg                                                                            >  3:68143/1‑149 (MQ=255)
                                                                     gTCAGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTACATGGTGGTTTTCCACAAGCCCCCGCGGGGCTGACGGGCGCAATCATACACAGCTAAAGCTGTCAGATCGCTTAAATTTTTTGCCAGAACaa                                                                                   >  1:9996/1‑142 (MQ=255)
                                                                     gTCAGATATTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGTTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAATTTTTGGCCAGAACaa                                                                                   <  2:9996/142‑1 (MQ=255)
                                                                       cAGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAATTTTTGGCCAGAACAATCGTTTgcg                                                                          >  6:98967/1‑149 (MQ=255)
                                                                          aTCTTTCTGTATGAATCGTGTTAGCGAGTATCTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATTGCTTAAATTTTTGGCCAGAACAATCGTTTGCGCGa                                                                       >  4:158947/1‑149 (MQ=255)
                                                                             tttCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAATTTTTGGCCAGAACAATCGttt                                                                             <  7:303913/140‑1 (MQ=255)
                                                                             tttCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAATTTTTGGCCAGAACAATCGttt                                                                             >  8:303913/1‑140 (MQ=255)
                                                                                 tGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAATTTTTGGCCAGAACAATCGTTTGCGCGATAAAGCt                                                                >  5:174157/1‑149 (MQ=255)
                                                                                   tATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAATTTTTGGCCAGAACAATCGTTTGCGCGATAAAGCTgg                                                              <  7:98623/149‑1 (MQ=255)
                                                                                     tGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAATTTTTGGCCAGAACAATCGTTTGCGCGATAAAGCTggg                                                             >  3:222380/1‑148 (MQ=255)
                                                                                     tGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAATTTTTGGCCAGAACAATCGTTTGCGCGATAAAGCTggg                                                             <  4:222380/148‑1 (MQ=255)
                                                                                     tGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAATTTTTGGCCAGAACAATCGTTTGCGCGATAAAGCTGGGa                                                            <  3:94854/149‑1 (MQ=255)
                                                                                      gAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAAttttt                                                                                              <  4:38843/114‑1 (MQ=255)
                                                                                       aaTCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAATTTTTGGCCAGAACAATCGTTTGCGCGATAAAGCTGGGAc                                                           >  3:133686/1‑148 (MQ=255)
                                                                                        aTCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTAGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCAAACACTGCTAAAGCTGTAAGATCGCTTAAATTTATGGCCAGAACAATCGTTTGCGCGATAAAGCTCGGACaa                                                         >  5:333524/1‑149 (MQ=255)
                                                                                           gtgtTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAATTTTTGGCCAGAACAATCGTTTGCGCGATAAAGCTGGGACAAcgc                                                      <  6:174157/149‑1 (MQ=255)
                                                                                               tAGCTAGTAACTACTGGGGATGTCCATGGTGTTTTTCCACAAGCCCCGGCGGTGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAATTTTTGGCCAGAACATTCGTTTGCGCTATAAAGCTGGGACAACGCGTg                                                   <  6:333524/148‑1 (MQ=255)
                                                                                                 gCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAATTTTTGGCCAGAACAATCGTTTGCGCGATAAAGCTGGGACAACGCGTGaa                                                 <  3:279945/148‑1 (MQ=255)
                                                                                                 gCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAATTTTTGGCCAGAACAATCGTTTGCGCGATAAAGCTGGGACAACGCGTGAAg                                                >  3:49734/1‑149 (MQ=255)
                                                                                                   gAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAATTTTTGGCCAGAACAATCGTTTGCGCGATAAAGCTGGGACAACGCGTGAAGc                                               >  7:164221/1‑148 (MQ=255)
                                                                                                     gTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAATTTTTGGCCAGAACAATCGTTTGCGCGATAAAGCTGGGACAACGCGTGAAGCTAt                                            <  4:68143/149‑1 (MQ=255)
                                                                                                          tactGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAATTTTTGGCCAGAACAATCGTTTGCGCGATAAAGCTGGGACAACGCGTGAAGCTATCAGTg                                       >  4:338286/1‑149 (MQ=255)
                                                                                                            ctGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAATTTTTGGCCAGAACAATCGTTTGCGCGATAAAGCTGGGACAACGCGTGAAGCTATCAGTg                                       <  4:133686/147‑1 (MQ=255)
                                                                                                               ggTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAATTTTTGGCCAGAACAATCGTTTGCGCGATAAAGCTGGGACAACGCGTGAAGCTATCAGTGcgac                                   >  8:4858/1‑148 (MQ=255)
                                                                                                                     tCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAATTTTTGTACAGAACAATAGTTTGCGAGATAAAGCTGGGACAACTCGTGAAGCTATAAGTGCTACGAttttt                            >  6:240681/1‑149 (MQ=255)
                                                                                                                     tCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAATTTTTGGCCAGAACAATCGTTTGCGCGATAAAGCTg                                                               >  8:153388/1‑114 (MQ=255)
                                                                                                                     tCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAATTTTTGGCCAGAACAATCGTTTGCGCGATAAAGCTGGGACAACGCGTGAAGCTATCAGTGCGACGAttgtt                            <  5:98967/149‑1 (MQ=255)
                                                                                                                       cATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAATTTTTGGCCAGAACAATCGTTTGCGCGATAAAGCTGGGACAACGCGTGAAGCTATCAGTGCGACGATtgttgt                          >  1:91478/1‑149 (MQ=255)
                                                                                                                       cATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAATTTTTGGCCAGAACAATCGTTTGCGCGATAAAGCTGGGACAACGCGTGAAGCTATCAGTGCGACGATtgttgt                          <  7:4858/149‑1 (MQ=255)
                                                                                                                        aTGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAATTTTTGGCCAGAACAATCGTTTGCGCGATAAAGCTGGGACAACGCGTGAAGCTATCAGTGCGACGATTGTtgtg                         >  1:275853/1‑149 (MQ=255)
                                                                                                                          ggtggtTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAATTTTTGGCCAGAACAATCGTTTGCGCGATAAAGCTGGGACAAcgc                                                      <  2:271622/118‑1 (MQ=255)
                                                                                                                          ggtggtTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAATTTTTGGCCAGAACAATCGTTTGCGCGATAAAGCTGGGACAAcgc                                                      >  1:271622/1‑118 (MQ=255)
                                                                                                                            tggtTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAATTTTTGGCCAGAACAATCGTTTGCGCGATAAAGCTGGGACAACGCGTGAAGCTATCAGTGcga                                    <  3:113696/134‑1 (MQ=255)
                                                                                                                            tggtTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAATTTTTGGCCAGAACAATCGTTTGCGCGATAAAGCTGGGACAACGCGTGAAGCTATCAGTGcga                                    >  4:113696/1‑134 (MQ=255)
                                                                                                                                  tCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAATTTTTGGCCAGAACAATCGTTTGCGCGATAAAGCTGGGACAACGCGTGAAGCTATCAGTGCGACGATTGTTGTGGATGGGaaa                <  4:49734/148‑1 (MQ=255)
                                                                                                                                  tCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAATTTTTGGCCAGAACAATCGTATGCGCGATAAAGCTGGGACAACGCGTGAAGCTATCAGTGCGACGATTGTTGTGGATGGGAAAg               >  7:149169/1‑149 (MQ=255)
                                                                                                                                         gCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAATTTTTGGCCAGAACAATCGTTTGCGCGATAAAGCTGGGACAACGCGTGAAGCTATCAGTGCGACGATTGTTGTGGATGGGAAAGCACAttt        <  2:91478/149‑1 (MQ=255)
                                                                                                                                               gcggTGATGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAATTTTTGGCCAGAACAATCGTTTGCGCGATAAAGCTGGGACAACGCGTGAAGCTATCAGTGCGACGATTGTTGTGGATGGGAAAGCACATTTGAtttt  <  5:240681/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
TTTTAACTTCTACTTCGGTCATAGGCGAACTCCACTCCCAGCAGCGGAGCACGCCGAGGGAATCAATCGGTCAGATCTTTCTGTATGAATCGTGTTAGCGAGTAACTACTGGGTAAGTCCATGGTGGTTTTCCACAAGCCCCGGCGGGGCTGACGGGCGCAATCATACACTGCTAAAGCTGTCAGATCGCTTAAATTTTTGGCCAGAACAATCGTTTGCGCGATAAAGCTGGGACAACGCGTGAAGCTATCAGTGCGACGATTGTTGTGGATGGGAAAGCACATTTGATTTT  >  NC_002947/5078734‑5079025

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: