Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,246,185 C→G 10.5% M30I (ATG→ATC)  alr ← alanine racemase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,246,1850CG10.5% 153.7 / 10.2 57M30I (ATG→ATC) alralanine racemase
Reads supporting (aligned to +/- strand):  ref base C (24/27);  new base G (4/2);  total (28/29)
Fisher's exact test for biased strand distribution p-value = 4.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.48e-01

CATAGGCATCGGCCTTGAGCACGGCGCACAGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTAACAGACGCAGGTGTGTTTGA  >  NC_002947/4246037‑4246331
                                                                                                                                                    |                                                                                                                                                  
cATAGGCATCGGCCGTGAGCACGGCGCACAGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTcc                                                                                                                                                    <  6:327727/149‑1 (MQ=255)
   aGGCATCGGCCTTGAGCACGGCGCACAGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATc                                                                                                                                                 >  4:190366/1‑149 (MQ=255)
       aTCGGCCTTGAGCACGGCGCACAGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGaa                                                                                                                                              >  3:320421/1‑147 (MQ=255)
       aTCGGCCTTGAGCACGGCGCACAGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGAc                                                                                                                                              >  7:197019/1‑148 (MQ=255)
            ccTTGAGCACGGCGCACAGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCGATGGACAACggc                                                                                                                                        <  5:163191/149‑2 (MQ=255)
                gAGCACGGCGCACAGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCGATGGACAACggcac                                                                                                                                      <  4:168358/147‑4 (MQ=255)
                gAGCACGGCGCACAGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGAc                                                                                                                                              >  6:386988/1‑139 (MQ=255)
                gAGCACGGCGCACAGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGAc                                                                                                                                              <  5:386988/139‑1 (MQ=255)
                gAGCACGGCGCACAGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACggtgg                                                                                                                                      <  5:328988/147‑1 (MQ=255)
                gAGCACGGCGCACAGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGgtggtg                                                                                                                                    >  1:346490/1‑149 (MQ=255)
                  gCACGGCGCACAGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGcc                                                                                                                                  <  6:271855/149‑1 (MQ=255)
                       gcgcACAGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATa                                                                                                                             <  1:372287/149‑1 (MQ=255)
                            caGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAggg                                                                                                                        <  1:156376/149‑1 (MQ=255)
                             aGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAggg                                                                                                                        <  4:57534/148‑1 (MQ=255)
                             aGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAggg                                                                                                                        >  3:57534/1‑148 (MQ=255)
                                  ggACTTGCCGGCCAGCTCGGCCTGCAGCGTGCTGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGc                                                                                                                    <  7:78780/147‑1 (MQ=255)
                                    aCTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGc                                                                                                                    <  5:348555/145‑1 (MQ=255)
                                    aCTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGc                                                                                                                    >  6:348555/1‑145 (MQ=255)
                                     cTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTc                                                                                                               >  8:93821/1‑149 (MQ=255)
                                         ccGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGt                                                                                                           <  7:293945/149‑1 (MQ=255)
                                          cGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTg                                                                                                          <  1:33583/149‑1 (MQ=255)
                                          cGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTg                                                                                                          <  3:190366/149‑1 (MQ=255)
                                            gCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGa                                                                                                         >  6:349799/1‑148 (MQ=255)
                                             ccAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATc                                                                                                       <  5:349799/149‑1 (MQ=255)
                                              cAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCa                                                                                                      <  4:221920/149‑1 (MQ=255)
                                                 ctgggCCTGCAGCGTGCGGATGTGGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTTCCGTTTTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGaa                                                                                                   <  6:352034/146‑1 (MQ=255)
                                                 cTCGGCCTGCAGCGTGCGGATGTGGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGaa                                                                                                   <  4:134017/149‑1 (MQ=255)
                                                  tCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCAGCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAg                                                                                                  >  8:107261/1‑149 (MQ=255)
                                                    ggCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCACCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAg                                                                                                  <  1:285672/147‑1 (MQ=255)
                                                              gTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGc                                                                                      <  3:116106/149‑1 (MQ=255)
                                                                gCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACA‑GGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGc                                                                                   <  7:6209/149‑1 (MQ=255)
                                                                     tgttgtGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAg                                                                                >  5:221813/1‑148 (MQ=255)
                                                                        tgtgCTGCAGAGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAg                                                                             <  7:134700/148‑1 (MQ=255)
                                                                         gtgCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAggg                                                                           >  3:232257/1‑149 (MQ=255)
                                                                          tgctgcAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAggg                                                                           <  7:94951/148‑1 (MQ=255)
                                                                            ctgcAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTa                                                                         <  1:87055/148‑1 (MQ=255)
                                                                                     gctggcgctgACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAgg                                                               >  7:28890/1‑149 (MQ=255)
                                                                                     gctggcgctgACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAgg                                                               >  1:89379/1‑149 (MQ=255)
                                                                                     gctggcgctgACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAgg                                                               <  3:240729/149‑1 (MQ=255)
                                                                                           gctgACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATgg                                                         >  8:313393/1‑149 (MQ=255)
                                                                                                      cccAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCgt                                              >  1:380985/1‑149 (MQ=255)
                                                                                                      cccAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCgt                                              >  7:113134/1‑149 (MQ=255)
                                                                                                      cccAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCgt                                              >  4:306017/1‑149 (MQ=255)
                                                                                                       ccAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCgtg                                             >  6:304054/1‑149 (MQ=255)
                                                                                                       ccAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCgtg                                             >  2:44150/1‑149 (MQ=255)
                                                                                                             aTTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTg                                       >  8:42267/1‑149 (MQ=255)
                                                                                                                gCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTc                                     <  6:221813/148‑1 (MQ=255)
                                                                                                                                aGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTaa                     <  7:93821/148‑1 (MQ=255)
                                                                                                                                 gggTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTAACa                   <  1:218866/149‑1 (MQ=255)
                                                                                                                                        ggTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTAACAGACGCAg            >  4:318165/1‑149 (MQ=255)
                                                                                                                                         gTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTAACAGACGCAgg           >  2:162268/1‑149 (MQ=255)
                                                                                                                                          caccGTTGTCGATGGACAACGGTGGTGCCGCCTACAGGGGGGCCGGTCCGGTGATCAGCAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGATAGGGCATGGTGTTTACTCGTGTTCTTGTCGAAGTGCGCGCCTTAACAGACGCAgct          >  6:110593/3‑147 (MQ=255)
                                                                                                                                          caccGTTGTCGATGGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTAACAGACGCAGgt          >  8:259916/3‑149 (MQ=255)
                                                                                                                                           accGTTGTCGATGGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTAACAGACGCAGgtg         >  3:284437/2‑149 (MQ=255)
                                                                                                                                           accGTTGTCGATGGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCAAGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTAACAGACGCAGgtg         >  6:336422/2‑149 (MQ=255)
                                                                                                                                              gTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTAACAGACGCAGgtgtgt      >  6:68557/1‑149 (MQ=255)
                                                                                                                                              gTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTAACAGACGCAGgtgtgt      <  1:44150/149‑1 (MQ=255)
                                                                                                                                                  tCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTAACAGACGCAGGTGTGTTTGa  <  2:380985/149‑1 (MQ=255)
                                                                                                                                                    cATCGACAACGGTGGTGCCGCATACAGGGGGGCCAGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTAGTAGTGTTCTTGTCGTAGTGCGAGCCTTAACAGACGCAGGTGTGttt    >  6:314018/1‑145 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CATAGGCATCGGCCTTGAGCACGGCGCACAGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTAACAGACGCAGGTGTGTTTGA  >  NC_002947/4246037‑4246331

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