Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,685,132 A→T 6.2% intergenic (‑10/+222) PP_1480 ← / ← patD hypothetical protein/medium chain aldehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,685,1320AT 6.2% 189.4 / 3.0 65intergenic (‑10/+222)PP_1480/patDhypothetical protein/medium chain aldehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (27/34);  new base T (2/2);  total (29/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.73e-01

CCCGCGATCCACGGCCTTCGATCGCCTTGGCGCCTATCGATATCAGAAATGACTCCAAACGCGGCCACTCCAGCGACTTTGGAACCGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTATCTGTAGGAGTGGTCTTGC‑GAGCTAGGGGTGCTGGGCGTCGACGTTAGCCGTGCATGGCGGTAGTCGTGTTGTGTGGCGAAAATCGAGC  >  NC_002947/1684990‑1685266
                                                                                                                                              |                                                                                                                                       
accGCGATCCACGGCCTTCGATCGCCTTGGCGCCTATCGATATCAGAAATGACTCCAAACGCGGCCACTCCAGCGACTTTGGAACCGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGAt                                                                                                                                   <  7:139659/148‑1 (MQ=255)
          aCGGCCTTCGATCGCCTTGGCGCCTATCGATATCAGAAATGACTCCAAACGCGGCCACTCCAGCGACTTTGGAACCGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGAt                                                                                                                         >  8:235538/1‑149 (MQ=255)
          aCGGCCTTCGATCGCCTTGGCGCCTATCGATATCAGAAATGACTCCAAACGCGGCCACTCCAGCGACTTTGGAACCGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGAt                                                                                                                         >  4:98840/1‑149 (MQ=255)
          aCGGCCTTCGATCGCCTTGGCGCCTATCGATATCAGAAATGACTCCAAACGCGGCCACTCCAGCGACTTTGGAACCGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGAt                                                                                                                         >  4:190318/1‑149 (MQ=255)
              ccTTCGATCGCCTTGGCGCCTATCGATATCAGAAATGACTCCAAACGCGGCCACTCCAGCGACTTTGGAACCGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGAttgtt                                                                                                                     <  7:157375/149‑1 (MQ=255)
                ttCGATCGCCTTGGCGCCTATCGATATCAGAAATGACTCCAAACGCGGCCACTCCAGCGACTTTGGAACCGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCATATTCTGATACCat                                                                                                                              <  5:326320/138‑2 (MQ=255)
                ttCGATCGCCTTGGCGCCTATCGATATCAGAAATGACTCCAAACGCGGCCACTCCAGCGACTTTGGAACCGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCATATTCTGATACCat                                                                                                                              >  6:326320/1‑137 (MQ=255)
                ttCGATCGCCTTGGCGCCTATCGATATCAGAAATGACTCCAAACGCGGCCACTCCAGCGACTTTGGAACCGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGATTGTTgg                                                                                                                   >  2:224683/1‑149 (MQ=255)
                  cGATCGCCTTGGCGCCTATCGATATCGGAAATGACTCCAAACGCGGCCACTCCAGAGACTTTGGAACAGGCTTTTTGAAAATGGATTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGATTGTTGGt                                                                                                                  <  6:316140/148‑1 (MQ=255)
                   gATCGCCTTGGCGCCTATCGATATCAGAAATGACTCCAAACGCGGCCACTCCAGCGACTTTGGAACCGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGATTGTTGGTAt                                                                                                                <  7:33133/149‑1 (MQ=255)
                       gCCTTGGCGCATATCGATATCTGAAATGAGTCCAAACGCCGCCACTCCAGCGACTTTTGAACGGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCTTTTTATCAGAATATGATACCAGCTGATTGTTGGTATGTa                                                                                                             <  6:146863/148‑1 (MQ=255)
                        ccTTGGCGCCTATCGATATCAGAAATGACTCCAAACGCGGCCACTCCAGCGACTTTGGAACCGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAg                                                                                                                              >  7:164605/1‑130 (MQ=255)
                        ccTTGGCGCCTATCGATATCAGAAATGACTCCAAACGCGGCCACTCCAGCGACTTTGGAACCGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAg                                                                                                                              <  8:164605/130‑1 (MQ=255)
                            ggCGCCTATCGATATCAGAAATGACTCCAAACGCGGCCACTCCAGCGACTTTGGAACCGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGATTGTTGGTATGTAGGAt                                                                                                         <  1:17175/147‑1 (MQ=255)
                                 cTATCGATATCAGAAATGACTCCAAACGCGGCCACTCCAGCGACTTTGGAACCGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGa                                                                                                  <  7:191967/149‑1 (MQ=255)
                                   aTCGATATCAGAAATGACTCCAAACGCGGCCACTCCAGCGACTTTGGAACCGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGc                                                                                                >  5:117488/1‑149 (MQ=255)
                                          tCAGAAATGACTCCAAACGCGGCCACTCCAGCGACTTTGGAACCGGCTTTTTGAAAATGGCTTCAAAGGTCTATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTATc                                                                                          <  8:165534/148‑1 (MQ=255)
                                            aGACATGACTCCAAACGCGGCCACTCCAGCGACTTTGGAACCGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTATCTGt                                                                                       <  8:139026/149‑1 (MQ=255)
                                            aGAAATGACTCCAAACGCGGCCACTCCAGCGACTTTGGAACCGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTATCTg                                                                                        >  5:148174/1‑148 (MQ=255)
                                             gAAATGACTCCAAACGCGGCCACTCCAGCGACTTTGGAACCGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTATCTGTa                                                                                      <  6:117488/149‑1 (MQ=255)
                                                aTGACTCCAAACGCGGCCACTCCAGCGACTTTGGAACCGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTATCTGTAGGa                                                                                   <  3:98840/149‑1 (MQ=255)
                                                     tCCAAACGCGGCCACTCCAGCGACTTTGGAACCGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCt                                                                                                                            <  6:153699/103‑1 (MQ=255)
                                                     tCCAAACGCGGCCACTCCAGCGACTTTGGAACCGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCt                                                                                                                            >  5:153699/1‑103 (MQ=255)
                                                     tCCAAACGCGGCCACTCCAGCGACTTTGGAACCGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTATCTGTAGGAGTGGt                                                                              >  2:128538/1‑149 (MQ=255)
                                                       cAAACGCGGCCACTCCAGCGACTTTGGAACCGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTATCTGTAGGAGTGGTCt                                                                            >  3:298583/1‑149 (MQ=255)
                                                        aaaCGCGGCCACTCCAGCGACTTTGGAACCGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTATCTGTAGGAGTGGTCtt                                                                           <  7:235538/149‑1 (MQ=255)
                                                              ggCCACTCCAGCGACTTTGGAACCGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTATCTGTAGGAGTGGTCTTGC‑GAGc                                                                    <  6:148174/149‑1 (MQ=255)
                                                                   cTCCAGCGACTTTGGAACCGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTATCTGTAGGAGTGGTCTTGC‑GAGCTAggg                                                               <  1:128538/149‑1 (MQ=255)
                                                                    tCCAGCGACTTTGGAACCGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCa                                                                                                                               <  8:29555/85‑1 (MQ=255)
                                                                    tCCAGCGACTTTGGAACCGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCa                                                                                                                               >  7:29555/1‑85 (MQ=255)
                                                                      cAGCGACTTTGGAACCGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTATCTGTAGGAGTGGTCTTGC‑GAGCTAgggg                                                              <  8:230946/147‑1 (MQ=255)
                                                                      cAGCGACTTTGGAACCGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTATCTGTAGGAGTGGTCTTGC‑GAGCTAGGGGTg                                                            >  3:306442/1‑149 (MQ=255)
                                                                        gCGACTTTGGAACCGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTATCTGTAGGa                                                                                   <  8:3685/125‑1 (MQ=255)
                                                                        gCGACTTTGGAACCGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTATCTGTAGGa                                                                                   >  7:3685/1‑125 (MQ=255)
                                                                                ggAACCGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATTATTATTCATGGTATCAGAATATGATACCAGCTGATTGTTTGTATGTAGTATTGGCGGAGCCCTATCTGTAGGAGTGGTCTTGC‑GAGCTAGGGGTGCTGGGCGTc                                                   <  4:119001/148‑1 (MQ=255)
                                                                                  aaCCGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTATCTGTAGGAGTGGTCtt                                                                           >  1:109326/1‑123 (MQ=255)
                                                                                  aaCCGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTATCTGTAGGAGTGGTCtt                                                                           <  2:109326/123‑1 (MQ=255)
                                                                                     cGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGATTGTTGGTATGTAGGATGTGCGGAGCCCTATCTGTAGGAGTGGTCTTGC‑GAGCTAGGGGTGCTGGGCGTCGACGt                                              >  2:199558/1‑148 (MQ=255)
                                                                                     cGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTATCTGTAGGAGTGGTCTTGc                                                                         <  2:254351/122‑1 (MQ=255)
                                                                                     cGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTATCTGTAGGAGTGGTCTTGc                                                                         >  1:254351/1‑122 (MQ=255)
                                                                                     cGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTATCTGTAGGAGTGGTCTTGC‑GAGCTAGGGGTGCTGGGCGTCGACGt                                              >  2:128951/1‑148 (MQ=255)
                                                                                      ggCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTATCTGTAGGAGTGGTCTTGC‑GAGCTAGGGGTGCTGGGCGTCGACGTTa                                            <  3:190318/149‑1 (MQ=255)
                                                                                         tttttGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTATCTGTAGGAGTGGTCTTGC‑GAGCTAGGGGTGCTGGGCGTCGACGTTAGcc                                         >  2:207976/1‑149 (MQ=255)
                                                                                           tttGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTATCTGTAGGAGTGGTCTTGC‑GAGCTAGGGGTGCTGGGCGTCGACGTTAGCCGt                                       >  6:112546/1‑149 (MQ=255)
                                                                                              gAAAATGGCTTCAACGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTATCTGTAGGAGTGGTCTTGC‑GAGCTAGGGGTGCTGGGCGTCGACGTTAGCCGTGc                                     <  4:251460/148‑1 (MQ=255)
                                                                                              gAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTATCTGTAGGAGTGGTCTTGC‑GAGCTAGGGGTGCTGGGCGTCGACGTTAGCCGTGc                                     <  2:34602/148‑1 (MQ=255)
                                                                                                 aaTGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTATCTGTAGGAGTGGTCTTGC‑GAGCTAGGGGTGCTGGGCGTCGACGTTAGCCGTGCATgg                                 <  4:298583/149‑1 (MQ=255)
                                                                                                    ggCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTATCTGTAGGAGTGGTCTTGC‑GAGCTAGGGGTGCTGGGCGTCGACGTTAGCCGTGCATGGCg                               >  6:209629/1‑148 (MQ=255)
                                                                                                    ggCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTATCTGTAGGAGTGGTCTTGC‑GAGCTAGGGGTGCTGGGCGTCGACGTTAGCCGTGCATGGCg                               <  5:209629/148‑1 (MQ=255)
                                                                                                       ttCAAAGGTCGATCTTTGATGATTATTCATGATATCAGAATATGATACCAGCTGATTGTTGGTATGTa                                                                                                             >  8:72546/1‑68 (MQ=255)
                                                                                                       ttCAAAGGTCGATCTTTGATGATTATTCATGATATCAGAATATGATACCAGCTGATTGTTGGTATGTa                                                                                                             <  7:72546/68‑1 (MQ=255)
                                                                                                         cAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTATCTGTAGGAGTGGTCTTGCAGAGCTAGGGGTGCTGGGCGTCGACGTTAGCCGTGCATGGCGGTa                            <  1:224683/147‑1 (MQ=255)
                                                                                                         cAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTATCTGTAGGAGTGGTCTTGC‑GAGCTAGGGGTGCTGGGCGTCGACGTTAGCCGTGCATGGCGGTAGTc                         >  5:5897/1‑149 (MQ=255)
                                                                                                                 gATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTATCTGTAGGAGTGGTCTTGC‑GAGCTAGGGGTGCTGGGCGTCGACGGTAGCCGTGCATGGCGGTAGTCGtgtt                    >  3:244690/1‑146 (MQ=255)
                                                                                                                 gATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTATCTGTAGGAGTGGTCTTGC‑GAGCTAGGGGTGCTGGGCGTCGACGGTAGCCGTGCATGGCGGTAGTCGtgtt                    <  4:244690/146‑1 (MQ=255)
                                                                                                                 gATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGATTGTTGGTATGCAGGATGGGCGGAGCCCTATCTGTAGGAGTGGTCTTGC‑GAGCTAGGGGTGCTGGGCGTCGACGTTAGCCGTGCATGGCGGTAGTCGTGTtgtg                 <  7:215611/149‑1 (MQ=255)
                                                                                                                       tgatgaTTATTCATGGTATCAGAATATGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTATCTGTAGGAGTGGTCTTGC‑GAGCTAGGGGTGCTGGGCGTCGACGTTAGCCGTGCATGGCGGTAGTCgtgt                     <  4:175978/139‑1 (MQ=255)
                                                                                                                       tgatgaTTATTCATGGTATCAGAATATGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTATCTGTAGGAGTGGTCTTGC‑GAGCTAGGGGTGCTGGGCGTCGACGTTAGCCGTGCATGGCGGTAGTCgtgt                     >  3:175978/1‑139 (MQ=255)
                                                                                                                       tgatgaTTATTCATGGTATCAGAATATGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTATCTGTAGGAGTGGTCTTGC‑GAGCTAGGGGTGCTGGGCGTCGACGTTAGCCGTGCATGGCGGTAGTCGTGTTGTGTgg              <  1:207976/146‑1 (MQ=255)
                                                                                                                           gatTATTCATGGTATCAGAATATGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTATCTGTAGGAGTGGTCTTGC‑GAGCTAGGGGTGCTGGGCGTCGACGTTAGCCGTGCATGGCGGTAGTCGTGTTGTGTGGCGAAAAt       >  1:8036/1‑149 (MQ=255)
                                                                                                                                ttCATGGTATCAGAATATGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTATCTGTAGGAGTGGTCTTGC‑GAGCTAGGGGTGCTGGGCGTCGACGTTAGCCGTGCATGGCGGTAGTCGTGTTGTGTGGCGAAAATCGAgc  >  7:240283/1‑149 (MQ=255)
                                                                                                                                  gctGGTATCATATTCTGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTATCTGTAGGAGTGGTCTTGC‑GAGCTAGGGGTGCTGGGCGTCGACGTTAGCCGTg                                      >  3:77593/3‑111 (MQ=255)
                                                                                                                                  gctGGTATCATATTCTGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTATCTGTAGGAGTGGTCTTGC‑GAGCTAGGGGTGCTGGGCGTCGACGTTAGCCGTg                                      <  4:77593/109‑1 (MQ=255)
                                                                                                                                         tCAGAATATGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTATCTGTAGGAGTGGTCTTGC‑GAGCTAGGGGTGCTGGGCGTAGACGTTAGCCGTGCATGGCGGTAGTCGTGTTGTGTGGCGAAAATCGa    <  6:203333/138‑1 (MQ=255)
                                                                                                                                         tCAGAATATGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTATCTGTAGGAGTGGTCTTGC‑GAGCTAGGGGTGCTGGGCGTAGACGTTAGCCGTGCATGGCGGTAGTCGTGTTGTGTGGCGAAAATCGa    >  5:203333/1‑138 (MQ=255)
                                                                                                                                              |                                                                                                                                       
CCCGCGATCCACGGCCTTCGATCGCCTTGGCGCCTATCGATATCAGAAATGACTCCAAACGCGGCCACTCCAGCGACTTTGGAACCGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTATCTGTAGGAGTGGTCTTGC‑GAGCTAGGGGTGCTGGGCGTCGACGTTAGCCGTGCATGGCGGTAGTCGTGTTGTGTGGCGAAAATCGAGC  >  NC_002947/1684990‑1685266

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑