Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 390,996 T→A 6.4% Y424N (TAC→AAC)  soxA → sarcosine oxidase subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947390,9960TA 6.4% 190.5 / 3.1 62Y424N (TAC→AAC) soxAsarcosine oxidase subunit alpha
Reads supporting (aligned to +/- strand):  ref base T (30/28);  new base A (2/2);  total (32/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.64e-01

GCATACTCACCTCCAGCGCCGTGATCGAGGCCAAAGGCAGCAAGCACGTCACCGGCGCCCGTGTGGCGGCGATCGATGTGCAGGCGCACAAAGTCACCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGCCGGGCGATGCGCCGCAGAAGCGTGAGTGCGTCGGTGGTATCAACGGCGTGTATGCCTTGGGCGATGTCATTG  >  NC_002947/390850‑391143
                                                                                                                                                  |                                                                                                                                                   
gCATACTCACCTGCAGCTCCTTGATCGAGGCCAAAGGCAGCAAGCACGTCACCGGCTCCCGTGTGGCGGCGATCGATGTGCAGGCGCACAAAGTCACGAGGCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTa                                                                                                                                                    <  3:243501/148‑1 (MQ=255)
   tACTCACCTCCAGCGCCGTGATCGAGGCCAAAGGCAGCAAGCACGTCACCGGCGCCCGTGTGGCGGCGATCGATGTGCAGGCGCACAAAGTCACCAGCCCAAGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGc                                                                                                                                                <  1:127842/149‑1 (MQ=255)
      tCACCTCCAGCGCCGTGATCGAGGCCAAAGGCAGCAAGCACGTCACCGGCGCCCGTGTGGCGGCGATCGATGTGCAGGCGCACAAAGTCACCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCTGTAACCGccc                                                                                                                                              >  6:288783/1‑148 (MQ=255)
      tCACCTCCAGCGCCGTGATCGAGGCCAAAGGCAGCAAGCACGTCACCGGCGCCCGTGTGGCGGCGATCGATGTGCAGGCGCACAAAGTCACCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGccc                                                                                                                                              >  5:116252/1‑148 (MQ=255)
        aCCTCCAGCGCCGTGATCGAGGCCAAAGGCAGCAAGCACGTCACCGGCGCCCGTGTGGCGGCGATCGATGTGCAGGCGCACAAAGTCACCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCTGTAACCGCCCGag                                                                                                                                           <  5:288783/149‑2 (MQ=255)
         ccTCCAGCGCCGTGATCGAGGCCAAAGGCAGCAAGCACGTCACCGGCGCCCGTGTGGCGGCGATCGATGTGCAGGCGCACAAAGTCACCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATc                                                                                                                                          <  6:116252/149‑1 (MQ=255)
             cAGCGCCGTGATCGAGGCCAAAGGCAGCAAGCACGTCACCGGCGCCCGTGTGGCGGCGATCGATGTGCAGGCGCACAAAGTCACCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTg                                                                                                                                       <  3:222825/148‑1 (MQ=255)
               gcgcCGTGATCGAGGCCAAAGGCAGCAAGCACGTCACCGGCGCCCGTGTGGCGGCGATCGATGTGCAGGCGCACAAAGTCACCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCAc                                                                                                                                    >  7:175511/1‑149 (MQ=255)
                       aTCGAGGCCAAAGGCAGCAAGCACGTCACCGGCGCCCGTGTGGCGGCGATCGATGTGCAGGCGCACAAAGTCACCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTc                                                                                                                            <  2:235677/149‑1 (MQ=255)
                        tCGAGGCCAAAGGCAGCAAGCACGTCACCGGCGCCCGTGTGGCGGCGATCGATGTGCAGGCGCACAAAGTCACCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCg                                                                                                                           <  1:286673/149‑1 (MQ=255)
                          gAGGCCAAAGGCAGCAAGCACGTCACCGGCGCCCGTGTGGCGGCGATCGATGTGCAGGCGCACAAAGTCACCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCa                                                                                                                         >  7:261267/1‑149 (MQ=255)
                                aaaGGCAGCAAGCACGTCACCGGCGCCCGTGTGGCGGCGATCGATGTGCAGGCGCACAAAGTCACCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCTGTAACCGCCCGag                                                                                                                                           >  4:135952/1‑124 (MQ=255)
                                aaaGGCAGCAAGCACGTCACCGGCGCCCGTGTGGCGGCGATCGATGTGCAGGCGCACAAAGTCACCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCTGTAACCGCCCGag                                                                                                                                           <  3:135952/125‑2 (MQ=255)
                                 aaGGCAGCAAGCACGTCACCGGCGCCCGTGTGGCGGCGATCGATGTGCAGGCGCACAAAGTCACCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGc                                                                                                                  <  3:219964/149‑1 (MQ=255)
                                  aGGCAGCAAGCACGTCACCGGCGCCCGTGTGGCGGCGATCGATGTGCAGGCGCACAAAGTCACCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCg                                                                                                                 >  6:88442/1‑149 (MQ=255)
                                        caagcaCGTCACCGGCGCCCGTGTGGCGGCGATCGATGTGCAGGCGCACAAAGTCACCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCACCt                                                                                                                                  <  7:226394/126‑1 (MQ=255)
                                        caagcaCGTCACCGGCGCCCGTGTGGCGGCGATCGATGTGCAGGCGCACAAAGTCACCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCACCt                                                                                                                                  >  8:226394/1‑126 (MQ=255)
                                         aagcaCGTCACCGGCGCCCGTGTGGCGGCGATCGATGTGCAGGCGCACAAAGTCACCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGccc                                                                                                          >  8:127391/1‑149 (MQ=255)
                                            caCGTCACCGGCGCCCGTGTGGCGGCGATCGATGTGCAGGCGCACAAAGTCACCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGt                                                                                                       <  5:88442/149‑1 (MQ=255)
                                             aCGTCACCGGCGCCCGTGTGGCGGCGATCGATGTGCAGGCGCACAAAGTCACCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTa                                                                                                      <  8:175511/149‑1 (MQ=255)
                                              cGTCACCGGCGCCCGTGTGGCGGCGATCGATGTGCAGGCGCACAAAGTCACCAGCCCAGGCGAAGTCCTTGAGTCCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTAt                                                                                                     >  8:274113/1‑149 (MQ=255)
                                                  aCCGGCGCCCGTGTGGCGGCGATCGATGTGCAGGCGCACAAAGTCACCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCAc                                                                                                                                    <  1:86814/114‑1 (MQ=255)
                                                  aCCGGCGCCCGTGTGGCGGCGATCGATGTGCAGGCGCACAAAGTCACCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCAc                                                                                                                                    >  2:86814/1‑114 (MQ=255)
                                                   ccGGCGCCCGTGTGGCGGCGATCGATGTGCAGGCGCACAAAGTCACCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATg                                                                                                    <  6:162154/145‑1 (MQ=255)
                                                   ccGGCGCCCGTGTGGCGGCGATCGATGTGCAGGCGCACAAAGTCACCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATg                                                                                                    >  5:162154/1‑145 (MQ=255)
                                                   ccGGCGCCCGTGTGGCGGCGATCGATGTGCAGGCGCACAAAGTCACCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGt                                                                                                <  6:363016/149‑1 (MQ=255)
                                                                      gatcgatGTGCAGGCGCACAAAGTCACCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCg                                                                             >  2:126109/1‑149 (MQ=255)
                                                                           atGTGCAGGCGCACAAAGTCACCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCAGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGCCg                                                                        >  8:141980/1‑149 (MQ=255)
                                                                            tgtgCAGGCGCACAAAGTCACCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGCCgg                                                                       <  8:261267/149‑1 (MQ=255)
                                                                                 aGGCGCACAAAGTCACCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGCCGGGCGAt                                                                  >  2:127823/1‑149 (MQ=255)
                                                                                           aGTCACCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGCCGGGCGATGCGCCGCAGa                                                        <  1:126109/149‑1 (MQ=255)
                                                                                                 cAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGcc                                                                                                                                               >  7:205540/1‑56 (MQ=255)
                                                                                                 cAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGcc                                                                                                                                               <  8:205540/56‑1 (MQ=255)
                                                                                                 cAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGcc                                                                         <  6:317810/126‑1 (MQ=255)
                                                                                                 cAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGCCGGGCGATGCGCCGCAGAAGCGt                                                   >  4:330989/1‑148 (MQ=255)
                                                                                                 cAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTGCGTGcc                                                                         >  5:317810/1‑126 (MQ=255)
                                                                                                  aGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCAGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGCCGGGCGATGCGCCGCAGAAGCGTGa                                                 <  7:141980/149‑1 (MQ=255)
                                                                                                    cccAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGCCGGGCGATGCGCCGCAGAAGCGTGAGt                                               <  7:138566/149‑1 (MQ=255)
                                                                                                       aGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGCCGGGCGATGCGCCGCAGAAGCGTGAGTGCg                                            >  5:255398/1‑149 (MQ=255)
                                                                                                         gCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATcc                                                                                      <  2:241628/105‑1 (MQ=255)
                                                                                                          cGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGCCGGGCGATGCGCCGCAGAAGCGTGAGTGCGTCg                                         >  5:90945/1‑149 (MQ=255)
                                                                                                             aGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGCCGGGCGATGCGCCGCAGAAGCGTGAGTGCGTCggtg                                      <  6:255398/149‑1 (MQ=255)
                                                                                                             aGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGCCGGGCGATGCGCCGCAGAAGCGTGAGTGCGTCggtg                                      >  2:207700/1‑149 (MQ=255)
                                                                                                                  ttGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGCCGGGCGATGCGCCGCAGAAGCGTGAGTGCGTCGGTGGTAt                                  >  6:73445/1‑148 (MQ=255)
                                                                                                                      gTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGCCGGGCGATGCGCCGCAGAAGCGTGAGTGCGTCGGTGGTATCAACg                             <  1:127823/149‑1 (MQ=255)
                                                                                                                               gATCGCCTCCTCGGGCGGTTAGAGCCGGTTCGTGCTCCTTGCTTCGCACCTGGGGGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGCCGGGCGTTGAGCGTCAGAAGCGTGAGTGCGTCGTTGGTATCAACGGCgt                         <  1:51806/144‑1 (MQ=255)
                                                                                                                               gATCGCCTCCTCGGGCGGTTACAGCGCGATCGTTCACCTCGCTTCGCACCTGGGCGGTCGCCCGGTATGGCCTGACGACATCCTTCGCTTCGTGCCGCGCGATGCGCCGCAGAAGCGTGACTCCGAGGGTGGTATCAACGGCgt                         >  2:51806/1‑144 (MQ=255)
                                                                                                                               gATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGCCGGGCGATGCGCCGCAGAAGCGTGAGTGCGTCGGTGGTATCAACGGc                           <  2:313852/142‑1 (MQ=255)
                                                                                                                               gATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGCCGGGCGATGCGCCGCAGAAGCGTGAGTGCGTCGGTGGTATCAACGGc                           >  1:313852/1‑142 (MQ=255)
                                                                                                                                aTCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGCCGGGCGATGCGCCGCAGAAGCGTGAGTGCGTCGGTGGTATCAACGGCGTGTATGc                   >  2:306870/1‑149 (MQ=255)
                                                                                                                                 tCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGCCGGGCGATGCGCCGCAGAAGCGTGAGTGCGTCgg                                        <  8:197103/127‑1 (MQ=255)
                                                                                                                                    cctcctCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGCCGGGCGATGCGCCGCAGAAGCGTGAGTGCGTCGGTGGTATCAACGGCGTGTATGCCTTg               >  3:258896/1‑149 (MQ=255)
                                                                                                                                    cctcctCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGCCGGGCGATGCGCCGCAGAAGCGTGAGTGCGTCGGTGGTATCAACGGCGTGTATGCCTTg               >  7:297503/1‑149 (MQ=255)
                                                                                                                                     ctcctcGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGCCGGGCGATGCGCCGCAGAAGCGTGAGTGCGTCGGTGGTATCAACGGCGTGTATGCCTTgg              >  1:334677/1‑149 (MQ=255)
                                                                                                                                      tcctcGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGCCGGGCGATGCGCCGCAGAAGCGTGAGTGCGTCGGTGGTATCAACGGCGTGTATGCCTTggg             >  8:149444/1‑149 (MQ=255)
                                                                                                                                      tcctcGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGCCGGGCGATGCGCCGCAGAAGCGTGAGTGCGTCGGTGGTATCAACGGCGTGTATGCCTTggg             >  2:237778/1‑149 (MQ=255)
                                                                                                                                        ctcGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGCCGGGCGATGCGCCGCAGAAGCGTGAGTGCGTCGGTGGTATCAACGGCGTGTATGCCTTggg             <  1:207700/147‑1 (MQ=255)
                                                                                                                                        ctcGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGCCGGGCGATGCGCCGCAGAAGCGTGAGTGCGTCGGTGGTATCAACGGCGTGTATGCCTTGGGc            >  3:245991/1‑148 (MQ=255)
                                                                                                                                           gggCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGCCGGGCGATGCGCCGCAGAAGCGTGAGTGCGTCGGTGGTATCAACGGCGTGTATGCCTTGGGCGATg        >  3:67328/1‑149 (MQ=255)
                                                                                                                                                gTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGCCGGGCGATGCGCCGCAGAAGCGTGAGTGCGTCGGTGGTATCAACGGCGTGTATGCCTTGGGCGATGTCAtt   <  7:127391/149‑1 (MQ=255)
                                                                                                                                                gTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGCCGGGCGATGCGCCGCAGAAGCGTGAGTGCGTCGGTGGTATCAACGGCGTGTATGCCTTGGGCGATGTCAtt   <  4:258896/149‑1 (MQ=255)
                                                                                                                                                gTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGCCGGGCGATGCGCCGCAGAAGCGTGAGTGCGTCGGTGGTATCAACGGCGTGTATGCCTTGGGCGATGTCAtt   >  2:274553/1‑149 (MQ=255)
                                                                                                                                                 ttACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGCCGGGCGATGCGCCGCAGAAGCGTGAGTGCGTCGGTGGTATCAACGGCGTGTATGCCTTGGGCGATGTCATTg  >  8:121191/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GCATACTCACCTCCAGCGCCGTGATCGAGGCCAAAGGCAGCAAGCACGTCACCGGCGCCCGTGTGGCGGCGATCGATGTGCAGGCGCACAAAGTCACCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGCCGGGCGATGCGCCGCAGAAGCGTGAGTGCGTCGGTGGTATCAACGGCGTGTATGCCTTGGGCGATGTCATTG  >  NC_002947/390850‑391143

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