Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,260,347 G→C 7.6% T158S (ACC→AGC) ‡ PP_2857 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,260,3470GC 7.6% 191.0 / 3.6 63T158S (ACC→AGC) ‡PP_2857hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (29/28);  new base C (2/3);  total (31/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 = 1.00e+00

CGTACGCTATTAGGTTAAGTGTCGTTCTTAACTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCGGTTGTGCACCAATGCTGATGCTCAGTGCGTTGGCTTGCTGGGTGTTGG  >  NC_002947/3260202‑3260494
                                                                                                                                                  |                                                                                                                                                   
cGTACGCTATTAGGTTAAGTGTCGTTCTTAACTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTgg                                                                                                                                                     <  7:115830/146‑1 (MQ=255)
cGTACGCTATTAGGTTAAGTGTCGTTCTTAACTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAg                                                                                                                                                  >  5:234431/1‑149 (MQ=255)
     gCTATTAGGTTAAGTGTCGTTCTTAACTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGAcc                                                                                                                                             <  6:128018/149‑1 (MQ=255)
       tATTAGGTTAAGTGTCGTTCTTAACTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCAGGTCACTACCAGTGAccag                                                                                                                                           >  2:341061/1‑147 (MQ=255)
       tATTAGGTTAAGTGTCGTTCTTAACTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTAGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCATCGGCTACCGGTGAAGGCCAGGTCACTACCAGTGAc                                                                                                                                              >  7:139465/1‑146 (MQ=255)
         ttAGGTTTAGTGTCGTTGTTAAATGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCAGGTCACTACCAGTGAcca                                                                                                                                            <  8:139465/146‑2 (MQ=255)
         ttAGGTTAAGTGTCGTTCTTAACTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTg                                                                                                                                           <  8:172541/147‑1 (MQ=255)
         ttAGGTTAAGTGTCGTTCTTAACTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCAGGTCACTACCAGTGAccag                                                                                                                                           <  2:279227/147‑3 (MQ=255)
           aGGTTAAGTGTCGTTCTTAACTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCAGGTCACTACCAGTGACCAGGcc                                                                                                                                        <  1:341061/148‑2 (MQ=255)
                     tCGTTCTTAACTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCctgctg                                                                                                                             <  7:134996/149‑1 (MQ=255)
                        ttCTTAACTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGaa                                                                                                                          <  3:209580/149‑1 (MQ=255)
                        ttCTTAACTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGaa                                                                                                                          <  6:297177/149‑1 (MQ=255)
                        tgattAACTGCCGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGaa                                                                                                                          <  6:234431/146‑1 (MQ=255)
                                tGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGtt                                                                                                                  >  1:266817/1‑149 (MQ=255)
                                  cAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGtt                                                                                                                  >  2:179846/1‑147 (MQ=255)
                                  cAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGtt                                                                                                                  >  2:174358/1‑147 (MQ=255)
                                  cAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGtt                                                                                                                  <  1:179846/147‑1 (MQ=255)
                                  cAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGtt                                                                                                                  <  1:174358/147‑1 (MQ=255)
                                    gggTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTc                                                                                                              <  1:356921/149‑1 (MQ=255)
                                               cgACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGttgttg                                                                                                   <  8:327371/149‑1 (MQ=255)
                                                           tCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACGTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGcc                                                                                                                                  >  1:238924/1‑107 (MQ=255)
                                                           tCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGcc                                                                                                                                  <  2:238924/106‑1 (MQ=255)
                                                           tCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACg                                                                                       >  2:25337/1‑149 (MQ=255)
                                                           tCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACg                                                                                       >  8:275096/1‑149 (MQ=255)
                                                           tCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACg                                                                                       >  8:217284/1‑149 (MQ=255)
                                                             aaGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGtt                                                                                     >  1:166224/1‑149 (MQ=255)
                                                                   ttGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGtt                                                                                     >  5:130443/1‑143 (MQ=255)
                                                                   ttGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGAGCCACGtt                                                                                     <  6:130443/143‑1 (MQ=255)
                                                                     ggCCATTCGGTTTCCCACCCCCGCACTCTCCTTCAACTGGAACACCCCTCGGCTACCGGTGAAGGCCAGGTCACTGGTAGTGACCTTGCGAAAGCATGCTGGATTGGCAGTTGCTCCTGAGTTGGTGATCAGCGACACGTTCTGTTGCa                                                                             >  7:9399/1‑149 (MQ=255)
                                                                         aTTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAg                                                                         >  3:142127/1‑149 (MQ=255)
                                                                         aTTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAg                                                                         >  7:118427/1‑149 (MQ=255)
                                                                         aTTCGGTTTCCCACCCCGGCACTCTGGTTAAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAg                                                                         >  2:316861/1‑149 (MQ=255)
                                                                             ggTTTCCCACCCCGGCACTCTTGTTCAACTGGATCACCCCTCGGCTACCTTTGAAGGCCTGGTCACTGTTAGTGACCTGGAGATAGCCTTCTTGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGAc                                                                      <  3:299929/148‑1 (MQ=255)
                                                                               tttCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGTCTACCGGTGAAGGCCTTGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCt                                                                   <  6:182170/149‑1 (MQ=255)
                                                                                  cccACCCCGGCAGTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCAGTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTgtc                                                                >  6:173983/1‑149 (MQ=255)
                                                                                   ccACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGAcc                                                                                                                                             <  5:94911/71‑1 (MQ=255)
                                                                                   ccACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGAcc                                                                                                                                             >  6:94911/1‑71 (MQ=255)
                                                                                    cACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTgtcgt                                                              >  1:314396/1‑149 (MQ=255)
                                                                                    cACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGAGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTGGTTGATCAGCGCCACCTTCAGTTGCATGAGGACGCTgtcgt                                                              >  7:79680/1‑149 (MQ=255)
                                                                                      ccccGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTgtcgtc                                                             >  4:28976/1‑148 (MQ=255)
                                                                                         cGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATg                                                         <  1:316861/149‑1 (MQ=255)
                                                                                         cGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATg                                                         <  8:302879/149‑1 (MQ=255)
                                                                                                  tGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTgtcgt                                                              >  1:197642/1‑136 (MQ=255)
                                                                                                  tGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTgtcgt                                                              <  2:197642/136‑1 (MQ=255)
                                                                                                  tGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCgg                                                 >  1:93676/1‑149 (MQ=255)
                                                                                                   ggTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCGGt                                                >  2:80381/1‑149 (MQ=255)
                                                                                                    gTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCg                                                  <  4:142127/146‑1 (MQ=255)
                                                                                                      tCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGt                                                                                 <  4:304187/113‑1 (MQ=255)
                                                                                                      tCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGt                                                                                 >  3:304187/1‑113 (MQ=255)
                                                                                                       cAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCGGTtgt                                             <  2:266817/148‑1 (MQ=255)
                                                                                                       cAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCGGTtgt                                             <  7:275096/148‑1 (MQ=255)
                                                                                                             gATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCGGTTGTGCACCaa                                      <  2:166224/149‑1 (MQ=255)
                                                                                                                   cccTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCGGtt                                               >  4:349469/1‑134 (MQ=255)
                                                                                                                   cccTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCGGtt                                               <  3:349469/134‑1 (MQ=255)
                                                                                                                    ccTCGTCTACCGGTGAAGGCCTGGTCTCTGGTAGTGACCGGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTGGTTGATCAGCGGCAGGTTCTGTTGCATGAGGAGGCTGTCGTCGATGCTTTGCGGTTGTGCAGCAATGCTGAt                               <  8:9399/149‑1 (MQ=255)
                                                                                                                    ccTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGATGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCGGTTGTGCACCAATGCTGAt                               >  1:108054/1‑149 (MQ=255)
                                                                                                                      tCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCGGTTGTGCACCAATGCTg                                 <  2:24304/145‑1 (MQ=255)
                                                                                                                      tCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCGGTTGTGCACCAATGCTg                                 >  1:24304/1‑145 (MQ=255)
                                                                                                                           tACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAg                                                                                                         <  7:61710/68‑1 (MQ=255)
                                                                                                                           tACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAg                                                                                                         >  8:61710/1‑68 (MQ=255)
                                                                                                                                gTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCGGTTGTGCACCAATGCTGATGCTCAGTGCGtt                   <  7:217284/149‑1 (MQ=255)
                                                                                                                                        cTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTGAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCGGTTGTGCACCTATGCTGATGCTGAGTGCGTTggctggct           >  4:220224/1‑149 (MQ=255)
                                                                                                                                        cTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCGGTtgt                                             >  6:9675/1‑115 (MQ=255)
                                                                                                                                        cTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCGGTtgt                                             <  5:9675/115‑1 (MQ=255)
                                                                                                                                                  gTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCGGTTGTGCACCAATGCTGATGCTCAGTGCGTTGGCTTGCTGGGTGTTg   <  8:118427/147‑1 (MQ=255)
                                                                                                                                                  gTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCGGTTGTGCACCAATGCTGATGCTCAGTGCGTTGGCTTGCTGGGTGTTg   <  2:93676/147‑1 (MQ=255)
                                                                                                                                                  gTAGTGACCTCGCGATAGCCTGCTGGAATGTCAGTTGCACCTGAGTTGTTGATCAGCGCCACGTTCAGTTGCATGAGGACGCTGTCGACGATGCTTTGCCGTTGTGCACCAATGCTGATGCGCACTCCCTTGGCTTGCTGGGAGTTgg  >  2:274306/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CGTACGCTATTAGGTTAAGTGTCGTTCTTAACTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCGGTTGTGCACCAATGCTGATGCTCAGTGCGTTGGCTTGCTGGGTGTTGG  >  NC_002947/3260202‑3260494

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: