Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,170,185 C→A 11.5% intergenic (‑293/‑94) katG ← / → PP_3669 catalase‑peroxidase/LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,170,1850CA11.5% 168.0 / 11.3 61intergenic (‑293/‑94)katG/PP_3669catalase‑peroxidase/LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (27/27);  new base A (3/4);  total (30/31)
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.09e-01

GGTGCTGACGCGCTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACCTGCTGGTGCTGTTCGAA  >  NC_002947/4170038‑4170332
                                                                                                                                                   |                                                                                                                                                   
ggTGCTGACGCGCTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCCGGCCTGTATCGGCCTCATCGCGGTCAAGcc                                                                                                                                                    <  4:84775/149‑1 (MQ=255)
 gTGCTGACGCGCTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGccc                                                                                                                                                   >  5:267670/1‑149 (MQ=255)
       aCGCGCTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGGCTTGACCGCGAt                                                                                                                                              >  8:92636/1‑148 (MQ=255)
         gcgcTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAg                                                                                                                                           <  5:110912/149‑1 (MQ=255)
         gcgcTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGGCTTGACCGCGATGAg                                                                                                                                           <  5:151440/149‑1 (MQ=255)
         gcgcTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGGCTTGACCGCGATGAg                                                                                                                                           <  7:92636/149‑1 (MQ=255)
         gcgcTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGGCTTGACCGCGATGAg                                                                                                                                           <  3:196180/149‑1 (MQ=255)
             tGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGCCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCg                                                                                                                                       >  6:277574/1‑149 (MQ=255)
               cccTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGc                                                                                                                                     >  7:120670/1‑149 (MQ=255)
                ccTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGca                                                                                                                                    <  7:172451/149‑1 (MQ=255)
                 cTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGca                                                                                                                                    >  4:101194/1‑148 (MQ=255)
                    tCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAg                                                                                                                                           <  3:138239/138‑1 (MQ=255)
                    tCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAg                                                                                                                                           >  4:138239/1‑138 (MQ=255)
                     cgcgGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATcc                                                                                                                               <  5:89860/149‑1 (MQ=255)
                      gcgGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCa                                                                                                                              >  7:40389/1‑149 (MQ=255)
                               gcgGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATacaaggca                                                                                                                     >  8:165938/1‑149 (MQ=255)
                                 ggggCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATacaaggcac                                                                                                                    >  8:315154/1‑148 (MQ=255)
                                   ggCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATacaaggcacaag                                                                                                                 <  8:40389/149‑1 (MQ=255)
                                        cGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGcccc                                                                                                            >  1:154845/1‑149 (MQ=255)
                                           gaAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCcgcg                                                                                                         <  3:101194/149‑1 (MQ=255)
                                                    gCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCtt                                                                                                 <  7:16982/148‑1 (MQ=255)
                                                    gCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCtt                                                                                                 <  7:165938/148‑1 (MQ=255)
                                                    gCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGCCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCtt                                                                                                 <  5:277574/148‑1 (MQ=255)
                                                        aGCTTGGGCAGATGTCTAATCGAGTGGCTGTTGGAGCCTGATAGGGTGGCTCTGTCAGGCGCAGGCCTGAATCGCCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATa                                                                                            <  7:178470/149‑1 (MQ=255)
                                                              ggCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAg                                                                                       <  6:267670/148‑1 (MQ=255)
                                                              ggCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCCCGCCATGCTTCAATACACAGc                                                                                      >  6:317650/1‑149 (MQ=255)
                                                                          aTAGAGTGGCTGTTGTAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCg                                                                                                                                       <  6:306629/88‑1 (MQ=255)
                                                                          aTAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCg                                                                                                                                       >  5:306629/1‑88 (MQ=255)
                                                                                  gCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACg                                                                  >  6:149400/1‑149 (MQ=255)
                                                                                  gCTGTTGGAGCCTGATAGGCTGGCTATGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGCGGCCGGCACATCCATACAAGGCACAAGCCCCCGCGGCCTGCTTCAATACACAGCCATACCCTTTGCCCAGGACg                                                                  >  8:208627/1‑149 (MQ=255)
                                                                                   cTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGc                                                                 >  8:32630/1‑149 (MQ=255)
                                                                                    tGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGcc                                                                <  7:315154/149‑1 (MQ=255)
                                                                                          aGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGccc                                                                        <  6:140018/135‑1 (MQ=255)
                                                                                          aGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGccc                                                                        >  5:140018/1‑135 (MQ=255)
                                                                                                aTAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCAt                                                      <  2:2330/147‑1 (MQ=255)
                                                                                                        gctCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCACGAACCgagaa                                            >  5:89900/1‑145 (MQ=255)
                                                                                                        gctCTGGCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCACGAAcc                                                 <  6:89900/144‑1 (MQ=255)
                                                                                                              tCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTg                                      >  8:4466/1‑149 (MQ=255)
                                                                                                                aGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTgc                                     <  5:149400/148‑1 (MQ=255)
                                                                                                                    gcAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGcgccgc                                >  2:89156/1‑149 (MQ=255)
                                                                                                                        gacTGTTTCGTCCTAATCGCGGTCACGCCTGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTc                             <  3:32983/146‑1 (MQ=255)
                                                                                                                        gCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTc                             >  6:276932/1‑148 (MQ=255)
                                                                                                                               tCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAAc                                            >  2:80203/1‑126 (MQ=255)
                                                                                                                               tCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAAc                                            <  1:80203/126‑1 (MQ=255)
                                                                                                                                   ccTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGc                                                                 >  3:57877/1‑101 (MQ=255)
                                                                                                                                   ccTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGc                                                                 <  4:57877/101‑1 (MQ=255)
                                                                                                                                   ccTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGa                           >  4:106087/1‑139 (MQ=255)
                                                                                                                                   ccTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGa                           <  3:106087/139‑1 (MQ=255)
                                                                                                                                    cTCATCGCGGGCTTGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGCCCTGCGCCGCGTCGACATGAACctgc                >  5:256965/1‑149 (MQ=255)
                                                                                                                                     tCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACctgct               >  8:77727/1‑149 (MQ=255)
                                                                                                                                     tCATCGCGGGCTTGACCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTgc                                     >  6:13858/1‑127 (MQ=255)
                                                                                                                                       aTCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACctgctg              <  5:317650/148‑1 (MQ=255)
                                                                                                                                       aTCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACCTGCTgg             >  3:206304/1‑149 (MQ=255)
                                                                                                                                       aTCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACCTGCTgg             >  8:293035/1‑149 (MQ=255)
                                                                                                                                       aTCGCGGGCTTGACCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGc                                               <  4:211990/115‑1 (MQ=255)
                                                                                                                                       aTCGCGGGCTTGACCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGc                                               >  3:211990/1‑115 (MQ=255)
                                                                                                                                        tCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACCTGCTGGt            <  7:32630/149‑1 (MQ=255)
                                                                                                                                         cgcgGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACCTGCTgg             >  2:175245/1‑147 (MQ=255)
                                                                                                                                         cgcgGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACCTGCTgg             <  1:175245/147‑1 (MQ=255)
                                                                                                                                             gTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACCTGCTGGTGCTGt       <  7:293035/149‑1 (MQ=255)
                                                                                                                                              tCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGa                           <  2:127426/128‑1 (MQ=255)
                                                                                                                                              tCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGa                           >  1:127426/1‑128 (MQ=255)
                                                                                                                                                  gCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACCTGCTGGTGCTGTTCGaa  <  7:77727/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GGTGCTGACGCGCTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACCTGCTGGTGCTGTTCGAA  >  NC_002947/4170038‑4170332

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