Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,137,445 A→G 5.6% C53R (TGC→CGC)  PP_4522 ← LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,137,4450AG 5.6% 205.8 / 2.1 71C53R (TGC→CGC) PP_4522LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (33/34);  new base G (2/2);  total (35/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 = 8.45e-01

CAAGGTTGAGGGTGCCACGCAGTTCGCCCTTCAGGGCGGCGGCATACTGCTCGAAACCGTCCAGCTCGGCCAGCAGGCGCAGGGTTTCCTGATGGAACAGCTCGCCCTTGCTGGTCAGGCTGAAGCCCCCCCGACCCCGGTGGCACAACACGATACCCAAGGCGCCCTCCAGCTGGCTCATGTAGGTACTGATGGCCGAGGTCGACAGGTTCAGCTCGCGCTGGGCGTTGGCAAAGCCCTGGTGGCGCACCACGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTC  >  NC_002947/5137301‑5137593
                                                                                                                                                |                                                                                                                                                     
cAAGGTTGAGGGTGCCACGCAGTTCGCCCTTCAGGGGGGCGGCATACTGCTCGAAACGGTCCAGCTCTTCCAGCAGGCGCAGGGTTTCCTGATGGAACAGCTCGCCCTTGCTGGTCAGGCTGAAGCCCCCCCGACCCCGGTGGCacaac                                                                                                                                                   <  6:221450/149‑1 (MQ=255)
   ggTTGAGGGTGCCACGCAGTTCGCCCTTCAGGGCGGCGGCATACTGCTCGAAACCGTCCAGCTCGGCCAGCAGGCGCAGGGTTTCCTGATGGAACAGCTCGCCCTTGCTGGTCAGGCTGAAGCCCCCCCGACCCCGGTGGCACAacac                                                                                                                                                 >  1:87481/1‑148 (MQ=255)
   ggTTGAGGGTGCCACGCAGTTCGCCCTTCAGGGCGGCGGCATACTGCTCGAAACCGTCCAGCTCGGCCAGCAGGCGCAGGGTTTCCTGATGGAACAGCTCGCCCTTGCTGGTCAGGCTGAAGCCCCCCCGACCCCGGTGGCACAACACg                                                                                                                                                <  4:187548/149‑1 (MQ=255)
               cACGCAGTTCGCCCTTCAGGGCGGCGGCATACTGCTCGAAACCGTCCAGCTCGGCCAGCAGGCGCAGGGTTTCCTGATGGAACAGCTCGCCCTTGCTGGTCAGGCTGAAGCCCCCCCGACCCCGGTGGCACAACACGATa                                                                                                                                             >  8:101656/1‑140 (MQ=255)
               cACGCAGTTCGCCCTTCAGGGCGGCGGCATACTGCTCGAAACCGTCCAGCTCGGCCAGCAGGCGCAGGGTTTCCTGATGGAACAGCTCGCCCTTGCTGGTCAGGCTGAAGCCCCCCCGACCCCGGTGGCACAACACGATa                                                                                                                                             <  7:101656/140‑1 (MQ=255)
                     cgccGCCCTGAAGGGCGGCGGCATACTGCTCGAAACCGTCCAGCTCGGCCAGCAGGCGCAGGGTTTCCTGATGGAACAGCTCGCCCTTGCTGGTCAGGCTGAAGCCCCCCCGACCCCGGTGGCACAACACGATACCCAAGGCGCCCTc                                                                                                                               >  1:26594/4‑148 (MQ=255)
                       ccGCCCTGAAGGGCGGCGGCATACTGCTCGAAACCGTCCAGCTCGGCCAGCAGGCGCAGGGTTTCCTGATGGAACAGCTCGCCCTTGCTGGTCAGGCTGAAGCCCCCCCGACCCCGGTGGCACAACACGATACCCAAGGCGCCCTCCAg                                                                                                                            >  2:1189/2‑149 (MQ=255)
                         gCCCTTCAGGGCGGCGGCATACTGCTCGAAACCGTCCAGCTCGGCCAGCAGGCGCAGGGTTTCCTGATGGAACAGCTCGCCCTTGCTGGTCAGGCTGAAGCCCCCCCGACCCCGGTGGCACAACACGATACCCAAGGCGCCCTCCAGCt                                                                                                                          >  8:170041/1‑149 (MQ=255)
                           ccTTCAGGGCGGCGGCATACTGCTCGAAACCGTCCAGCTCGGCCAGCAGGCGCAGGGTTTCCTGATGGAACAGCTCGCCCTTGCTGGTCAGGCTGAAGCCCCCCCGACCCCGGTGGCACAACACGATACCCAAGGCGCCCTCCAGCTgg                                                                                                                        <  2:229382/149‑1 (MQ=255)
                           ccTTCAGGGCGGCGGCATACTGCTCGAAACCGTCCAGCTCGGCCAGCAGGCGCAGGGTTTCCTGATGGAACAGCTCGCCCTTGCTGGTCAGGCTGAAGCCCCCCCGACCCCGGTGGCACAACACGATACCCAAGGCGCCCTCCAGCTgg                                                                                                                        >  3:71622/1‑149 (MQ=255)
                           ccTTCAGGGCGGCGGCATACTGCTCGAAACCGTCCAGCTCGGCCAGCAGGCGCAGGGTTTCCTGATGGAACAGCTCGCCCTTGCTGGTCAGGCTGAAGCCCCCCCGACCCCGGTGGCACAACACGATACCCAAGGCGCCCTCCAGCTgg                                                                                                                        >  3:93031/1‑149 (MQ=255)
                           ccTTCAGGGCGGCGGCATACTGCTCGAAACCGTCCAGCTCGGCCAGCAGGCGCAGGGTTTCCTGATGGAACAGCTCGCCCTTGCTGGTCAGGCTGAAGCCCCCCCGACCCCGGTGGCACAACACGATACCCAAGGCGCCCTCCAGCTgg                                                                                                                        <  4:210404/149‑1 (MQ=255)
                           ccTTCAGGGCGGCGGCATACTGCTCGAAACCGTCCAGCTCGGCCAGCAGGCGCAGGGTTTCCTGATGGAACAGCTCGCCCTTGCTGGTCAGGCTGAAGCCCCCCCGACCCCGGTGGCACAACACGATACCCAAGGCGCCCTCCAGCTgg                                                                                                                        <  5:163562/149‑1 (MQ=255)
                           ccTTCAGGGCGGCGGCATACTGCTCGAAACCGTCCAGCTCGGCCAGCAGGCGCAGGGTTTCCTGATGGAACAGCTCGCCCTTGCTGGTCAGGCTGAAGCCCCCCCGACCCCGGTGGCACAACACGATACCCAAGGCGCCCTCCAGCTgg                                                                                                                        <  2:87481/149‑1 (MQ=255)
                           ccTTCAGGGCGGCGGCATACTGCTCGAAACCGTCCAGCTCGGCCAGCAGGCGCAGGGTTTCCTGATGGAACAGCTCGCCCTTGCTGGTCAGGCTGAAGCCCCCCCGACCCCGGTGGCACAACACGATACCCAAGGCGCCCTCCAGCTgg                                                                                                                        >  6:73390/1‑149 (MQ=255)
                           ccTTCAGGGCGGCGGCATACTGCTCGAAACCGTCCAGCTCGGCCAGCAGGCGCAGGGTTTCCTGATGGAACAGCTCGCCCTTGCTGGTCAGGCTGAAGCCCCCCCGACCCCGGTGGCACAACACGATACCCAAGGCGCCCTCCAGCTgg                                                                                                                        <  8:25020/149‑1 (MQ=255)
                           ccTTCAGGGCGGCGGCATACTGCTCGAAACCGTCCAGCTCGGCCAGCAGGCGCAGGGTTTCCTGATGGAACAGCTCGCCCTTGCTGGTCAGGCTGAAGCCCCCCCGACCCCGGTGGCACAACACGATACCCAAGGCGCCCTCCAGCTgg                                                                                                                        >  1:96497/1‑149 (MQ=255)
                            cTTCAGGGCGGCGGCATACTGCTCGAAACCGTCCAGCTCGGCCAGCAGGCGCAGGGTTTCCTGATGGAACAGCTCGCCCTTGCTGGTCAGGCTGAAGCCCCCCCGACCCCGGTGGCACAACACGATACCCAAGGCGCCCTCCAGCTGGc                                                                                                                       <  7:281133/149‑1 (MQ=255)
                                         gcATACTGCTCGAAACCGTCCAGCTCGGCCAGCAGGCGCAGGGTTTCCTGATGGAACAGCTCGCCCTTGCTGGTCAGGCTGAAGCCCCCCCGACCCCGGTGGCACAACACGATACCCAAGGCGCCCTCCAGCTGGCTCATGTAGGTACt                                                                                                         >  1:405075/1‑149 (MQ=255)
                                            tACTGCTCGAAACCGTCCAGCTCGGCCAGCAGGCGCAGGGTTTCCTGATGGAACAGCTCGCCCTTGCTGGTCAGGCTGAAGCCCCCCCGACCCCGGTGGCACAACACGATACCCAAGGCGCCCTCCAGCTGGCTCATGTAGGTACTGAt                                                                                                      <  7:41493/149‑1 (MQ=255)
                                                gCTCGAAACCGTCCAGCTCGGCCAGCAGGCGCAGGGTTTCCTGATGGAACAGCTCGCCCTTGCTGGTCAGGCTGAAGCCCCCCCGCCCCCGGTGGCACAACACGATACCCAAGGCGCCCTCCAGCTGGCTCATGTAGGTACTGATGGcc                                                                                                  >  2:112062/1‑149 (MQ=255)
                                                         cctCCAGCTCGGCCAGCAGGCGCAGGGTTTCCTGATGGAACAGCTCGCCCTTGCTGGTCAGGCTGAAGCCCCCCCGACCCCGGTGGCACAACACGATACCCAAGGCGCCCTCCAGCTGGCTCATGTAGGTACTGATGGCCGAGGTCGAc                                                                                         <  8:366880/147‑1 (MQ=255)
                                                          gTCCAGCTCGGCCAGCAGGCGCAGGGTTTCCTGATGGAACAGCTCGCCCTTGCTGGTCAGGCTGAAGCCCCCCCGACCCCGGTGGCACAACACGATACCCAAGGCGCCCTCCAGCTGGCTCATGTAGGTACTGATGGCCGAGGTCGa                                                                                          >  2:288206/1‑147 (MQ=255)
                                                          gTCCAGCTCGGCCAGCAGGCGCAGGGTTTCCTGATGGAACAGCTCGCCCTTGCTGGTCAGGCTGAAGCCCCCCCGACCCCGGTGGCACAACACGATACCCAAGGCGCCCTCCAGCTGGCTCATGTAGGTACTGATGGCCGAGGTCGACa                                                                                        <  4:93031/149‑1 (MQ=255)
                                                          gTCCAGCTCGGCCAGCAGGCGCAGGGTTTCCTGATGGAACAGCTCGCCCTTGCTGGTCAGGCTGAAGCCCCCCCGACCCCGGTGGCACAACACGATACCCAAGGCGCCCTCCAGCTGGCTCATGTAGGTACTGATGGCCGAGGTCGACa                                                                                        <  7:156532/149‑1 (MQ=255)
                                                            ccAGCTCGGCCAGCAGGCGCAGGGTTTCCTGATGGAACAGCTCGCCCTTGCTGGTCAGGCTGAAGCCCCCCCGACCCCGGTGGCACAACACGATACCCAAGGCGCCCt                                                                                                                                <  8:163435/108‑1 (MQ=255)
                                                            ccAGCTCGGCCAGCAGGCGCAGGGTTTCCTGATGGAACAGCTCGCCCTTGCTGGTCAGGCTGAAGCCCCCCCGACCCCGGTGGCACAACACGATACCCAAGGCGCCCt                                                                                                                                >  7:163435/1‑108 (MQ=255)
                                                              aGCTCGGCCAGCAGGCGCAGGGTTTCCTGATGGAACAGCTCGCCCTTGCTGGTCAGGCTGAAGCCCCCCCGACCCCGGTGGCACAACACGATACCCAAGGCGCCCTCCAGCTGGCTCATGTAGGTACTGATGGCCGAGGTCGACAGGtt                                                                                    <  5:73390/149‑1 (MQ=255)
                                                                  cGGCCAGCAGGCGCAGGGTTTCCTGATGGAACAGCTCGCCCTTGCTGGTCAGGCTGAAGCCCCCCCGACCCCGGTGGCACAACACGATACCCAAGGCGCCCTCCAGCTGGCTCATGTAGGTACTGATGGCCGAGGTCGACAGGtt                                                                                    >  2:407142/1‑145 (MQ=255)
                                                                  cGGCCAGCAGGCGCAGGGTTTCCTGATGGAACAGCTCGCCCTTGCTGGTCAGGCTGAAGCCCCCCCGACCCCGGTGGCACAACACGATACCCAAGGCGCCCTCCAGCTGGCTCATGTAGGTACTGATGGCCGAGGTCGACAGGtt                                                                                    <  1:407142/145‑1 (MQ=255)
                                                                   ggCCAGCAGGCGCAGGGTTTCCTGATGGAACAGCTCGCCCTTGCTGGTCAGGCTGAAGCCCCCCCGACCCCGGTGGCACAACACGATACCCAAGGCGCCCTCCAGCTGGCTCATGTAGGTACTGATGGCCGAGGTCGACAGGTTCAg                                                                                 <  2:405075/147‑1 (MQ=255)
                                                                   ggCCAGCAGGCGCAGGGTTTCCTGATGGAACAGCTCGCCCTTGCTGGTCAGGCTGAAGCCCCCCCGACCCCGGTGGCACAACACGATACCCAAGGCGCCCTCCAGCTGGCTCATGTAGGTACTGATGGCCGAGGTCGACAGGTTCAGc                                                                                <  5:22173/148‑1 (MQ=255)
                                                                   ggCCAGCAGGCGCAGGGTTTCCTGATGGAACAGCTCGCCCTTGCTGGTCAGGCTGAAGCCCCCCCGACCCCGGTGGCACAACACGATACCCAAGGCGCCCTCCAGCTGGCTCATGTAGGTACTGATGGCCGAGGTCGACAGGTTCAGc                                                                                <  2:273233/148‑1 (MQ=255)
                                                                      cagcagGCGCAGGGTTTCCTGATGGAACAGCTCGCCCTTGCTGGTCAGGCTGAAGCCCCCCCGACCCCGGTGGCACAACACGATACCCAAGGCGCCCTCCAGCTGGCTCATGTAGGTACTGATGGCCGAGGTCGACAGGTTCAGCTcgc                                                                            <  1:1189/149‑1 (MQ=255)
                                                                       agcagGCGCAGGGTTTCCTGATGGAACAGCTCGCCCTTGCTGGTCAGGCTGAAGCCCCCCCGACCCCGGTGGCACAACACGATACCCAAGGCGCCCTCCAGCTGGCTCATGTAGGTACTGATGGCCGAGGTCGACAGGTTCAGCTcgcg                                                                           >  4:342403/1‑149 (MQ=255)
                                                                                 tgggttCCTGATGGAACATCTCGCCCTTGCCGGTGATGCCTAAGGCCCCCCGACCCCGGTGGCAAAACCCGATACCCAAGGCCCCCTCCATCTGGCTCATGTTCCGTACTGATGGCCGAGGTCGACAGGTTCAGCTCGCGCTGGGCGtt                                                                  <  1:288206/145‑1 (MQ=255)
                                                                                  gctttCCTGATGGAACACCTCGCCCTTGCTGGTCAGGCTGAAGCCCCCCCGACCCCGGTGGCACACCACGAGACCCAAGGCCACCTCCAGCTGGAGCATGTAGGTACTGAGGGCCGAGCTCGCCAGGTACAGCTCGCGCGGGGCGTTgc                                                                >  4:254571/3‑148 (MQ=255)
                                                                                    tttCCTGATGGAACAGCTCGCCCTTGCTGGTCAGGCTGAAGCCCCCCCGACCCCGGTGGCGCAACACGATACCCAAGGCGCCCTCCAGCTGGCTCATGTAGGTACTGATGGCCGAGGTCGACAGGTTCAGCTCGCGCTGGGCGTTGGc                                                               >  1:180938/1‑148 (MQ=255)
                                                                                    tttCCTGATGGAACAGCTCGCCCTTGCTGGTCAGGCTGAAGCCCCCCCGACCCCGGTGGCACAACACGATACCCAAGGCGCCCTCCAGCTGGCTCATGTAGGTACTGATGGCCGAGGTCGACAGGTTCAGCTCGCGCTGGGCGTTGGc                                                               >  4:119060/1‑148 (MQ=255)
                                                                                    tttCCTGATGGAACAGCTCGCCCTTGCTGGTCAGGCTGAAGCCCCCCCGACCCCGGTGGCACAACACGATACCCAAGGCGCCCTCCAGCTGGCTCATGTAGGTACTGATGGCCGAGGTCGACAGGTTCAGCTCGCGCTGGGCGTTGGc                                                               >  5:187855/1‑148 (MQ=255)
                                                                                      tCCTGATGGAACAGCTCGCCCTTGCTGGTCAGGCTGAAGCCCCCCCGACCCAGGTGGCACAACACGATACCCAAGGCGCCCTGCAGCTGGCTCATGTAGGTACTGAGGGCCGAGGTCGACAGGTTCAGCTCGCGCTGGGCGTTGGCa                                                              >  6:78148/1‑147 (MQ=255)
                                                                                       ccTGATGGAACAGCTCGCCCTTGCTGGTCAGGCTGAAGCCCCCCCGACCCCGGTGGCACAACACGATACCCAAGGCGCCCTCCAGCTGGCTCATGTAGGTACTGATGGCCGAGGTCGACAGGTTCAGCTCGCGCTGGGCGTTGGCaaa                                                            <  2:96497/148‑1 (MQ=255)
                                                                                        cTGATGGAACAGCTCGCCCTTGCTGGTCAGGCTGAAGCCCCCCCGACCCCGGTGGCACAACACGATACCCAAGGCGCCCTCCAGCTGGCTCATGTAGGTACTGATGGCCGAGGTCGACAGGTTCAGCTCGCGCTGGGCGTTGGCAAAGc                                                          >  4:116610/1‑149 (MQ=255)
                                                                                           aTGGAACAGCTCGCCCTTTCTGGTCAGGCTGAAGCCCCCCCGACCCCGGTGGCACAACACGATACCCAAGGCGCCCTCCAGCTGGCTCATGTAGGTACTGATGGCCGAGGTCGACAGGTTCAGCTCGCGCTGGGCGTTGGCAAAGcc                                                         <  6:187855/147‑1 (MQ=255)
                                                                                            tGGAACAGCTCGCCCTTGCTGGTCAGGCTGAAGCCCCCCCGACCCCGGTGGCACAACACGATACCCAAGGCGCCCTCCAGCTGGCTCATGTAGGTACTGATGGCCGAGGTCGACAGGTTCAGCTCGCGCTGGGCTTTGCCAACGccca                                                       <  4:71622/148‑2 (MQ=255)
                                                                                              gAACAGCTCGCCCTTGCTGGTCAGGCTGAAGCCCCCCCGACCCCGGTGGCACAACACGATACCCAAGGCGCCCTCCAGCTGGCTCATGTAGGTACTGATGGACGAGGTCGACAGGTTCAGCTCGCGCTGGGCGTTGGCAAAGCCCtggt                                                    <  1:293152/149‑1 (MQ=255)
                                                                                               aaCAGCTCGCCCTTGCTGGTCAGGCTGAAGCCCCCCCGACCCCGGTGGCACAACACGATACCCAAGGCGCCCTCCAGCTGGCTCATGTAGGTACTGATGGCCGAGGTCGACAGGTTCAGCTCGCGCTGGGCGTTGGCAAAGCCCtggtg                                                   >  6:39837/1‑149 (MQ=255)
                                                                                                  aGCTCGCCCTTGCTGGTCAGGCTGAAGCCCCCCCGACCCCGGTGGCACAACACGATACCCAAGGCGCCCTCCAGCTGGCTCATGTAGGTACTGATGGCCGAGGTCGAc                                                                                         >  4:361939/1‑108 (MQ=255)
                                                                                                  aGCTCGCCCTTGCTGGTCAGGCTGAAGCCCCCCCGACCCCGGTGGCACAACACGATACCCAAGGCGCCCTCCAGCTGGCTCATGTAGGTACTGATGGCCGAGGTCGAc                                                                                         <  3:361939/108‑1 (MQ=255)
                                                                                                       gCCCTTGCTGGTCAGGCTGAAGCCCCCCCGACCCCGGTGGCACAACACGATACCCAAGGCGCCCTCCAGCTGGCTCATGTAGGTACTGATGGCCGAGGTCGACAGGTTCAGCTCGCGCTGGGCGTTGGCAAAGCCCTGGTGGCGcacca                                           >  1:42483/1‑149 (MQ=255)
                                                                                                              cTGGTCAGGCTGAAGCCCGCCCGACCCCGGTGGCACAACACGATACCCAAGGCGCCCTCCAGCTGGCTCATGTAGGTACTGATGGCCGAGGTCGACAGGTTCAGCTCGCGCTGGGCGTTGGCAAAGCCCTGGTGGCGCACCACGCtga                                     <  2:26594/148‑1 (MQ=255)
                                                                                                                ggTCAGGCTGAAGCCCCCCCGACCCCGGTGGCACAACACGATACCCAAGGCGCCCTCCAGCTGGCTCATGTAGGTACTGATGGCCGAGGTCGACAGGTTCAGCTCGCGCTGGGCGTTg                                                                 >  3:151555/1‑118 (MQ=255)
                                                                                                                  tCAGGCTGAAGCCCCCCCGACCCCGGTGGCACAACACGATACCCTAGGCGCCCTCCAGCTGGCTCATGt                                                                                                                >  4:295291/1‑69 (MQ=255)
                                                                                                                  tCAGGCTGAAGCCCCCCCGACCCCGGTGGCACAACACGATACCCTAGGCGCCCTCCAGCTGGCTCATGt                                                                                                                <  3:295291/69‑1 (MQ=255)
                                                                                                                  tCAGGCTGAAGCCCCCCCGACCCCGGTGGCACAACACGATACCCAAGGCGCCCTCCAGCTGGCTCATGTAGGTACTGATGGCCGAGGTCGACAGGTTCAGCTCGCGCTGGGCGTTGGCAAAGCCCTGGTGGCGCACCACGCtgatg                                   >  4:91752/1‑146 (MQ=255)
                                                                                                                   cAGGCTGAAGCCCCCCCGACCCCGGTGGCACAACACGATACCCAAGGCGCCCTCCAGCTGGCTCATGTAGGTACTGATGGCCGAGGTCGACAGGTTCAGCTCGCGCTGGGCGTTGGCAAAGCCCTGGTGGCGCACCACGCtgatg                                   >  3:142360/1‑145 (MQ=255)
                                                                                                                    aGGCTGAAGCCCCCCCGACCCCGGTGGCACAACACGATACCCAAGGCGCCCTCCAGCTGGCTCATGTAGGTACTGATGGCCGAGGTCGACAGGTTCAGCTCGCGCTGGGCGTTGGCAAAGCCCTGGTGGCGCACCACGCTGATGAAGa                               <  2:42483/148‑1 (MQ=255)
                                                                                                                    aGGCTGAAGCCCCCCCGACCCCGGTGGCACAACACGATACCCAAGGCGCCCTCCAGCTGGCTCATGTAGGTACTGATGGCCGAGGTCGACAGGTTCAGCTCGCGCTGGGCGTTGGCAAAGCCCTGGTGGCGCACCACGCTGATGAAGAt                              <  3:91752/149‑1 (MQ=255)
                                                                                                                      gCTGAAGCCCCCCCGACCCCGGTGGCACAACACGATACCCAAGGCGCCCTCCAGCTGGCTCATGTAGGTACTGATGGCCGAGGTCGACAGGTTCAGCTCgcgc                                                                          >  6:86886/1‑103 (MQ=255)
                                                                                                                      gCTGAAGCCCCCCCGACCCCGGTGGCACAACACGATACCCAAGGCGCCCTCCAGCTGGCTCATGTAGGTACTGATGGCCGAGGTCGACAGGTTCAGCTCgcgc                                                                          <  5:86886/103‑1 (MQ=255)
                                                                                                                        tGAAGCCCCCCCAACCCCGGTGGCGCAACACGATACCCAAGGCGCCCTCCAGCTGGCTCATGTAGGTACTGATGGCCGAGGTCGACAGGTTCAGCTCGCGCTGGGCGTTGGCAAAGCCCTGGTGGCGCACCACGCtgat                                    <  7:357467/139‑1 (MQ=255)
                                                                                                                        tGAAGCCCCCCCAACCCCGGTGGCGCAACACGATACCCAAGGCGCCCTCCAGCTGGCTCATGTAGGTACTGATGGCCGAGGTCGACAGGTTCAGCTCGCGCTGGGCGTTGGCAAAGCCCTGGTGGCGCACCACGCtgat                                    >  8:357467/1‑139 (MQ=255)
                                                                                                                            gCCCCCCCGACCCCGGTGGCACAACACGATACCCAAGGCGCCCTCCAGCTGGCTCATGTAGGTACTGATGGCCGAGGTCGACAGGTTCAGCTCGCGCTGGGCGTTGGCAAAGCCCTGGTGGCGCACCACGCTGATGAAGATGCgcag                        >  2:119229/1‑147 (MQ=255)
                                                                                                                             cccccccGACCCCGGTGGCACAACACGATACCCAAGGCGCCCTCCAGCTGGCTCTTGTAGGTACTGATGGCCGAGGTCGACAGGTTCAGCTCGCGCTGGGCGTTGGCAAAGCCCTGGTGGCGCACCAGGCTGATGAAGCTGCGcagcag                     >  3:210365/1‑149 (MQ=255)
                                                                                                                                ccccGACCCCGGTGGCACAACACGATACCCAAGGCGCCCTCCAGCTGGCTCATGTAGGTACTGATGGCCGAGGTCGACAGGTTCAGCTCGCGCTGGGCGTTGGCAAAGCCCTGGTGGCGCACCACGCTGATGAAGATGCGCAGCAGtt                   >  3:241097/1‑148 (MQ=255)
                                                                                                                                 cccGACCCCGGTGGCACAACACGATACCCAAGGCGCCCTCCAGCTGGCTCATGTAGGTACTGATGGCCGAGGTCGACAGGTTCAGCTCGCGCTGGGCGTTGGCAAAGCCCTGGTGGCGCACCACGCTGATGAAGATGCGCAGCAGttt                  >  6:27980/1‑148 (MQ=255)
                                                                                                                                 cccGACCCCGGTGGCACAACACGATACCCAAGGCGCCCTCCAGCTGGCTCATGTAGGTACTGATGGCCGAGGTCGACAGGTTCAGCTCGCGCTGGGCGTTGGCAAAGCCCTGGTGGCGCACCACGCTGATGAAGATGCGCAGCAGttt                  >  5:282419/1‑148 (MQ=255)
                                                                                                                                      ccccGGTGGCGCAACACGATACCCAAGGCGCCCTCCAGCTGGCTCATGTAGGTACTGATGGCCGAGGTCGACAGGTTCAGCTCGCGCTGGGCGTTGGCAAAGCCCTGGTGGCGCACCACGCTGATGAAGATGCGCAGCAGTTTCAGGTc            <  2:180938/149‑1 (MQ=255)
                                                                                                                                        gcGGTGGCACACCACGATACCCCAGGCGCCCTCCAGCTGGCTCATGTAGGTACTGATGGCCGAGGGCGACAGGTTCAGCTCGCGCTGGGCGTTGGCACAGCCCTGGTGGCGCACGACGCTGATGAAGATGCGCAGCAGTTTCAGGTCgg          <  3:254571/148‑1 (MQ=255)
                                                                                                                                         cGGTGGCACAACACGATACCCAAGGCGCCCTCCAGCTGGCTCATGTAGGTACTGATGGCCGAGGTCGACAGGTTCAGCTCGCGCTGGGCGTTGGCAAAGCCCTGGTGGCGCACCACGCTGATGAAGATGCGCAGCAGTTTCAGGTCggg         <  1:112062/149‑1 (MQ=255)
                                                                                                                                          ggTGGCACAACACGATACCCAAGGCGCCCTCCAGCTGGCTCATGTAGGTACTGATGGCCGAGGTCGACAGGTTCAGCTCGCGCTGGGCGTTGGCAAAGCCCTGGTGGCGCACCACGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGc        <  6:258523/149‑1 (MQ=255)
                                                                                                                                             ggCACAACACGATACCCAAGGCGCCCTCCAGCTGGCTCATGTAGGTACTGATGGCCGAGGTCGACAGGTTCAGCTCGCGCTGGGCGTTGGCAAAGCCCTGGTGGCGCACCACGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAg     <  4:142360/149‑1 (MQ=255)
                                                                                                                                                acaacaCGATACCCAAGGCGCCCTCCAGCTGGCTCATGTAGGTACTGATGGCCGAGGTCGACAGGTTCAGCTCGCGCTGGGCGTTGGCAAAGCCCTGGTGGCGCACCACGCTGATGAAGATGCgcagca                      >  2:362301/1‑129 (MQ=255)
                                                                                                                                                acaacaCGATACCCAAGGCGCCCTCCAGCTGGCTCATGTAGGTACTGATGGCCGAGGTCGACAGGTTCAGCTCGCGCTGGGCGTTGGCAAAGCCCTGGTGGCGCACCACGCTGATGAAGATGCgcagca                      <  1:362301/129‑1 (MQ=255)
                                                                                                                                                acaacaCGATACCCAAGGCGCCCTCCAGCTGGCTCATGTAGGTACTGATGGCCGAGGTCGACAGGTTCAGCTCGCGCTGGGCGTTGGCAAAGCCCTGGTGGCGCACCACGCTGATGAAGATGCGCAGCAGTTTGAGGTCGGGCAAGGTc  >  1:401689/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
CAAGGTTGAGGGTGCCACGCAGTTCGCCCTTCAGGGCGGCGGCATACTGCTCGAAACCGTCCAGCTCGGCCAGCAGGCGCAGGGTTTCCTGATGGAACAGCTCGCCCTTGCTGGTCAGGCTGAAGCCCCCCCGACCCCGGTGGCACAACACGATACCCAAGGCGCCCTCCAGCTGGCTCATGTAGGTACTGATGGCCGAGGTCGACAGGTTCAGCTCGCGCTGGGCGTTGGCAAAGCCCTGGTGGCGCACCACGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTC  >  NC_002947/5137301‑5137593

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: