Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 609,779 A→C 9.5% L356R (CTG→CGG) ‡ PP_0525 ← B12 family TonB‑dependent receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947609,7790AC 9.5% 183.0 / 10.9 63L356R (CTG→CGG) ‡PP_0525B12 family TonB‑dependent receptor
Reads supporting (aligned to +/- strand):  ref base A (29/28);  new base C (3/3);  total (32/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.63e-01

CTTGTCGTGGCGCAGGCCAAGCTCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTGGGCCAGCTCAAACCGAGTGCTCCACTGCTC  >  NC_002947/609634‑609927
                                                                                                                                                 |                                                                                                                                                    
cTTGTCGTGGCGCAGGCCAAGCTCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAgg                                                                                                                                                    <  6:235182/148‑1 (MQ=255)
      gTGGCGCAGGCTAAGCTCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCtgttgt                                                                                                                                             <  4:96649/149‑1 (MQ=255)
      gTGGCGCAGGCCAAGCTCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCCTGCtgttgg                                                                                                                                             <  4:20942/149‑2 (MQ=255)
         gcgcAGGCCAAGCTCGGTGTAAACGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCCTGCTGTTgggta                                                                                                                                          <  4:34008/149‑2 (MQ=255)
         gcgcAGGCCAAGCTCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGACGCTACCCAACAGCAGGCTGTTGTGtt                                                                                                                                          <  2:242939/149‑1 (MQ=255)
           gcAGGCCAAGCTCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCa                                                                                                                                        <  3:116666/149‑1 (MQ=255)
                ccAAGCTCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAg                                                                                                                                   <  2:113861/149‑1 (MQ=255)
                  aaGCTCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAggg                                                                                                                                 >  1:189358/1‑149 (MQ=255)
                    gCTCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGc                                                                                                                                >  1:248303/1‑148 (MQ=255)
                          tgtAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGgttg                                                                                                                         <  2:248303/149‑1 (MQ=255)
                              aaaGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCt                                                                                                                     <  4:134556/149‑1 (MQ=255)
                               aaGTATTCGTCATGAAAGCGATGCTTGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGgtt                                                                                                                          <  6:28906/143‑1 (MQ=255)
                               aaGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGgtt                                                                                                                          >  5:28906/1‑143 (MQ=255)
                                 gTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGGCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGcc                                                                                                                  <  5:333952/149‑1 (MQ=255)
                                   aTTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGcc                                                                                                                  >  4:274062/1‑147 (MQ=255)
                                          aTGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCATGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCtt                                                                                                         >  3:32348/1‑149 (MQ=255)
                                               aGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTcg                                                                                                    >  6:101412/1‑149 (MQ=255)
                                                  gATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCgg                                                                                                 <  5:101412/149‑1 (MQ=255)
                                                   aTGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCgg                                                                                                 >  7:183824/1‑148 (MQ=255)
                                                     gCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCggtag                                                                                              >  7:42853/1‑149 (MQ=255)
                                                      cTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCggtagg                                                                                             >  2:90822/1‑149 (MQ=255)
                                                       tGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCggtaggt                                                                                            >  7:9999/1‑149 (MQ=255)
                                                        ggACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCggtaggta                                                                                           >  7:244135/1‑149 (MQ=255)
                                                         gACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTGGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTAt                                                                                          <  2:143198/149‑1 (MQ=255)
                                                              aCAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCggtag                                                                                              <  6:117103/140‑1 (MQ=255)
                                                              aCAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCggtag                                                                                              >  5:117103/1‑140 (MQ=255)
                                                                    cATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTgcg                                                                               <  8:244135/149‑1 (MQ=255)
                                                                     aTGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGAAACGGTCTTCGTACAAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTCCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTgcgc                                                                              >  4:127090/1‑149 (MQ=255)
                                                                       gATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTACCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCtt                                                                            >  1:164518/1‑149 (MQ=255)
                                                                        aTTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCtt                                                                            <  3:274062/148‑1 (MQ=255)
                                                                               cGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCt                                                                    <  8:42853/149‑1 (MQ=255)
                                                                                     tCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCtt                                                              <  1:20757/149‑1 (MQ=255)
                                                                                     tCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCtt                                                              <  2:189358/149‑1 (MQ=255)
                                                                                       ttCCGTGAAGTCCTAACTCGCATGGACACGGTCTTCGTACCAGACGCTACCCAACAGCAGGCAGTTGTGTTCACCCAGGGCTAGGTTGTTCTGACAAGTTGATTGGTCGCGGAAGGTATTGAAAGGTGCGCTTGAGTGGCAGAGCATGt                                                            >  4:100743/1‑149 (MQ=255)
                                                                                            tGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCggtag                                                                                              >  4:248004/1‑110 (MQ=255)
                                                                                            tGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCggtag                                                                                              <  3:248004/110‑1 (MQ=255)
                                                                                                gTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCttt                                                   >  7:217143/1‑149 (MQ=255)
                                                                                                 tCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTc                                                  <  6:249126/149‑1 (MQ=255)
                                                                                                      aCTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGAGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATcgc                                             >  1:237157/1‑149 (MQ=255)
                                                                                                      aCTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTCTGTTCACCCAGGCCTAGGTTCTTCTGCCAAGTTCCTTGGTCGCGCTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTACATcgc                                             >  3:94319/1‑149 (MQ=255)
                                                                                                          gCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCtt                                                              <  3:267576/128‑1 (MQ=255)
                                                                                                          gCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCtt                                                              >  4:267576/1‑128 (MQ=255)
                                                                                                                tcacGGTCTTCGTACCAGTCGCTACCCTACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTTCGCTTGAGt                                                                        <  7:205411/111‑1 (MQ=255)
                                                                                                                  acGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCtgt                                 >  2:252639/1‑149 (MQ=255)
                                                                                                                    ggTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTgg                               >  3:268658/1‑149 (MQ=255)
                                                                                                                    ggTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTgg                               <  8:183824/149‑1 (MQ=255)
                                                                                                                      tCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTGGGc                             <  8:217143/149‑1 (MQ=255)
                                                                                                                        ttCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTGGGc                             >  3:18883/1‑147 (MQ=255)
                                                                                                                               cAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTGGGCCAGCTCaaa                    >  5:319690/1‑149 (MQ=255)
                                                                                                                                aGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGGTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTGGGCCAGCTCAAAc                   <  5:122159/149‑1 (MQ=255)
                                                                                                                                 gTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTGGGCCAGCTCAAAcc                  >  5:214423/1‑149 (MQ=255)
                                                                                                                                   cGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTGCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTGGGCCAGCTCAAACCGa                <  8:9999/149‑1 (MQ=255)
                                                                                                                                    gCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTGGGCCAGCTCAAACCGAg               >  7:198823/1‑149 (MQ=255)
                                                                                                                                       acacAACAGCCTGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTg                                                                                 <  3:251769/77‑1 (MQ=255)
                                                                                                                                       acacAACAGCCTGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTg                                                                                 >  4:251769/4‑80 (MQ=255)
                                                                                                                                        cacAACAGCCTGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTGGGCCAGCTCAAACCGAGTGCt           >  3:321024/3‑149 (MQ=255)
                                                                                                                                         ccAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTGGGCCAGCTCAAACCGa                <  6:216180/143‑1 (MQ=255)
                                                                                                                                         ccAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTGGGCCAGCTCAAACCGa                >  5:216180/1‑143 (MQ=255)
                                                                                                                                         ccAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGAGGGCCAGCTCAAACCGAGAGCTc          >  8:320371/1‑149 (MQ=255)
                                                                                                                                          cAACAGCCTGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTGGGCCAGCTCAAACCGAGTGCTc          >  8:168391/1‑148 (MQ=255)
                                                                                                                                           aaCAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGt                                                                        >  3:248175/1‑85 (MQ=255)
                                                                                                                                           aaCAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGt                                                                        <  4:248175/85‑1 (MQ=255)
                                                                                                                                            aCAGGAGGCTGTTGTGTTCTCCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTGGGCCAGCTCAAACCGAGTGCTcctc       <  4:94319/149‑3 (MQ=255)
                                                                                                                                             cagcagGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTGGGCCAGCTCAAACCGAGTGCTCCACt      <  4:18883/149‑1 (MQ=255)
                                                                                                                                                cagGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTGGGCCAGCTCAAACCGAGTGCTCCACTGCt   <  8:198823/149‑1 (MQ=255)
                                                                                                                                                 agGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTGGGCCAGCTCAAACCGAGTGCTCCACTGCTc  >  5:119203/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
CTTGTCGTGGCGCAGGCCAAGCTCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTGGGCCAGCTCAAACCGAGTGCTCCACTGCTC  >  NC_002947/609634‑609927

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: