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.C98.4% 207.5 / ‑3.3 62intergenic (+880/‑989)hslO/PP_0254molecular chaperone HslO/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (28/33);  total (28/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.27e-01

AAGCCATGTTCTCGGTGCAGAACTTCCTGCTGCCGGCTGCCGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATCTGGAAAGCCATCAAGCACGGCGCAGTAC  >  NC_002947/307718‑308007
                                                                                                                                                    |                                                                                                                                              
aaGCCATGTTCTCGGTGCAGAACTTCCTGCTGCCGGCTGCCGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGc                                                                                                                                               >  5:349660/1‑149 (MQ=255)
 aGCCATGTTCTCGGTGCAGAACTTCCTGCTGCCGGCTGCCGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGcc                                                                                                                                              <  5:308434/149‑2 (MQ=255)
   ccATGTTCTCGGTGCAGAACTTCCTGCTGCCGGCTGCCGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGTCAGCCGAAAGGGta                                                                                                                                            <  7:272521/149‑2 (MQ=255)
       gTTCTCGGTGGAGATCTTCCTGCTGGCGGCTTCCGATTTACTGCCTATGCACTGCTCGGCCAACATCGGCGAAGATGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGccgtta                                                                                                                                          <  6:165793/147‑6 (MQ=255)
             ggTGCAGAACTTCCTGCTGCCGGCTGCCGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCg                                                                                                                                   >  7:30582/1‑149 (MQ=255)
              gTGCAGAACTTCCTGCTGCCGGCTGCCGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCgg                                                                                                                                  <  8:316108/149‑1 (MQ=255)
                  aGAACTTCCTGCTGCCGGCTGCCGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTTGGCTGACGAAAGCCGTTACCTGATCGGTgac                                                                                                                              <  6:307423/149‑1 (MQ=255)
                   gAACTTCCTGCTGCCGGCTGCCGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTgac                                                                                                                              >  1:161056/1‑148 (MQ=255)
                         ccTGCTGCCGGCTGCCGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGAcgacga                                                                                                                         <  3:347920/147‑1 (MQ=255)
                         ccTGCTGCCGGCTGCCGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGAcgacga                                                                                                                         <  6:359240/147‑1 (MQ=255)
                            gctgccggctgccgATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACg                                                                                                                    >  4:242655/1‑149 (MQ=255)
                             ctgccggctgccgATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATc                                                                                                                                    <  1:259725/132‑1 (MQ=255)
                             ctgccggctgccgATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATc                                                                                                                                    >  2:259725/1‑132 (MQ=255)
                             ctgccggctgccgATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCAc                                                                                                                     <  3:248856/147‑1 (MQ=255)
                                 cggctgccgATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTgg                                                                                                               >  3:181023/1‑149 (MQ=255)
                                  ggctgccgATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTggg                                                                                                              <  5:115000/149‑1 (MQ=255)
                                       ccgATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCg                                                                                                           >  7:62531/1‑147 (MQ=255)
                                        cgATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAg                                                                                                        >  2:374636/1‑149 (MQ=255)
                                         gATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAgg                                                                                                       >  4:36399/1‑149 (MQ=255)
                                            gTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGt                                                                                                    >  7:208805/1‑149 (MQ=255)
                                             tACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGtt                                                                                                   <  1:374636/149‑1 (MQ=255)
                                                 gCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCt                                                                                               >  7:83511/1‑149 (MQ=255)
                                                  ccTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCtt                                                                                              <  8:83511/149‑1 (MQ=255)
                                                    tATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGt                                                                                                    <  1:270437/141‑1 (MQ=255)
                                                    tATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGAGt                                                                                                    >  2:270437/1‑141 (MQ=255)
                                                          cTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATgg                                                                                      <  2:82894/149‑1 (MQ=255)
                                                           tGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGtt                                                                                                   <  1:384587/135‑1 (MQ=255)
                                                           tGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGtt                                                                                                   >  2:384587/1‑135 (MQ=255)
                                                                ggCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAggcg                                                                                >  1:387663/1‑149 (MQ=255)
                                                                  ccAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAggcggc                                                                              <  4:181023/149‑1 (MQ=255)
                                                                   cAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGgct                                                                             >  4:39993/1‑149 (MQ=255)
                                                                    aaCATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGgctg                                                                            <  8:62531/149‑1 (MQ=255)
                                                                      cATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCAATGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGgctgc                                                                           >  2:70030/1‑148 (MQ=255)
                                                                       aTCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTa                                                                         >  2:354215/1‑149 (MQ=255)
                                                                           gCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGc                                                                      >  2:340346/1‑148 (MQ=255)
                                                                           gCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGc                                                                      >  5:26916/1‑148 (MQ=255)
                                                                           gCGAAGAAGGCGATGTGACCTTGTTCTTCAGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGgctg                                                                            <  3:91580/142‑1 (MQ=255)
                                                                           gCGAAGAAGGCGATGTGACCTTGTTCTTCAGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGgctg                                                                            >  4:91580/1‑142 (MQ=255)
                                                                               agaagGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGt                                                                 <  3:36399/149‑1 (MQ=255)
                                                                               agaagGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGt                                                                 <  3:39993/149‑1 (MQ=255)
                                                                                 aagGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGc                                                               <  2:387663/149‑1 (MQ=255)
                                                                                 aagGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGc                                                               <  3:242655/149‑1 (MQ=255)
                                                                                   ggCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATgg                                                                                      >  5:12710/1‑124 (MQ=255)
                                                                                   ggCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATgg                                                                                      <  6:12710/124‑1 (MQ=255)
                                                                                       aTGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGAc                                                         <  1:354215/149‑1 (MQ=255)
                                                                                       aTGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGAc                                                         <  3:366964/149‑1 (MQ=255)
                                                                                            aCCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTc                                                    <  2:161056/149‑1 (MQ=255)
                                                                                              cTTGTTCTTCGGCCTGTCCGGCAATGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTcc                                                   <  1:70030/148‑1 (MQ=255)
                                                                                                     ttcGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAgaagaa                                           >  1:373820/1‑149 (MQ=255)
                                                                                                        ggCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGa                                        <  8:30582/149‑1 (MQ=255)
                                                                                                          ccTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGc                                      <  2:64439/149‑1 (MQ=255)
                                                                                                              tCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCg                                    <  8:100400/147‑1 (MQ=255)
                                                                                                                 ggCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCAt                               >  6:54254/1‑149 (MQ=255)
                                                                                                                   cACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATGGGTGAGGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGTAAGGCGGCTGCTACGCCAAGTGCATCTACCTTTACGAGAAGAACGAGACGGTCATCt                             <  6:349660/149‑1 (MQ=255)
                                                                                                                      tGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATCTgg                           >  6:117648/1‑148 (MQ=255)
                                                                                                                      tGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATCTgg                           >  6:48625/1‑148 (MQ=255)
                                                                                                                        gCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATCTGGa                          >  2:133550/1‑147 (MQ=255)
                                                                                                                              cTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATCTGGAAAGCCAt                   <  5:54254/148‑1 (MQ=255)
                                                                                                                                  ccTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATCTGGAAAGCCATCAAGc              >  7:325188/1‑149 (MQ=255)
                                                                                                                                   cTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATCTGGAAAGCCATCAAGCa             <  8:208805/149‑1 (MQ=255)
                                                                                                                                         gCTGACGAAAGCCGTTACCTGATCGGTGATGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATc                              <  3:378533/126‑1 (MQ=255)
                                                                                                                                         gCTGACGAAAGCCGTTACCTCATCGGTGATGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACCCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATc                              >  4:378533/1‑126 (MQ=255)
                                                                                                                                           tGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATCTGGAAAGCCATCAAGCACGGCGCAg     >  7:220490/1‑149 (MQ=255)
                                                                                                                                              cGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATCTGGAAAGCCATCAATCACGGCGCAGTa   <  6:26916/148‑1 (MQ=255)
                                                                                                                                              cGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATCTGGAAAGCCATCAAGCACGGCGCAGTAc  <  8:220490/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
AAGCCATGTTCTCGGTGCAGAACTTCCTGCTGCCGGCTGCCGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATCTGGAAAGCCATCAAGCACGGCGCAGTAC  >  NC_002947/307718‑308007

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: