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.5% 229.5 / ‑3.2 66intergenic (+880/‑989)hslO/PP_0254molecular chaperone HslO/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (32/33);  total (32/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.96e-01

AGCCATGTTCTCGGTGCAGAACTTCCTGCTGCCGGCTGCCGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTAC‑‑CCTTTCGGCTGACGAAAGCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATCTGGAAAGCCATCAAGCACGGCGCAGTACT  >  NC_002947/307719‑308008
                                                                                                                                                      |                                                                                                                                               
aGCCATGTTCTCGGTGCAGAACTTCCTGCTGCCGGCTGCCGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTAC‑‑CCTTTCGGCTGACGAAAGcc                                                                                                                                               <  3:187835/149‑2 (MQ=255)
   cATGTTCTCGGTGCAGAACTTCCTGCTGCCGGCTGCCGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTAC‑‑CCTTTCGGCTGACGAAAGccgtt                                                                                                                                            <  4:291528/149‑5 (MQ=255)
            ggTGCAGAACTTCCTGCTGCCGGCTGCCGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTAC‑‑CCTTTCGGCTGACGAAAGCCGTTACCTGATc                                                                                                                                     <  5:340485/148‑1 (MQ=255)
                    aCTTCCTGCTGCCGGCTGCCGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTAC‑‑CCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTgacgac                                                                                                                            <  7:308477/149‑1 (MQ=255)
                                cggctgccgATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTAC‑‑CCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTgg                                                                                                                <  2:303294/149‑1 (MQ=255)
                                  gctgccgATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTAC‑‑CCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTgggg                                                                                                              <  2:204948/149‑1 (MQ=255)
                                  gctgccgATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTAC‑‑CCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTgggg                                                                                                              >  6:191589/1‑149 (MQ=255)
                                      ccgATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTAC‑‑CCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGaa                                                                                                          <  5:191589/149‑1 (MQ=255)
                                       cgATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTAC‑‑CCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAg                                                                                                         >  3:114656/1‑149 (MQ=255)
                                       cgATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTAC‑‑CCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTGGGGGGGagg                                                                                                         >  2:327663/1‑147 (MQ=255)
                                        gATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTAC‑‑CCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAgg                                                                                                        >  5:39888/1‑149 (MQ=255)
                                        gATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTAC‑‑CCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAgg                                                                                                        >  5:45303/1‑149 (MQ=255)
                                         aTGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTAC‑‑CCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGc                                                                                                       <  3:217742/149‑1 (MQ=255)
                                            tACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTAC‑‑CCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGtt                                                                                                    <  6:243040/149‑1 (MQ=255)
                                             aCTGCCTATGCACTGCGCGGCCAACATGGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTAC‑‑CCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGtt                                                                                                    <  8:72129/148‑1 (MQ=255)
                                              cTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTAC‑‑CCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTga                                                                                                                                <  8:106786/119‑1 (MQ=255)
                                              cTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTAC‑‑CCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTga                                                                                                                                >  7:106786/1‑119 (MQ=255)
                                              cTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTAC‑‑CCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGt                                                                                                  >  7:147293/1‑149 (MQ=255)
                                                 ccTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTAC‑‑CCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCtt                                                                                               >  5:283537/1‑149 (MQ=255)
                                                  cTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTAC‑‑CCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCtt                                                                                               >  4:89838/1‑148 (MQ=255)
                                                  cTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTAC‑‑CCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCtt                                                                                               >  4:84722/1‑148 (MQ=255)
                                                     tGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTAC‑‑CCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCaa                                                                                            <  8:169918/148‑1 (MQ=255)
                                                     tGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTAC‑‑CCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCa                                                                                             <  2:139696/147‑1 (MQ=255)
                                                        aCTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTAC‑‑CCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATg                                                                                        <  2:129748/149‑1 (MQ=255)
                                                            cgcgGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTAC‑‑CCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGaa                                                                                     <  7:151906/148‑1 (MQ=255)
                                                             gcgGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTAC‑‑CCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAgg                                                                                   >  3:35509/1‑149 (MQ=255)
                                                                       tCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTAC‑‑CCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTAc                                                                         <  4:114656/149‑1 (MQ=255)
                                                                              agaagGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTAC‑‑CCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAg                                                                   <  6:283537/148‑1 (MQ=255)
                                                                               gaagGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCTCCCTTCGGCTGTTGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAg                                                                   <  1:120616/149‑1 (MQ=255)
                                                                                     gATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTAC‑‑CCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGa                                                           >  6:45662/1‑149 (MQ=255)
                                                                                      aTGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTAC‑‑CCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGAc                                                          <  4:35509/149‑1 (MQ=255)
                                                                                      aTGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTAC‑‑CCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGAc                                                          <  3:89838/149‑1 (MQ=255)
                                                                                      aTGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTAC‑‑CCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGAc                                                          <  3:84722/149‑1 (MQ=255)
                                                                                           aCCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTAC‑‑CCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTc                                                     >  7:200490/1‑149 (MQ=255)
                                                                                            ccTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTAC‑‑CCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTcc                                                    <  6:43731/149‑1 (MQ=255)
                                                                                             cTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTAC‑‑CCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCAt                                                              <  7:263746/138‑1 (MQ=255)
                                                                                             cTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTAC‑‑CCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCAt                                                              >  8:263746/1‑138 (MQ=255)
                                                                                               tGTTCTTCGGCCTGTCCGGCACTGGCAAGACTAC‑‑CCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCaa                                                                    <  4:229134/130‑1 (MQ=255)
                                                                                               tGTTCTTCGGCCTGTCCGGCACTGGCAAGACTAC‑‑CCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCaa                                                                    >  3:229134/1‑130 (MQ=255)
                                                                                               tGTTCTTCGGCCTGTCCGGCACTGGCAAGACTAC‑‑CCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTcc                                                    >  7:10690/1‑146 (MQ=255)
                                                                                                gTTCTTCGGCCTGTCCGGCACTGGCAAGACTAC‑‑CCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCgag                                                 >  4:178394/1‑148 (MQ=255)
                                                                                                gTTCTTCGGCCTGTCCGGCACTGGCAAGACTAC‑‑CCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCgag                                                 >  7:327915/1‑148 (MQ=255)
                                                                                                   cttcGGCCTGTCCGGCACTGGCAAGACTAC‑‑CCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCgaga                                                <  2:4450/146‑1 (MQ=255)
                                                                                                   cttcGGCCTGTCCGGCACTGGCAAGACTAC‑‑CCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCgaga                                                >  1:4450/1‑146 (MQ=255)
                                                                                                      cGGCCTGTCCGGCACTGGCAAGACTAC‑‑CCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAAc                                           <  8:147293/148‑1 (MQ=255)
                                                                                                            gTCCGGCACTGGCAAGACTAC‑‑CCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCgg                                    >  8:209794/1‑149 (MQ=255)
                                                                                                                ggCACTGGCAAGACTAC‑‑CCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCa                                 >  3:274879/1‑148 (MQ=255)
                                                                                                                 gCACTGGCAAGACTAC‑‑CCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCg                                                                                                      <  7:5734/78‑1 (MQ=255)
                                                                                                                 gCACTGGCAAGACTAC‑‑CCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCg                                                                                                      >  8:5734/1‑78 (MQ=255)
                                                                                                                 gCACTGGCAAGACTAC‑‑CCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATc                               <  8:188436/149‑1 (MQ=255)
                                                                                                                 gCACTGGCAAGACTAC‑‑CCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATc                               >  2:246845/1‑149 (MQ=255)
                                                                                                                   aCTGGCAAGACTAC‑‑CCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATCTg                             >  1:184002/1‑149 (MQ=255)
                                                                                                                     tGGCAAGACTAC‑‑CCTTTCGGCTGACGAAAGCCGTTGCCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATCTgg                            >  8:265235/1‑148 (MQ=255)
                                                                                                                       gCAAGACTAC‑‑CCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAg                                                                   >  7:179944/1‑107 (MQ=255)
                                                                                                                       gCAAGACTAC‑‑CCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAg                                                                   <  8:179944/107‑1 (MQ=255)
                                                                                                                        cAAGACTAC‑‑CCTTTCTGCTGACGAAATCCGTTACCTGATCTGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATCTGGAAAg                        <  8:327915/149‑1 (MQ=255)
                                                                                                                        cAAGACTAC‑‑CCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGATCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGTCTGCTACGCCAAGTGCATCGACCTGTCAGAGTTGAACGAGCCGGTCATCTGGaaa                         <  7:57288/148‑1 (MQ=255)
                                                                                                                        cAAGACTAC‑‑CCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATCTGGAAAg                        <  1:327663/149‑1 (MQ=255)
                                                                                                                             cTAC‑‑CCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATCTGGAAAGCCATc                   <  6:45303/149‑1 (MQ=255)
                                                                                                                             cTAC‑‑CCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATCTGGAAAGCCATc                   <  6:39888/149‑1 (MQ=255)
                                                                                                                                     tttCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATCTGGAAAGCCATCAAGCAc             >  7:31516/1‑149 (MQ=255)
                                                                                                                                         ggCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATCTGGAAAGCCATCAAGCACGgcg         >  5:346275/1‑149 (MQ=255)
                                                                                                                                         ggCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATCTGGAAAGCCATCAAGCACGgcg         <  3:178394/149‑1 (MQ=255)
                                                                                                                                          gCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATCTGGAAAGCCATCAAGCACGgcg         >  5:156415/1‑148 (MQ=255)
                                                                                                                                                gAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTcc                                                    >  6:211753/1‑99 (MQ=255)
                                                                                                                                                gAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTcc                                                    <  5:211753/99‑1 (MQ=255)
                                                                                                                                                gAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATCTGGAAAGCCATCAAGCACGGCGCAGTACt  >  4:51350/1‑149 (MQ=255)
                                                                                                                                                gAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATCTGGAAAGCCATCAAGCACGGCGCAGTACt  >  1:348633/1‑149 (MQ=255)
                                                                                                                                                      |                                                                                                                                               
AGCCATGTTCTCGGTGCAGAACTTCCTGCTGCCGGCTGCCGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTAC‑‑CCTTTCGGCTGACGAAAGCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATCTGGAAAGCCATCAAGCACGGCGCAGTACT  >  NC_002947/307719‑308008

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: