Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 390,993 G→T 6.9% G423C (GGT→TGT)  soxA → sarcosine oxidase subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947390,9930GT 6.9% 174.7 / 3.3 59G423C (GGT→TGT) soxAsarcosine oxidase subunit alpha
Reads supporting (aligned to +/- strand):  ref base G (28/25);  new base T (2/2);  total (31/28)
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.36e-01

GCATACTCACCTCCAGCGCCGTGATCGAGGCCAAAGGCAGCAAGCACGTCACCGGCGCCCGTGTGGCGGCGATCGATGTGCAGGCGCACAAAGTCACCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGCCGGGCGATGCGCCGCAGAAGCGTGAGTGCGTCGGTGGTATCAACGGCGTGTATGCCTTGGGCGATG  >  NC_002947/390850‑391137
                                                                                                                                               |                                                                                                                                                
gCATACTCACCTGCAGCTCCTTGATCGAGGCCAAAGGCAGCAAGCACGTCACCGGCTCCCGTGTGGCGGCGATCGATGTGCAGGCGCACAAAGTCACGAGGCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTa                                                                                                                                              <  3:243501/148‑1 (MQ=255)
   tACTCACCTCCAGCGCCGTGATCGAGGCCAAAGGCAGCAAGCACGTCACCGGCGCCCGTGTGGCGGCGATCGATGTGCAGGCGCACAAAGTCACCAGCCCAAGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGc                                                                                                                                          <  1:127842/149‑1 (MQ=255)
      tCACCTCCAGCGCCGTGATCGAGGCCAAAGGCAGCAAGCACGTCACCGGCGCCCGTGTGGCGGCGATCGATGTGCAGGCGCACAAAGTCACCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCgg                                                                                                                                                 <  8:305841/139‑1 (MQ=255)
      tCACCTCCAGCGCCGTGATCGAGGCCAAAGGCAGCAAGCACGTCACCGGCGCCCGTGTGGCGGCGATCGATGTGCAGGCGCACAAAGTCACCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCgg                                                                                                                                                 >  7:305841/1‑139 (MQ=255)
      tCACCTCCAGCGCCGTGATCGAGGCCAAAGGCAGCAAGCACGTCACCGGCGCCCGTGTGGCGGCGATCGATGTGCAGGCGCACAAAGTCACCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCTGTAACCGccc                                                                                                                                        >  6:288783/1‑148 (MQ=255)
      tCACCTCCAGCGCCGTGATCGAGGCCAAAGGCAGCAAGCACGTCACCGGCGCCCGTGTGGCGGCGATCGATGTGCAGGCGCACAAAGTCACCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGccc                                                                                                                                        >  5:116252/1‑148 (MQ=255)
        aCCTCCAGCGCCGTGATCGAGGCCAAAGGCAGCAAGCACGTCACCGGCGCCCGTGTGGCGGCGATCGATGTGCAGGCGCACAAAGTCACCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCgg                                                                                                                                                 >  1:175004/1‑137 (MQ=255)
        aCCTCCAGCGCCGTGATCGAGGCCAAAGGCAGCAAGCACGTCACCGGCGCCCGTGTGGCGGCGATCGATGTGCAGGCGCACAAAGTCACCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCgg                                                                                                                                                 <  2:175004/137‑1 (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)
                     tGATCGAGGCCAAAGGCAGCAAGCACGTCACCGGCGCCCGTGTGGCGGCGATCGATGTGCAGGCGCACAAAGTCACCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCg                                                                                                                                                  >  4:75697/1‑123 (MQ=255)
                     tGATCGAGGCCAAAGGCAGCAAGCACGTCACCGGCGCCCGTGTGGCGGCGATCGATGTGCAGGCGCACAAAGTCACCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCg                                                                                                                                                  <  3:75697/123‑1 (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                                                                                                                              >  2:86814/1‑114 (MQ=255)
                                                  aCCGGCGCCCGTGTGGCGGCGATCGATGTGCAGGCGCACAAAGTCACCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCAc                                                                                                                              <  1:86814/114‑1 (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         >  7:297503/1‑149 (MQ=255)
                                                                                                                                    cctcctCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGCCGGGCGATGCGCCGCAGAAGCGTGAGTGCGTCGGTGGTATCAACGGCGTGTATGCCTTg         >  3:258896/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)
                                                                                                                                               |                                                                                                                                                
GCATACTCACCTCCAGCGCCGTGATCGAGGCCAAAGGCAGCAAGCACGTCACCGGCGCCCGTGTGGCGGCGATCGATGTGCAGGCGCACAAAGTCACCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGCCGGGCGATGCGCCGCAGAAGCGTGAGTGCGTCGGTGGTATCAACGGCGTGTATGCCTTGGGCGATG  >  NC_002947/390850‑391137

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