Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,532,634 C→A 9.5% W183C (TGG→TGT)  PP_3121 ← transcriptional regulator, LysR family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,532,6340CA 9.5% 185.4 / 7.0 61W183C (TGG→TGT) PP_3121transcriptional regulator, LysR family
Reads supporting (aligned to +/- strand):  ref base C (26/29);  new base A (3/3);  total (29/32)
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.98e-01

GCCACCATGCGGCCGATGGTCGTCAGCTGGTGGCTTTCGAACGCGACTGCCAGCTTGCCGTGGCGGGCCGCCACATTCTGCTCCATCAGCAGGCGCACCGCCGAAGGGCGCTGCAGGGCCACGAAATCTTCGGCCATCAGCTGTGCCCAGCTGACCTGGGGCAGCAGTGCCAAGGGCGAGTCGGCGGCCACTACGGCAACGAAGCGGTCCATGTACAACGGGTGGAAGTCCAGCCCATCGATGTTATCCGGCTCGAAGCCGATGCCCAGTTCGACGCGGCGATGGCGTACCAGTT  >  NC_002947/3532488‑3532782
                                                                                                                                                  |                                                                                                                                                    
gCCACCATGCGGCCGATGGTCGTCAGCTGGTGGCTTTCGAACGCGACTGCCAGCTTGCCGTGGCGGGCCGCCACATTCTGCTCCATCAGCAGGCGCACCGCCGAAGGGCGCTGCAGGGCCACGAAATCTTCGGCCATCAGCTGTGcc                                                                                                                                                      <  2:86766/147‑1 (MQ=255)
      aTGCGGCCGATGGTCGTCAGCTGGTGGCTTTCGAACGCGACTGCCAGCTTGCCGTGGCGGGCCGCCACATTCTGCTCCATCAGCAGGCGCACCGCCGAAAGGCGCTGCAGGGCCACGAAATCTTCGGCCATCAGCTGTGCCCAGCTGAc                                                                                                                                              <  1:193465/149‑1 (MQ=255)
         cGGCCGATGGTCGTCAGCTGGTGGCTTTCGAACGCGACTGCCAGCTTGCCGTGGCGGGCCGCCACATTCTGCTCCATCAGCAGGCGCACCGCCGAAGGGCGCTGCAGGGCCACGAAATCTTCGGCCATCAGCTGTGCCCAGCTGACCTg                                                                                                                                           <  7:157865/149‑1 (MQ=255)
         cGGCCGATGGTCGTCAGCTGGTGGCTTTCGAACGCGACTGCCAGCTTGCCGTGGCGGGCCGCCACATTCTGCTCCATCAGCAGGCGCACCGCCGAAGGGCGCTGCAGGGCCACGAAATCTTCGGCCATCAGCTGGGCACAGCTGatgg                                                                                                                                            <  1:40681/148‑4 (MQ=255)
                     gtcAGCTGGTGGCTTTCGAACGCGACTGCCAGCTTGCCGTGGCGGGCCGCCACATTCTGCTCCATCAGCAGGCGCACCGCCGAAGGGCGCTGCAGGGCCACGAAATCTTCGGCCATCAGCTGGGCACAGCTGatg                                                                                                                                             <  4:325801/135‑3 (MQ=255)
                     gtcAGCTGGTGGCTTTCGAACGCGACTGCCAGCTTGCCGTGGCGGGCCGCCACATTCTGCTCCATCAGCAGGCGCACCGCCGAAGGGCGCTGCAGGGCCACGAAATCTTCGGCCATCAGCTGGGCACAGCTGatg                                                                                                                                             >  3:325801/1‑133 (MQ=255)
                        aGCTGGTGGCTTTCGAACGCGACTGCCAGCTTGCCGTGGCGGGCCGCCACATTCTGCTCCATCAGCAGGCGCACCGCCGAAGGGCGCTGCAGGGCCACGAAATCTTCGGCCATCAGCTGTGCCCAGCTGACCTGGGGCAGCAGTGcc                                                                                                                              >  1:257747/1‑147 (MQ=255)
                        aGCTGGTGGCTTTCGAACGCGACTGCCAGCTTGCCGTGGCGGGCCGCCACATTCTGCTCCATCAGCAGGCGCACCGCCGAAGGGCGCTGCAGGGCCACGAAATCTTCGGCCATCAGCTGTGCCCAGCTGACCTGGGGCAGCAGTGcc                                                                                                                              >  3:248414/1‑147 (MQ=255)
                        aGCTGGTGGCTTTCGAACGCGACTGCCAGCTTGCCGTGGCGGGCCGCCACATTCTGCTCCATCAGCAGGCGCACCGCCGAAGGGCGCTGCAGGGCCACGAAATCTTCGGCCATCAGCTGTGCCCAGCTGACCTGGGGCAGCAGTGcc                                                                                                                              >  3:231983/1‑147 (MQ=255)
                           tggtggCTTTCGAACGCGACTGCCAGCTTGCCGTGGCGGGCCGCCACATTCTGCTCCATCAGCAGGCGCACCGCCGAAGGGCGCTGCAGGGCCACGAAATCTTCGGCCATCAGCTGTGCCCAGCTGACCTGGGGCAGCAGTGCCAAggg                                                                                                                         >  3:70475/1‑149 (MQ=255)
                               gtcTTTCGAACGAGACTGCCATCTTGACGTGGCGGGCCGCCACATTCTGATCACTCAGCTGGCGCACCGCCGAAGGGCGCTGCAGGGCCACGAAATCTTCGGCCATCAGCTGTGCCCAGCTGACCTg                                                                                                                                           <  1:160834/125‑1 (MQ=255)
                               ggCTTTCGAACGAGACTGCCAGATTGCCGTGGCGGGCCGCCACCTTCTGCTCCATCAGCAGGCGCACCGCCGAAGGGCGCTGCAGGGAAACGAAATCTTCGGCCATCAACTGTGCCAATCTGAccag                                                                                                                                           >  2:160834/1‑125 (MQ=255)
                                     cGAACGCGACTGCCAGCTTGCCGTGGCGGGCCGCCACATTCTGCTCCATCAGCAGGCGCACCGCCGAAGGGCGCTGCAGGGCCACGAAATCTTCGGCCATCAGCTGTGCCCAGCTGACCTGGGGCAGCAGTGCCAAGGGCGAGTcgac                                                                                                                <  8:6867/148‑3 (MQ=255)
                                      gAACGCGACTGCCAGCTTGCCGTGGCGGGCCGCCACATTCTGCTCCATCAGCAGGCGCACCGCCGAAGGGCGCTGCAGGGCCACGAAATCTTCGGCCATCAGCTGTGCCCAGCTGACCTGGGGCAGCAGTGCCAAGGGCGAGTcggcgg                                                                                                              >  3:198937/1‑149 (MQ=255)
                                      gAACGCGACTGCCAGCTTGCCGTGGCGGGCCGCCACATTCTGCTCCATCAGCAGGCGCACCGCCGAAGGGCGCTGCAGGGCCACGAAATCTTCGGCCATCAGCTGTGCCCAGCTGACCTGGGGCAGCAGTGCCAAGGGCGAGTcggcg                                                                                                               <  1:17505/148‑1 (MQ=255)
                                      gAACGCGACTGCCAGCTTGCCGTGGCGGGCCGCCACATTCTGCTCCATCAGCAGGCGCACCGCCGAAGGGCGCTGCAGGGCCACGAAATCTTCGGCCATCAGCTGTGCCCAGCTGACCTGGGGCAGCAGTGCCAAGGGCGAGTcggcg                                                                                                               <  1:181459/148‑1 (MQ=255)
                                              cTGCCAGCTTGCCGTGGCGGGCCGCCACATTCTGCTCCATCAGCAGGCGCACCGCCGAAGGGCGCTGCAGGGCCACGAAATCTTCGGCCATCAGCTGTGCCCAGCTGACCTGGGGCAGCAGTGCCAAGGGCGAGTCGGCGGCCACTACg                                                                                                      <  4:248414/149‑1 (MQ=255)
                                              cTGCCAGCTTGCCGTGGCGGGCCGCCACATTCTGCTCCATCAGCAGGCGCACCGCCGAAGGGCGCTGCAGGGCCACGAAATCTTCGGCCATCAGCTGTGCCCAGCTGACCTGGGGCAGCAGTGCCAAGGGCGAGTCGGCGGCCACTACg                                                                                                      <  4:231983/149‑1 (MQ=255)
                                              cTGCCAGCTTGCCGTGGCGGGCCGCCACATTCTGCCCCATCAGCAGGCGCACCGCCGAAGGGCGCTGCAGGGCCACGAAATCTTCGGCCATCAGCTGTGCCCAGCTGACCTGGGGCAGCAGTGCCAAGGGCGAGTCGGCGGCCACTACg                                                                                                      <  2:163932/149‑1 (MQ=255)
                                               tGCCAGCTTGCCGTGGCGGGCCGCCACATTCTGCTCCATCAGCAGGCGCACCGCCGAAGGGCGCTGCAGGGCCACGAAATCTTCGGCCATCAGCTGTGCCCAGCTGACCTGGGGCAGCAGTGCCAAGGGCGAGTCGGCGGCCACTACgg                                                                                                     >  2:2337/1‑149 (MQ=255)
                                                        gCCGTGGCGGGCCGCCACATTCTGCTCCATCAGCAGGCGCACCGCCGAAGGGCGCTGCAGGGCCACGAAATCTTCGGCCATCAGCTGGGCACAGCTGatgg                                                                                                                                            <  6:364068/101‑4 (MQ=255)
                                                        gCCGTGGCGGGCCGCCACATTCTGCTCCATCAGCAGGCGCACCGCCGAAGGGCGCTGCAGGGCCACGAAATCTTCGGCCATCAGCTGGGCACAGCTGatgg                                                                                                                                            >  5:364068/1‑98 (MQ=255)
                                                         ccGTGGCGGGCCGCCACATTCTGCTCCATCAGCAGGCGCACCGCCGAAGGGCGCTGCAGGGCCACGAAATCTTCGGCCATCAGCTGTGCCCAGCTGACCTGGGGCAGCAGTGCCAAGGGCGAGTCGGCGGCCACTACGGCAACGAAGc                                                                                            <  2:257747/148‑1 (MQ=255)
                                                         ccGTGGCGGGCCGCCACATTCTGCTCCATCAGCAGGCGCACCGCCGAAGGGCGCTGCAGGGCCACGAAATCTTCGGCCATCAGCTGTGCCCAGCTGACCTGGGGCAGCAGTGCCAAGGGCGAGTCGGCGGCCACTACGGCAACGAAGc                                                                                            <  1:2337/148‑1 (MQ=255)
                                                         ccGTGGCGGGCCGCCACATTCTGCTCCATCAGCAGGCGCACCGCCGAAGGGCGCTGCAGGGCCACGAAATCTTCGGCCATCAGCTGTGCCCAGCTGACCTGGGGCAGCAGTGCCAAGGGCGAGTCGGCGGCCACTACGGCAACGAAGCg                                                                                           >  2:331495/1‑149 (MQ=255)
                                                          cGTGGCGGGCCGCCACATTCTGCTCCATCAGCAGGCGCACCGCCGAAGGGCGCTGCAGGGCCACGAAATCTTCGGCCATCAGCTGTGCCCAGCTGAc                                                                                                                                              <  8:123330/97‑1 (MQ=255)
                                                          cGTGGCGGGCCGCCACATTCTGCTCCATCAGCAGGCGCACCGCCGAAGGGCGCTGCAGGGCCACGAAATCTTCGGCCATCAGCTGTGCCCAGCTGAc                                                                                                                                              >  7:123330/1‑97 (MQ=255)
                                                                gggCCGCCACATTCTGCTCCATCAGCAGGCGCACCGCCGAAGGGCGCTGCAGGGCCACGAAATCTTCGGCCATCAGCTGTGCCCAGCTGACCt                                                                                                                                            <  2:315562/93‑1 (MQ=255)
                                                                gggCCGCCACATTCTGCTCCATCAGCAGGCGCACCGCCGAAGGGCGCTGCAGGGCCACGAAATCTTCGGCCATCAGCTGTGCCCAGCTGACCt                                                                                                                                            >  1:315562/1‑93 (MQ=255)
                                                                   ccgccACATTCTGCTCCATCAGCAGGCGCACCGCCGAAGGGCGCTGCAGGGCCACGAAATCTTCGGCCATCAGCTGTGCCCAGCTGACCTGGGGCAGCAGTGCCAAGGGCGAGTCGGCGGCCACTACGGCAACGAAGCAGTCCATGTa                                                                                  >  4:115685/1‑148 (MQ=255)
                                                                    cgccACATTCTGCTCCATCAGCAGGCGCACCGCCGAAGGGCGCTGCAGGGCCACGAAATCTTCGGCCATCAGCTGTGCCCAGCTGACCTGGGGCAGCAGTGCCAAGGGCGAGTCGGCGGCCACTACGGCAACGAAGCGGTCCATg                                                                                    >  6:296434/1‑145 (MQ=255)
                                                                    cgccACATTCTGCTCCATCAGCAGGCGCACCGCCGAAGGGCGCTGCAGGGCCACGAAATCTTCGGCCATCAGCTGTGCCCAGCTGACCTGGGGCAGCAGTGCCAAGGGCGAGTCGGCGGCCACTACGGCAACGAAGCGGTCCATg                                                                                    <  5:296434/145‑1 (MQ=255)
                                                                          aTTCTGCTCCATCAGCAGGCGCACCGCCGAAGGGCGCTGCAGGGCCACGAAATCTTCGGCCATCAGCTGTGCCCAGCTGACCTGGGGCAGCAGTGCCAAGGGCGAGTCGGCGGCCACTACGGCAACGAAGCGGTCCATGTACAACGGGt                                                                          <  4:198937/149‑1 (MQ=255)
                                                                          aTTCTGCTCCATCAGCAGGCGCACCGCCGAAGGGCGCTGCAGGGCCACGAAATCTTCGGCCATCAGCGGTGCCCAGCTGACCTGGGGCAGCAGTGCCAAGGGCGAGTCGGCGGCCACTACGGCAACGAAGCGGTCCATGTACAACGGGt                                                                          <  4:283967/149‑1 (MQ=255)
                                                                               gCTCCATCAGCAGGCGCACCGCCGAAGGGCGCTGCAGGGCCACGAAATCTTCGGCCATCAGCTGTGCCCAGCTGACCTGGGGCAGCAGTGCCAAGGGCGAGTCGGCGGCCACTACGGCAACGAAg                                                                                             <  3:228549/125‑1 (MQ=255)
                                                                               gCTCCATCAGCAGGCGCACCGCCGAAGGGCGCTGCAGGGCCACGAAATCTTCGGCCATCAGCTGTGCCCAGCTGACCTGGGGCAGCAGTGCCAAGGGCGAGTCGGCGGCCACTACGGCAACGAAg                                                                                             >  4:228549/1‑125 (MQ=255)
                                                                                cTCCATCAGCAGGCGCACCGCCGAAGGGCGCTGCAGGGCCACGAAATCTTCGGCCATCAGCTGTGCCCAGCTGACCTGGGGCAGCAGTGCCAAGGGCGAGTCGGCGGCCACTACGGCAACGAAGCAGTCCATGTACAACGGGTGGaa                                                                      <  3:115685/147‑1 (MQ=255)
                                                                                   cATCAGCAGGCGCACCGCCGAAGGGCGCTGCAGGGCCTCGAAATCTTCGGCCATCAGCTGTGCCCAGGTGACCTGGGGCAGCAGTGCCAAGGGCGAGTCGGCGGCCACTACGGCAACGAAGCGGTCCATGTACAACGGGTGGAAGTCCa                                                                 <  5:143743/149‑1 (MQ=255)
                                                                                   cATCAGCAGGCGCACCGCCGAAGGGCGCTGCAGGGCCACGAAATCTTCGGCCATCAGCTGTGCCCAGCTGACCTGGGGCAGCAGTGCCAAg                                                                                                                           <  5:356047/91‑1 (MQ=255)
                                                                                   cATCAGCAGGCGCACCGCCGAAGGGCGCTGCAGGGCCACGAAATCTTCGGCCATCAGCTGTGCCCAGCTGACCTGGGGCAGCAGTGCCAAg                                                                                                                           >  6:356047/1‑91 (MQ=255)
                                                                                   cATCAGCAGGCGCACCGCCGAAGGGCGCTGCAGGGCCACGAAATCTTCGGCCATCAGCTGTGCCCAGCTGACCTGGGGCAGCAGTGCCAAGGGCGAGTCGGCGGCCACTACGGCAACGAAGCGGTCCATGTa                                                                                  >  2:230504/1‑132 (MQ=255)
                                                                                   cATCAGCAGGCGCACCGCCGAAGGGCGCTGCAGGGCCACGAAATCTTCGGCCATCAGCTGTGCCCAGCTGACCTGGGGCAGCAGTGCCAAGGGCGAGTCGGCGGCCACTACGGCAACGAAGCGGTCCATGTa                                                                                  <  1:230504/132‑1 (MQ=255)
                                                                                   cATCAGCAGGCGCACCGCCGAAGGGCGCTGCAGGGCCACGAAATCTTCGGCCATCAGCTGTGCCCAGCTGACCTGGGGCAGCAGTGCCAAGGGCGAGTCGGCGGCCACTACGGCAACGAAGCGGTCCATGTACAACGGGTGGAAGGCCa                                                                 >  2:124336/1‑149 (MQ=255)
                                                                                        gcagGCGCACCGCCGAAGGGCGCTGCAGGGCCACGAAATCTTCGGCCATCAGCTGTGCCCAGCTGACCTGGGGCAGCAGTGCCAAGGGCGAGTCGGCGGCCACTACGGCAACGAAg                                                                                             <  2:307332/116‑1 (MQ=255)
                                                                                        gcagGCGCACCGCCGAAGGGCGCTGCAGGGCCACGAAATCTTCGGCCATCAGCTGTGCCCAGCTGACCTGGGGCAGCAGTGCCAAGGGCGAGTCGGCGGCCACTACGGCAACGAAg                                                                                             >  1:307332/1‑116 (MQ=255)
                                                                                        gcagGCGCACCGCCGAAGGGCGCTGCAGGGCCACGAAATCTTCGGCCATCAGCTGTGCCCAGCTGACCTGGGGCAGCAGTGCCAAGGGCGAGTCGGCGGCCACTACGGCAACGAAGCGGTCCATGTACAACGGGTGGAAGTCCAGCCCa                                                            >  5:335489/1‑149 (MQ=255)
                                                                                        gcagGCGCACCGCCGAAGGGCGCTGCAGGGCCACGAAATCTTCGGCCATCAGCTGTGCCCAGCTGACCTGGGGCAGCAGTGCCAAGGGCGAGTCGGCGGCCACTACGGCAACGAAGCGGTCCATGTACAACGGGTGGAAGTCCAGCCCa                                                            >  4:66430/1‑149 (MQ=255)
                                                                                            gcgcACCGCCGAAGGGCGCTGCAGGGCCACGAAATCTTCGGCCATCAGCTGTGCCCAGCTGACCTGGGGCAGCAGTGCCAAGGGCGAGTCGGCGGCCACTACGGCAACg                                                                                                <  2:31265/109‑1 (MQ=255)
                                                                                            gcgcACCGCCGAAGGGCGCTGCAGGGCCACGAAATCTTCGGCCATCAGCTGTGCCCAGCTGACCTGGGGCAGCAGTGCCAAGGGCGAGTCGGCGGCCACTACGGCAACg                                                                                                >  1:31265/1‑109 (MQ=255)
                                                                                            gcgcACCGCCGAAGGGCGCTGCAGGGCCACGAAATCTTCGGCCATCAGCTGTGCCCAGCTGACCTGGGGCAGCAGTGCCAAGGGCGAGTCGGCGGCCACTACGGCAACGAAGCGGTCCATGTACAACGGGTGGAAGTCCAGCCCATCGa                                                        <  3:66430/149‑1 (MQ=255)
                                                                                                 ccgccgAAGGGCGCTGCAGGGCCACGAAATCTTCGGCCATCAGCTGTGCCCAGCTGACCTGGGGCAGCAGTGCCAAGGGCGAGTCGGCGGCCACTACGGCAACGAAGCGGTCCATGTACAACGGGTGGAAGTCCAGCCCATCGATGTTa                                                   >  6:332208/1‑149 (MQ=255)
                                                                                                     cgAAGGGCGCTGCAGGGCCACGAAATCTTCGGCCATCAGCTGTGCCCAGCTGACCTGGGGCAGCAGTGCCAAGGGCGAGTCGGCGGCCACTACGGCAACGAAGCGGTCCATGTACAACGGGTGGAAGTCCAGCCCATCGATGTTATCCg                                               >  2:40283/1‑149 (MQ=255)
                                                                                                              cTGCAGGGCCACGAAATCTTCGGCCATCAGCTGTGCCCAGCTGACCTGGGGCAGCAGTGCCAAGGGCGAGTCGGCGGCCACTACGGCAACGAAGCGGTCCATGTACAACGGGTGGAAGTCCAGCCCATCGATGTTATCCGGCTCGAAGc                                      <  6:335489/149‑1 (MQ=255)
                                                                                                                     gCCACGAAATCTTCGGCCATCAGCTGTGCCCAGCTGACCTGGGGCAGCAGTGCCAAGGGCGAGTCGGCGGCCACTACGGCAACGAAGCGGTcc                                                                                       <  4:39829/93‑1 (MQ=255)
                                                                                                                     gCCACGAAATCTTCGGCCATCAGCTGTGCCCAGCTGACCTGGGGCAGCAGTGCCAAGGGCGAGTCGGCGGCCACTACGGCAACGAAGCGGTcc                                                                                       >  3:39829/1‑93 (MQ=255)
                                                                                                                       cACGAAATCTTCGGCCATCAGCTGTGCCCAGCTGACCTGGGGCAGCAGTGCCAAGGGCGAGTCGGCGGCCACTACGGCAACGAAGCGGTCCATGTACAACGGGTGGAAGTCCAGCCCATCGATGTTATCCGGCTCGAAGCCGATGcc                               >  2:292175/1‑147 (MQ=255)
                                                                                                                       cACGAAATCTTCGGCCATCAGCTGTGCCCAGCTGACCTGGGGCAGCAGTGCCAAGGGCGAGTCGGCGGCCACTACGGCAACGAAGCGGTCCATGTACAACGGGTGGAAGTCCAGCCCATCGATGTTATCCGGCTCGAAGCCGATGcc                               <  1:292175/147‑1 (MQ=255)
                                                                                                                                  cGGCCATCAGCTGTGCCCAGCTGACCTGGGGCAGCAGTGCCAAGGGCGAGTCGGCGGCCACTACGGCAACGAAGCGGTCCATGTACAACGGGTGGAAGTCCAGCCCATCGATGTTATCCGGCTCGAAGCCGATGCCCAGTTCGAc                      <  1:40283/145‑1 (MQ=255)
                                                                                                                                        tCAGCTGGGCACAGCTGACCTGGGGCAGCAGTGCCAAGGGCGAGTCGGCGGCCACTACGGCAACGAAGCGGTCCATGTACAACGGGTGGAAGTCCAGCCCATCGATGTTATCCGGCTCGAAGCCGATGCCCAGTTCGACGCGGCGATgg            >  1:256137/1‑149 (MQ=255)
                                                                                                                                          aGCTGTGCCCAGCTGACCTGGGGCAGCAGTGCCAAGGGCGAGTCGGCGGCCACTACGGCAACGAAGCGGTCCATGTACAACGGGTGGAAGTCCAGCCCATCGATGTTATCCGGCTCGAAGCCGATGCCCAGTTCGACGCGGCGATg             <  2:201505/146‑1 (MQ=255)
                                                                                                                                          aGCTGTGCCCAGCTGACCTGGGGCAGCAGTGCCAAGGGCGAGTCGGCGGCCACTACGGCAACGAAGCGGTCCATGTACAACGGGTGGAAGTCCAGCCCATCGATGTTATCCGGCTCGAAGCCGATGCCCAGTTCGACGCGGCGATg             >  1:201505/1‑146 (MQ=255)
                                                                                                                                           gCTGTGCCCAGCTGACCTGGGGCAGCAGTGCCAAGGGCGAGTCGGCGGCCACTACGGCAACGAAGCGGTCCATGTACAACGGGTGGAAGTCCAGCCCATCGATGTTATCCGGCTCGAAGCCGATGCCCAGTTCGACGCGGCGATGGCg          >  4:280099/1‑148 (MQ=255)
                                                                                                                                                 cccAGCTGACCTGGGGCAGCAGTGCCAAGGGCGAGTCGGCGGCCACTACGGCAACGAAGCGGTCCATGTACAACGGGTGGAAGTCCAGCCCATCGATGTTATCCGGCTCGAAGCCGATGCCCAGTTCGACGCGGCGATGGCGTACCAGt   >  2:158210/1‑149 (MQ=255)
                                                                                                                                                  ccAGCTGACCTGGGGCAGCAGTGCCAAGGGCGAGTCGGCGGCCACTACGGCAACGAAGCGGTCCATGTACAACGGGTGGAAGTCCAGCCCATCGATGTTATCCGGCTCGAAGCCGATGCCCAGTTCGACGCGGCGATGGCGTACCAGtt  >  6:330414/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GCCACCATGCGGCCGATGGTCGTCAGCTGGTGGCTTTCGAACGCGACTGCCAGCTTGCCGTGGCGGGCCGCCACATTCTGCTCCATCAGCAGGCGCACCGCCGAAGGGCGCTGCAGGGCCACGAAATCTTCGGCCATCAGCTGTGCCCAGCTGACCTGGGGCAGCAGTGCCAAGGGCGAGTCGGCGGCCACTACGGCAACGAAGCGGTCCATGTACAACGGGTGGAAGTCCAGCCCATCGATGTTATCCGGCTCGAAGCCGATGCCCAGTTCGACGCGGCGATGGCGTACCAGTT  >  NC_002947/3532488‑3532782

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