Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 307,866 +C intergenic (+881/‑988) hslO → / → PP_0254 molecular chaperone HslO/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947307,8651.C94.8% 187.8 / 3.8 58intergenic (+880/‑989)hslO/PP_0254molecular chaperone HslO/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (2/1);  new base C (26/29);  total (28/30)
Fisher's exact test for biased strand distribution p-value = 6.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

AAGCCATGTTCTCGGTGCAGAACTTCCTGCTGCCGGCTGCCGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAG‑C‑GTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATCTGGAAAGCCATCAAGCACGGCGCAG  >  NC_002947/307718‑308004
                                                                                                                                                    |                                                                                                                                            
aaGCCATGTTCTCGGTGCAGAACTTCCTGCTGCCGGCTGCCGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAG‑c                                                                                                                                             <  2:307858/149‑1 (MQ=255)
  gCCATGTTCTCGGTGCAGAACTTCCTGCTGCCGGCTGCCGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGTCAGCCGAAAG‑G‑Gt                                                                                                                                          <  6:109940/149‑1 (MQ=255)
   ccATGTTCTCGGTGCAGAACTTCCTGCTGCCGGCTGCCGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAG‑ccgt                                                                                                                                          >  8:107561/1‑146 (MQ=255)
    cATGTTCTCGGTGCAGAACTTCCTGCTGCCGGCTGCCGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAG‑ccgt                                                                                                                                          >  8:332940/1‑145 (MQ=255)
        ttCTCGGTGCAGAACTTCCTGCTGCCGGCTGCCGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCC‑GTTACCt                                                                                                                                     <  7:160425/149‑1 (MQ=255)
             ggTGCAGAACTTCCTGCTGCCGGCTGCCGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCC‑GTTACCTGATCg                                                                                                                                <  5:7518/149‑1 (MQ=255)
               tGCAGAACTTCCTGCTGCCGGCTGCCGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCC‑GTTACCTGATCGGt                                                                                                                              >  5:276703/1‑149 (MQ=255)
               tGCAGAACTTCCTGCTGCCGGCTGCCGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCC‑GTTACCTGATCGGt                                                                                                                              <  7:332940/149‑1 (MQ=255)
                   gAACTTCCTGCTGCCGGCTGCCGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCC‑GTTACCTGATCGGTgac                                                                                                                           >  4:44027/1‑148 (MQ=255)
                    aaCTTCCTGCTGCCGGCTGCCGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCC‑GTTACCTGATCGGTgacga                                                                                                                         >  6:262999/1‑149 (MQ=255)
                      cTTCCTGCTGCCGGCTGCCGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCC‑GTTACCTGATCGGTgacgac                                                                                                                        >  6:11867/1‑148 (MQ=255)
                         ccTGCTGCCGGCTGCCGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCC‑GTTACCTGATCGGTGAcgacga                                                                                                                      <  2:80349/147‑1 (MQ=255)
                         ccTGCTGCCGGCTGCCGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCC‑GTTACCTGATCGGTGAcgacga                                                                                                                      <  2:348836/147‑1 (MQ=255)
                           tgctgcCGGCTGCCGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCC‑GTTACCTGATCGGTGACGACGAGCAc                                                                                                                  >  1:293287/1‑149 (MQ=255)
                           tgctgcCGGCTGCCGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTCTTCTTCGGCCTGTCCGGCACTGGCAACACTACCCTTTCGGCGGCCGCAAGCC‑GTTACATGCTCGGTGACGACGGGctc                                                                                                                  >  6:329773/1‑147 (MQ=255)
                                ccggctgccgATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCC‑GTTACCTGATCGGTGACGACGAGCACGGTTg                                                                                                             <  1:126942/149‑1 (MQ=255)
                                 cggctgccgATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCC‑GTTACCTGATCGGTGACGACGAGCACGGTTgg                                                                                                            >  7:27395/1‑149 (MQ=255)
                                       ccgATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCC‑GTTACCTGATCGGTGACGACGAGCACGGTTGGGGCg                                                                                                        >  4:31192/1‑147 (MQ=255)
                                       ccgATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCC‑GTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGaa                                                                                                      <  6:276703/149‑1 (MQ=255)
                                         gATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAG‑CAGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAgg                                                                                                    >  4:262725/1‑149 (MQ=255)
                                            gTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCC‑GTTAcc                                                                                                                                      >  4:107426/1‑112 (MQ=255)
                                            gTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCC‑GTTAcc                                                                                                                                      <  3:107426/112‑1 (MQ=255)
                                             tACTGCCTATGCACTGCGCGGCCAACATCAGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCC‑GTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGtt                                                                                                <  1:77186/149‑1 (MQ=255)
                                               cTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCC‑GTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGtt                                                                                                >  8:166900/1‑147 (MQ=255)
                                               cTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCC‑GTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGtt                                                                                                <  7:166900/147‑1 (MQ=255)
                                                  ccTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCC‑GTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTc                                                                                             <  6:268525/147‑1 (MQ=255)
                                                   cTATGCACTGCGCGGCCAACCTCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCC‑GTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGGCTTc                                                                                          <  4:51667/149‑1 (MQ=255)
                                                      tGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCC‑GTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCaa                                                                                        >  6:338982/1‑148 (MQ=255)
                                                         aCTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCC‑GTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATg                                                                                    >  1:363789/1‑149 (MQ=255)
                                                          cTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCC‑GTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATgg                                                                                   <  7:107561/149‑1 (MQ=255)
                                                                  ccAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCC‑GTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGtt                                                                                                <  4:152835/128‑1 (MQ=255)
                                                                  ccAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCC‑GTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGtt                                                                                                >  3:152835/1‑128 (MQ=255)
                                                                   cAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCC‑GTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGgct                                                                          >  5:132958/1‑149 (MQ=255)
                                                                   cAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTCTCCGGCACTGGCAAGACTACCCTTTAGGCTGACGAAAGCC‑GTTACCTGAACGGTGACTACGAGCACGGATCGGGCGAAGGCGTTGTCTTCAACATGTAAGGCGgct                                                                          >  8:328155/1‑149 (MQ=255)
                                                                    aaCATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCTGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCC‑GTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGgctg                                                                         >  5:87806/1‑149 (MQ=255)
                                                                        tCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCC‑GTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTAc                                                                     >  6:365806/1‑149 (MQ=255)
                                                                           gCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAG‑CAGTTACCTGATCGGTGACGACGAGCACGGATGGGGCGAAGCCGTTGTCATCAACATGGAATGCGGCTGCTActc                                                                   >  7:363618/1‑146 (MQ=255)
                                                                             gaagaaGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCC‑GTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCaa                                                                <  8:27395/149‑1 (MQ=255)
                                                                             gaagaaGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCC‑GTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCaa                                                                <  8:284056/149‑1 (MQ=255)
                                                                                   ggCGATGTGACCTTGTTCTTCGGCCTGTCCTGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCC‑GTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCAt                                                          <  6:87806/149‑1 (MQ=255)
                                                                                    gCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCC‑GTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATc                                                         <  5:262999/149‑1 (MQ=255)
                                                                                       aTGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCC‑GTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGAc                                                      <  2:363789/149‑1 (MQ=255)
                                                                                             ccTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCC‑GTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTcc                                                >  2:126949/1‑149 (MQ=255)
                                                                                              cTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCC‑GTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTcc                                                >  2:107356/1‑148 (MQ=255)
                                                                                               ttGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCC‑GTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGt                                                              >  8:92726/1‑133 (MQ=255)
                                                                                               ttGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCC‑GTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGt                                                              <  7:92726/133‑1 (MQ=255)
                                                                                                 tttcttcGGCCTGTCCGTCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCC‑GTTACCTGTTCGGTGACGACGAGGACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCgag                                             <  5:329773/147‑1 (MQ=255)
                                                                                                  ttcttcGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCC‑GTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCgaga                                            <  3:31192/148‑1 (MQ=255)
                                                                                                              tCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCC‑GTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCgg                                <  3:262725/148‑1 (MQ=255)
                                                                                                              tCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCC‑GTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGt                               >  5:356880/1‑149 (MQ=255)
                                                                                                                    aCTGGCAAGACTACCCTTTCGGCTGACGAAAGCC‑GTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATCTg                         <  3:182131/149‑1 (MQ=255)
                                                                                                                        gCAAGACTACCCTTTCGGCTGACGAAAGCC‑GTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATCTGGaaa                     <  5:11867/149‑1 (MQ=255)
                                                                                                                          aaGACTACCCTTTCGGCTGACGAAAGCC‑GTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATCTGGAAAGc                   >  3:253189/1‑149 (MQ=255)
                                                                                                                           aGACTACCCTTTCGGCTGACGAAAGCC‑GTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATCTGGAAAGcc                  <  6:356880/149‑1 (MQ=255)
                                                                                                                                 cccTTTCGGCTGACGAAAGCC‑GTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATCTGGAAAGCCATCa              <  5:365806/147‑1 (MQ=255)
                                                                                                                                 cccTTTCGGCTGACGAAAGCC‑GTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATCTGGAAAGCCATCAAg            >  6:274263/1‑149 (MQ=255)
                                                                                                                                 cccTTTCGGCTGACGAAAGCC‑GTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATCTGGAAAGCCATCAAg            >  4:336926/1‑149 (MQ=255)
                                                                                                                                  ccTTTCGGCTGACGAAAGCC‑GTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCAt                                                          <  6:307018/102‑1 (MQ=255)
                                                                                                                                  ccTTTCGGCTGACGAAAGCC‑GTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCAt                                                          >  5:307018/1‑102 (MQ=255)
                                                                                                                                        ggCTGACGAAAGCC‑GTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATCTGGAAAGCCATCAAGCACGgcg     <  6:46704/149‑1 (MQ=255)
                                                                                                                                            gACGAAAGCC‑GTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATCTGGAAAGCCATCAAGAACGGCGAAg  >  4:230926/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                            
AAGCCATGTTCTCGGTGCAGAACTTCCTGCTGCCGGCTGCCGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAG‑C‑GTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATCTGGAAAGCCATCAAGCACGGCGCAG  >  NC_002947/307718‑308004

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