Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,149,454 C→G 9.2% V356L (GTT→CTT)  rne ← ribonuclease E

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,149,4540CG 9.2% 195.8 / 6.7 65V356L (GTT→CTT) rneribonuclease E
Reads supporting (aligned to +/- strand):  ref base C (29/30);  new base G (3/3);  total (32/33)
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.25e-01

GGCAGACGATGCCGCTGCTTTCGCCCAGCGATGGGCGCAGACGCTGACGGGACATTTCCAGCAGGCCGAAGCGCGAGATGCGGCCAACCTGCACGCGGGCGCGGTCGGCTTCCAGGCATTCACGCACGCGCTCTTCAACAGCGCGCTGGTTCTTCGCCGGGGTCATGTCGATGAAGTCGATGACG‑ATCAGGCCGCCGATGTCACGCAGGCGCAGCTGGCGGGCGATTTCCTCGGCCGCTTCCAGGTTGGTCTGCAGGGCGGTTTCTTCGATGTCGCTGCCTTTGG  >  NC_002947/2149316‑2149600
                                                                                                                                          |                                                                                                                                                   
ggCAGACGATGCCGCTGCTTTCGCCCAGCGATGGGCGCAGACGCTGACGGGACATTTCCAGCAGGCCGAAGCGCGAGATGCGGCCAACCTGCACGCGGGCGCGGTCGGCTTCCAGGCATTCACGCACGCGCTCTTCAACAGCGCGCTgg                                                                                                                                           <  3:118399/149‑1 (MQ=255)
   aGACGATGCCGCTGCTTTCGCCCAGCGATGGGCGCAGACGCTGACGGGACATTTCCAGCAGGCCGAAGCGCGAGATGCGGCCAACCTGCACGCGGGCGCGGTCGGCTTCCAGGCATTCACGCACGCGCTCTTCAACAGCGCGCTGGttc                                                                                                                                        <  7:206642/149‑1 (MQ=255)
   aGACGATGCCGCTGCTTTCGCCCAGCGATGGGCGCAGACGCTGACGGGACATTTCCAGCAGGCCGAAGCGCGAGATGCGGCCAACCTGCACGCGGGCGCGGTCGGCTTCCAGGCATTCACGCACGCGCTCTTCAACAGCGCGCTGGttc                                                                                                                                        <  7:155359/149‑1 (MQ=255)
       gATGCCGCTGCTTTCGCCCAGCGATGGGCGCAGACGCTGACGGGACATTTCCAGCAGGCCGAAGCGCGAGATGCGGCCAACCTGCACGCGGGCGCGGTCGGCTTCCAGGCATTCACGCACGCGCTCTTCAACAGCGCGCTGGTTCTTCg                                                                                                                                    >  3:288209/1‑149 (MQ=255)
           ccGCTGCTTTCGCCCAGCGATGGGCGCAGACGCTGACGGGACATTTCCAGCAGGCCGAAGCGCGAGATGCGGCCAACCTGCACGCGGGCGCGGTCGGCTTCCAGGCATTCACGCACGCGCTGTTGAAGAgcgcgtgcgtga                                                                                                                                        >  1:271206/1‑134 (MQ=255)
           ccGCTGCTTTCGCCCAGCGATGGGCGCAGACGCTGACGGGACATTTCCAGCAGGCCGAAGCGCGAGATGCGGCCAACCTGCACGCGGGCGCGGTCGGCTTCCAGGCATTCACGCACGCGCTGTTGAAGAgcgcgtgcgtga                                                                                                                                        <  2:271206/141‑8 (MQ=255)
             gctgctTTCGCCCAGCGATGGGCGCAGACGCTGACGGGACATTTCCAGCAGGCCGAAGCGCGAGATGCGGCCAACCTGCACGCGGGCGCGGTCGGCTTCCAGGCATTCACGCACGCGCTCTTCAACAGCGCGCTGGTTCTTCGCCgggg                                                                                                                              <  4:288209/149‑1 (MQ=255)
                 tctgcGCCCATCGCTGGGCGCAGACGCTGACGGGACATTTCCAGCAGGCCGAAGCGCGAGATGCGGCCAACCTGCACGCGGGCGCGGTCGGCTTCCAGGCATTCACGCACGCGCTCTTCAACAGCGCGCTGGTTCTTCGCCGGGGTCAt                                                                                                                          >  8:196224/5‑149 (MQ=255)
                 ctTTCGCCCAGCGATGGGCGCAGACGCTGACGGGACATTTCCAGCAGGCCGAAGCGCGAGATGCGGCCAACCTGCACGCGGGCGCGGTCGGCTTCCAGGCATTCACGCACGCGCTCTTCAACAGCGCGCTGGTTCTTCGCCGGGGTc                                                                                                                            <  3:198605/147‑1 (MQ=255)
                  tttCGCCCAGCGATGGGCGCAGACGCTGACGGGACATTTCCAGCAGGCCGAAGCGCGAGATGCGGCCAACCTGCACGCGGGCGCGGTCGGCTTCCAGGCATTCACGCACGCGCTCTTCAACAGCGCGCTGGTTCTTCGCCGGGGTCATg                                                                                                                         >  4:206416/1‑149 (MQ=255)
                  tttCGCCCAGCGATGGGCGCAGACGCGGACGGGACATTTCCAGCAGGCCGAAGCGCGAGATGCGGCCAACCTGCACGCGGGCGCGGTCGGCTTCCAGGCATGCACGCACGCGCTCTTCAACAGCGCGCTGGTTCTTCGCCGGGGTCATg                                                                                                                         >  2:203959/1‑149 (MQ=255)
                    tCGCCCAGCGATGGGCGCAGACGCTGACGGGACATTTCCAGCAGGCCGAAGCGCGAGATGCGGCCAACCTGCACGCGGGCGCGGTCGGCTTCCAGGCATTCACGCACGCGCTCTTCAACAGCGCGCTGGTTCTTCGCCGGGGTCATGt                                                                                                                        <  3:206416/148‑1 (MQ=255)
                    tCGCCCAGCGATGGGCGCAGACGCTGACGGGACATTTCCAGCAGGCCGAAGCGCGAGATGCGGCCAACCTGCACGCGGGCGCGGTCGGCTTCCAGGCATTCACGCACGCGCTCTTCAACAGCGCGCTGGTTCTTCGCCGGGGTCATGt                                                                                                                        <  4:16725/148‑1 (MQ=255)
                      gCCCAGCGATGGGCGCAGACGCTGACGGGACATTTCCAGCAGGCCGAAGCGCGAGATGCGGCCAACCTGCACGCGGGGGCAATAAATATCCAAACATTAACGCACGCGCTCTTAAACAGCGCGCTGGGTCTTAGCCCGGGTCATGTCGa                                                                                                                     >  1:115158/1‑149 (MQ=255)
                      gCCCAGCGATGGGCGCAGACGCTGACGGGACATTTCCAGCAGGCCGAAGCGCGAGATGCGGCCAACCTGCACGCGGGCGCGGTCGGCTTCCAGGCATTCACGCACGCGCTCTTCAACAGCGCGCTGGTTCTTCGCCGGGGTCATGTCGa                                                                                                                     >  2:200214/1‑149 (MQ=255)
                          aGCGATGGGCGCAGACGCTGACGGGACATTTCCAGCAGGCCGAAGCGCGAGATGCGGCCAACCTGCACGCGGGCGCGGTCGGCTTCCAGGCATTCACGCACGCGCTCTTCAACAGCGCGCTGGTTCTTCGCCGGGGTCATGTCGATGaa                                                                                                                 <  4:162423/149‑1 (MQ=255)
                               tGGGCGCAGACGCTGACGGGACATTTCCAGCAGGCCGAAGCGCGAGATGCGGCCAACCTGCACGCGGGCGCGGTCGGCTTCCAGGCATTCACGCACGCGCTCTTCAACAGCGCGCTGGTTCTTCGCCGGGGTCATGTCGATGAAGTCGa                                                                                                            <  8:237479/149‑1 (MQ=255)
                                  gcgcAGACGCTGACGGGACATTTCCAGCAGGCCGAAGCGCGAGATGCGGCCAACCTGCACGCGGGCGCGGTCGGCTTCCAGGCATTCACGCACGCGCTCTTCAACAGCGCGCTGGTTCTTCGCCGGGGTCATGTCGATGAAGTCGATGa                                                                                                         >  1:3125/1‑149 (MQ=255)
                                  gcgcAGACGCTGACGGGACATTTCCAGCAGGCCGAAGCGCGAGATGCGGCCAACCTGCACGCGGGCGCGGTCGGCTTCCAGGCATTCACGCACGCGCTCTGCAACAGCGCGCTGGTTCTTCGCCGGGGTCATGTCGATGAGGTCGATGa                                                                                                         >  8:226220/1‑149 (MQ=255)
                                       gACGCTGACGGGACATATCCAGCAGGCCGAAGCGCGAGATGCGGCCAACCTGCACGCGGGCGCGGTCGGCTTCCAGGCATTCACGCACGCGCTCTTCAACAGCGCGCTGGTTCTTCGCCGGGGTCATGTCGATGAAGTCGATGACG‑ATc                                                                                                   <  2:36685/149‑1 (MQ=255)
                                         cGCTGACGGGACATTTCCAGCAGGCCGAAGCGCGAGATGCGGCCAACCTGCACGCGGGCGCGGTCGGCTTCCAGGCATTCACGCACGCGCTCTTCAACAGCGCGCTGGTTCTTCGCCGGGGTCATGTCGATGAAGTCGATGACG‑ATCAg                                                                                                 <  2:215629/149‑1 (MQ=255)
                                                   aCATTTCCAGCAGGCCGAAGCGCGAGATGCGGCCAACCTGCACGCGGGCGCGGTCGGCTTCCAGGCATTCACGCACGCGCTCTTCAACAGCGCGCTGGTTCTTCGCCGGGGTCATGTCGATGAAGTCGATGACG‑ATCAGGCCGCCGATg                                                                                       <  4:155914/149‑1 (MQ=255)
                                                        tCCAGCAGGCCGAAGCGCGAGATGCGGCCAACCTGCACGCGGGCGCGGTCGGCTTCCAGGCATTCACGCACGCGCTCTTCAACAGCGCGCTGGTTCTTCGCCGGGGTCATGTCGATGAAGTCGATGACG‑ATCAGGCCGCCGATGTCAc                                                                                   <  1:200214/148‑1 (MQ=255)
                                                                   gAAGCGCGAGATGCGGCCAACCTGCACGCGGGCGCGGTCGGCTTCCAGGCATTCACGCACGCGCTCTTCAACAGCGCGCTGGTTCTTCGCCGGGGTCATGTCGATGAAGTCGATGACG‑ATCAGGCCGCCGATGTCACGCAGGCGCAGCt                                                                       >  6:71108/1‑149 (MQ=255)
                                                                   gAAGCGCGAGATGCGGCCAACCTGCACGCGGGCGCGGTCGGCTTCCAGGCATTCACGCACGCGCTCTTCAACAGCGCGCTGGTTCTTCGCCGGGGTCATGTCGATGAAGTCGATGACG‑ATCAGGCCGCCGATGTCACGCAGGCGCAGCt                                                                       >  4:4515/1‑149 (MQ=255)
                                                                        gcgAGATGCGGCCAACCTGCACGCGGGCGCGGTCGGCTTCCAGGCATTCACGCACGCGCTCTTCAACAGCGCGCTGGTTCTTCGCCGGGGTCATGTCGATGAAGTCGATGACG‑ATCAGGCCGCCGATGTCACGCAGGCGCAGCTggcgg                                                                  <  7:196224/149‑1 (MQ=255)
                                                                            gaTGCGGCCAACCTGCACGCGGGCGCGGTCGGCTTCCAGGCATTCACGCACGCGCTCTTCAACAGCGCGCTGGTTCTTCGCCGGGGTCATGTCGATGAAGTCGATGACG‑ATCAGGCCGCCGATGTCACGCAGGCGCAGCTGGCGGGCGa                                                              >  3:187948/1‑149 (MQ=255)
                                                                              tGCGGCCAACCTGCACGCGGGCGCGGTCGGCTTCCAGGCATTCACGCACGCGCTCTTCAACAGCGCGCTGGTTCTTCGCCGGGGTCATGTCGATGAAGTCGATGACG‑ATCAGGCCGCCGATGTCACGCAGGCGCAGCTGGCGGGCGAtt                                                            <  5:71108/149‑1 (MQ=255)
                                                                                 ggCCAACCTGCACGCGGGCGCGGTCGGCTTCCAGGCATTCACGCACGCGCTCTTCAACAGCGCGCTGGTTCTTCGCCGGGGTCATGTCGATGAAGTCGATGACG‑ATCAGgccgcc                                                                                           <  7:157456/115‑1 (MQ=255)
                                                                                 ggCCAACCTGCACGCGGGCGCGGTCGGCTTCCAGGCATTCACGCACGCGCTCTTCAACAGCGCGCTGGTTCTTCGCCGGGGTCATGTCGATGAAGTCGATGACG‑ATCAGgccgcc                                                                                           >  8:157456/1‑115 (MQ=255)
                                                                                   ccAACCTGCACGCGGGCGCGGTCGGCTTCCAGTCATACACGCACGCGCTCTTCAACAGCGAGCTGGTTCTTCGAAGGCGTCATGAGGATGAAGTCGATGAAT‑ATCAGGACGCCGATGTCACGCAGGCGCAGCTGGCGGGCGATATCCTc                                                       >  4:211391/1‑149 (MQ=255)
                                                                                    cAACCTGCACGCGGGCGCGGTCGGCTTCCAGGCATTCACGCACGCGCTCTTCAACAGCGCGCTGGTTCTTCGCCGGGGTCATGTCGATGAAGTCGATGACG‑ATCAGGCCGCCGATGTCACGCAGGCGCAGCTGGCGGGCGATTTCCTCg                                                      <  2:3125/149‑1 (MQ=255)
                                                                                      aCCTGCACGCGGGCGCGGTCGGCTTCCAGGCATTCACGCACGCGCTCTTCAACAGCGCGCTGGTTCTTCGCCGGGGTCATGTCGATGAAGTCGATGACG‑ATCAGGCCGCCGATGTCACGCAGGCGCAGCTGGCGGGCGATTTCCTCGGc                                                    >  1:58323/1‑149 (MQ=255)
                                                                                         tGCACGCGGGCGCGGTCGGCTTCCAGGCATTCACGCACGCGCTCTTCAACAGCGCGCTGGTTCTTCGCCGGGGTCATGTCGATGAAGTCGATGACG‑ATCAGGCCGCCGATGTCACGCAGGCGCAGCTGGCGGGCGATTTCCTCg                                                      <  3:235362/144‑1 (MQ=255)
                                                                                         tGCACGCGGGCGCGGTCGGCTTCCAGGCATTCACGCACGCGCTCTTCAACAGCGCGCTGGTTCTTCGCCGGGGTCATGTCGATGAAGTCGATGACG‑ATCAGGCCGCCGATGTCACGCAGGCGCAGCTGGCGGGCGATTTCCTCg                                                      >  4:235362/1‑144 (MQ=255)
                                                                                          gCACGCGGGCGCGGTCGGCTTCCAGGCATTCACGCACGCGCTCTTCAACAGCGCGCTGGTTCTTCGCCGGGGTCATGTCGATGAAGTCGATGACG‑ATCAGGCCGCCGATGTCACGCAGGCGCAGCTGGCGGGCGATTTCCTCGGCCGCt                                                >  4:114695/1‑149 (MQ=255)
                                                                                                 ggCGCGGTCGGCTTCCAGGCATTCACGCACGCGCTCTTCAACAGCGCGCTGGTTCTTCGCCGGGGTCATGTCGATGAAGTCGATGACG‑ATCAGGCCGCCGATGTCACGCAGGCGCAGCTGGCGGGCGATTTCCTCGGCCGCTTCCAGGt                                         <  2:115158/149‑1 (MQ=255)
                                                                                                   cgcgGTCGGCTTCCAGGCATTCACGCACGCGCTCTTCAACAGCGCGCTGGTTCTTCGCCGGGGTCATGTCGATGAAGTCGATGACG‑ATCAGGCCGCCGATGTCACGCAGGCGCAGCTGGCGGGCGATTTCCTCGGCCGCTTCCAGGTTg                                       >  7:20791/1‑149 (MQ=255)
                                                                                                   cgcgGTCGGCTTCCAGGCATTCACGCACGCGCTCTTCAACAGCGCGCTGGTTCTTCGCCGGGGTCATGTCGATGAAGTCGATGACG‑ATCAGGCCGCCGATGTCACGCAGGCGCAGCTGGCGGGCGATTTCCTCGGCCGCTTCCAGGTTg                                       <  3:211391/149‑1 (MQ=255)
                                                                                                       gTCGGCTTCCAGGCATTCACGCACGCGCTCTTCAACAGCGCGCTGGTTCTTCGCCGGGGTCATGTCGATGAAGTCGATGACG‑ATCAGGCCGCCGATGTCACGCAGGCGCAGCTGGCGGGCGATTTCCTCGGCCGCTTCCAGGTTGGTCt                                   >  2:98112/1‑149 (MQ=255)
                                                                                                        tCGGCTTCCAGGCATTCACGCACGCGCTCTTCAACAGCGCGCTGGTTCTTCGCCGGGGTCATGTCGATGAAGTCGATGACG‑ATCAGGCCGCCGATGTCACGCAGGCGCAGCTGGCGGGCGATTTCCTCGGCCGCTTCCAGGTTGGTCTg                                  <  4:187948/149‑1 (MQ=255)
                                                                                                          ggCTTCCAGGCATTCACGCACGCGCTCTTCAACAGCGCGCTGGTTCTTCGCCGGGGTCATGTCGATGAAGTCGATGACG‑ATCAGGCCGCCGATGTCACGCAGGCGCAGCTGGCGGGCGATTTCCTCGGCCGCTTCCAGGTTGGTCTg                                  >  2:10350/1‑147 (MQ=255)
                                                                                                             ttCCAGGCATTCACGCACGCGCTCTTCAACAGCGCGCTGGTTCTTCGCCGGGGTCATGTCGATGAAGTCGATGACG‑ATCAGGCCGCCGATGTCACGCAGGCGCAGCTGGCGGGCGATTTCCTCGGCCGCTTCCAGGTTGGTCTGCAggg                             >  7:163113/1‑149 (MQ=255)
                                                                                                                 aGGCATTCACGCACGCGCTCTTCAACAGCGCGCTGGTTCTTCGCCGGGGTCATGTCGATGAAGTCGATGACG‑ATCAGGCCGCCGATGTCACGCAGGCGCAGCTGGCGGGCGATTTCCTCGGCCGCTTCCAGGTTGGTCTGCAGGGCGGt                         <  3:114695/149‑1 (MQ=255)
                                                                                                                 aGGCATTCACGCACGCGCTCTTCAACAGCGCGCTGGTTCTTCGCCGGGGTCATGTCGATGAAGTCGATGACG‑ATCAGGCCGCCGATGTCACGCAGGCGCAGCTGGCGGGCGATTTCCTCGGCCGCTTCCAGGTTGGTCTGCAGGGCGGt                         >  7:161922/1‑149 (MQ=255)
                                                                                                                      ttCACGCACGCGCTCTTCAACAGCGCGCTGGTTCTTCGCCGGGGTCATGTCGATGAAGTCGATGACG‑ATCAGGCCGCCGATGTCACGCAGGCGCAGCTGGCGGGCGATTTCCTCGGCCGCTTCCAGGTTGGTCTGCAGGGCGGTttctt                    <  8:163113/149‑1 (MQ=255)
                                                                                                                          cgcacgcGCTCTTCAACAGCGCGCTGGTTCTTCGCCGGGGTCATGTCGATGAAGTCGATGACG‑ATCAGGCCGCCGATGTCACGCAGGCGCAGCTGGCGGGCGATTTCCTCGGCCGCTTCCAGGTTGGTCTGCAGGGCGGTTTCTTCGAt                >  8:58116/1‑149 (MQ=255)
                                                                                                                          ccagcgcgCTGTTGAAGAGCGCGCTGGTTCTTCGCCGGGGTCATGTCGATGAAGTCGATGACG‑ATCAGGCCGCCGATGTCACGCAGGCGCAGCTGGCGGGCGATTTCCTCGGCCGCTTCCAGGTTg                                       >  7:192564/5‑126 (MQ=38)
                                                                                                                          ccagcgcgCTGTTGAAGAGCGCGCTGGTTCTTCGCCGGGGTCATGTCGATGAAGTCGATGACG‑ATCAGGCCGCCGATGTCACGCAGGCGCAGCTGGCGGGCGATTTCCTCGGCCGCTTCCAGGTTg                                       <  8:192564/122‑1 (MQ=38)
                                                                                                                            agcgcgCTGTTGAAGAGCGCGCTGGTTCTTCGCCGGGGTCATGTCGATGAAGTCGATGACG‑ATCAGGCCGCCGATGTCACGCAGGCGCAGCTGGCGGGCGATTTCCTCGGCCGCTTCCAGGTTGGTCTGc                                 <  3:261041/128‑1 (MQ=255)
                                                                                                                            agcgcgCTGTTGAAGAGCGCGCTGGTTCTTCGCCGGGGTCATGTCGATGAAGTCGATGACG‑ATCAGGCCGCCGATGTCACGCAGGCGCAGCTGGCGGGCGATTTCCTCGGCCGCTTCCAGGTTGGTCTGc                                 >  4:261041/3‑130 (MQ=255)
                                                                                                                             acgcGCTCTTCAACAGCGCGCTGGTTCTTCGCCGGGGTCATGTCGATGAAGTCGATGACG‑ATCAGGCCGCCGATGTCACGCAGGCGCAGCTggcgg                                                                  >  2:86269/1‑96 (MQ=255)
                                                                                                                             acgcGCTCTTCAACAGCGCGCTGGTTCTTCGCCGGGGTCATGTCGATGAAGTCGATGACG‑ATCAGGCCGCCGATGTCACGCAGGCGCAGCTggcgg                                                                  >  2:81426/1‑96 (MQ=255)
                                                                                                                             acgcGCTCTTCAACAGCGCGCTGGTTCTTCGCCGGGGTCATGTCGATGAAGTCGATGACG‑ATCAGGCCGCCGATGTCACGCAGGCGCAGCTggcgg                                                                  <  1:86269/96‑1 (MQ=255)
                                                                                                                             acgcGCTCTTCAACAGCGCGCTGGTTCTTCGCCGGGGTCATGTCGATGAAGTCGATGACG‑ATCAGGCCGCCGATGTCACGCAGGCGCAGCTggcgg                                                                  <  1:81426/96‑1 (MQ=255)
                                                                                                                             acgcGCTCTTCAACAGCGCGCTGGTTCTTCGCCGGGGTCATGTCGATGAAGTCGATGACG‑ATCAGGCCGCCGATGTCACGCAGGCGCAGCTGGCGGGCGATTTCCTCGGCCGCTTCCAGGTTGGTCTGCAGGGCGGTTTCTTCGATGTc             >  4:201880/1‑149 (MQ=255)
                                                                                                                             acgcGCTCTTCAACAGCGCGCTGGTTCTTCGCCGGGGTCATGTCGATGAAGTCGATGACG‑ATCAGGCCGCCGATGTCACGCAGGCGCAGCTGGCGGGCGATTTCCTCGGCCGCTTCCAGGTTGGTCTGCAGGGCGGTTTCTTCGATGTc             <  3:4515/149‑1 (MQ=255)
                                                                                                                                      tCAACAGCGCGCTGGTTCTTCGCCGGGGTCATGTCGATGAAGTCGATGACG‑ATCAGGCCGCCGATGTCACGCAGGCGCAGCTGGCGGGCGATTTCCTCGGCCGCTTCCAGGTTGGTCTGCAGGGCGGTttcttc                   <  2:322336/134‑1 (MQ=255)
                                                                                                                                      tCAACAGCGCGCTGGTTCTTCGCCGGGGTCATGTCGATGAAGTCGATGACG‑ATCAGGCCGCCGATGTCACGCAGGCGCAGCTGGCGGGCGATTTCCTCGGCCGCTTCCAGGTTGGTCTGCAGGGCGGTttcttc                   >  1:322336/1‑134 (MQ=255)
                                                                                                                                       cAACAGCGCGCTGGTTCTTCGCCGGGGTCATGTCGATGAAGTCGATGCCG‑ATGAGTCCGCGGATGTCTCGCAGGCGCAGCTGGCGGGCGATTTCCTCGGCCGCTTCCAGGTTGGTCTGCTGGGCGGTTTCTTCGATGTCGCTGCCttt    <  1:203959/148‑1 (MQ=255)
                                                                                                                                       cAACAGCGCGCTGGTTCTTCGCCGGGGTCATGTCGATGAAGTCGATGACG‑ATCAGGCCGCCGATGTCACGCAGGCGCAGCTGGCGGGCGATTTCCTCGGCCGCTTCCAGGTTGGTCTGc                                 <  4:17463/119‑1 (MQ=255)
                                                                                                                                       cAACAGCGCGCTGGTTCTTCGCCGGGGTCATGTCGATGAAGTCGATGACG‑ATCAGGCCGCCGATGTCACGCAGGCGCAGCTGGCGGGCGATTTCCTCGGCCGCTTCCAGGTTGGTCTGc                                 >  3:17463/1‑119 (MQ=255)
                                                                                                                                         aCAGCGCGCTGGTTCTTCGCCGGGGTCATGTCGATGAAGTCGATGACGAATCAGGCCGCAGATGTCACGCAGGCGCAGCTGGCGGGCGATTTCCTCGGCCGCTTCCAGGTTGGTCTGCAGGGCGGTTTCTTCGATGTCGCTGCCTTtgg  >  1:85207/1‑149 (MQ=255)
                                                                                                                                         aCAGCGCGCTGGTTCTTCGCCGGGGTCATGTCGATGAAGTCGATGACG‑ATCAGGCCGCCGATGTCACGCAGGCGCAGCTGGCGGGCGATTTCCTCGGCCGCTTCCAGGTTGGTCTg                                  <  2:187940/116‑1 (MQ=255)
                                                                                                                                         aCAGCGCGCTGGTTCTTCGCCGGGGTCATGTCGATGAAGTCGATGACG‑ATCAGGCCGCCGATGTCACGCAGGCGCAGCTGGCGGGCGATTTCCTCGGCCGCTTCCAGGTTGGTCTg                                  >  1:187940/1‑116 (MQ=255)
                                                                                                                                          |                                                                                                                                                   
GGCAGACGATGCCGCTGCTTTCGCCCAGCGATGGGCGCAGACGCTGACGGGACATTTCCAGCAGGCCGAAGCGCGAGATGCGGCCAACCTGCACGCGGGCGCGGTCGGCTTCCAGGCATTCACGCACGCGCTCTTCAACAGCGCGCTGGTTCTTCGCCGGGGTCATGTCGATGAAGTCGATGACG‑ATCAGGCCGCCGATGTCACGCAGGCGCAGCTGGCGGGCGATTTCCTCGGCCGCTTCCAGGTTGGTCTGCAGGGCGGTTTCTTCGATGTCGCTGCCTTTGG  >  NC_002947/2149316‑2149600

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