Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,989,565 A→T 6.4% E272V (GAG→GTG)  PP_1780 → mannosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,989,5650AT 6.4% 223.4 / 3.3 73E272V (GAG→GTG) PP_1780mannosyltransferase
Reads supporting (aligned to +/- strand):  ref base A (33/35);  new base T (2/3);  total (35/38)
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

GCCTTCGCCAGCGCCTTCGATGCTGAGTATGGGGTGACCCTTGAGGACCTGGCTGAGCGCTATGAAGAAAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTCGAGCGTCATCGCGCCGAACTGGAGAGCCTCAAGCAAGTGGTGAACGGCCTGAACCACCATCTGTCCGATACCCA‑GCAACGGCTGGCTCTTGCTGATGAGC  >  NC_002947/1989420‑1989711
                                                                                                                                                 |                                                                                                                                                    
gCCTTCGCCAGCGCCTTCGATGCTGAGTATGGGGTGACCCTTGAGGACCTGGCTGAGCGCTATGAAGAAAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCt                                                                                                                                                   >  6:70763/1‑149 (MQ=255)
 ccTTCGCCAGCGCCTTCGATGCTGAGTATGGGGTGACCCTTGAGGACCTGGCTGAGCGCTATGAAGAAAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTg                                                                                                                                                  >  4:233850/1‑149 (MQ=255)
  cTTCGCCAGCTCCTTCGATGCTGAGTATGGGGTGACCCTTGAGGACCTGGCTGAGCGCTATGAAGAAAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTgg                                                                                                                                                 <  5:212087/149‑1 (MQ=255)
  cTTCGCCAGCGCCTTCGATGCTTAGTATGGGGTGACCCTTGAGGACCTGGCTGAGCGCTATGAAGAAAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTTGAGCTgg                                                                                                                                                 <  3:309417/149‑1 (MQ=255)
  cTTCGCCAGCGCCTTCGATGCTGAGTATGGGGTGACCCTTGAGGACCTGGCTGAGCGCTATGAAGAAAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTgg                                                                                                                                                 <  4:307384/149‑1 (MQ=255)
  cTTCGCCAGCGCCTTCGATGCTGAGTATGGGGTGACCCTTGAGGACCTGGCTGAGCGCTATGAAGAAAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTgg                                                                                                                                                 >  6:77363/1‑149 (MQ=255)
   ttCGCCAGCGCCTTCGATGCTGAGTATGGGGTGACCCTTGAGGACCTGGCTGAGCGCTATGAAGAAAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGTGCTgac                                                                                                                                                >  6:133389/1‑147 (MQ=255)
   ttCGCCAGCGCCTTCGATGCTGAGTATGGGGTGACCCTTGAGGACCTGGCTGAGCGCTATGAAGAAAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTggg                                                                                                                                                >  6:70996/1‑149 (MQ=255)
      gCCAGCGCCTTCGATGCTGAGTATGGGGTGACCCTTGAGGACCTGGCTGAGCGCTATGAAGCAAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGt                                                                                                                                               >  6:38649/1‑147 (MQ=255)
       ccAGCGCCTTCGATGCTGAGTATGGGGTGACCCTTGAGGACCTGGCTGAGCGCTATGAAGAAAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAg                                                                                                                                            >  4:320909/1‑149 (MQ=255)
       ccAGCGCCTTCGATGCTGAGTATGGGGTGACCCTTGAGGACCTGGCTGAGCGCTATGAAGAAAGGCTGGCAAGCAACTTCAAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATACGGCATTGTAGCAAGTTAACCTGGAGCTGGCTCAg                                                                                                                                            >  6:11253/1‑149 (MQ=255)
             ccTTCGATGCTGAGTATGGGGTGACCCTTGAGGACCTGGCTGAGCGCTATGAAGAAAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGTGCTGACCCAGCTCCAg                                                                                                                                      <  7:267210/149‑1 (MQ=255)
              cTTCGATGCTGAGTATGGGGTGACCCTTGAGGACCTGGCTGAGCGCTATGAAGAAAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGTGCTGACCCAGCTCCAgg                                                                                                                                     <  5:133389/149‑1 (MQ=255)
              cTTCGATGCGGAGTATGGGGTGACCCTTGAGGACCTGGCTGAGCGCTATGAAGAAAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGTGCTGACCCAGCTCCAgg                                                                                                                                     <  5:366297/149‑1 (MQ=255)
               tgcGATGATGAGTATGCGGTGACCCTTGAGGACCTGGCTGAGCGCTATGAAGAAAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGACGAAAGCTGCCATTAGGCATTGTCGCAAGTTAACCTGGACCTGGGTCAGCACCAgg                                                                                                                                     >  6:210825/3‑148 (MQ=255)
                    tGCTGAGTATGGGGGGACCCTTGAGGACCTGGCTGAGCGCTATGAAGAAAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGccc                                                                                                                               <  7:49866/149‑1 (MQ=255)
                        gAGTATGGGGTGACCCTTGAGGACCTGGCTGAGCGCTATGAAGAAAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTtgt                                                                                                                           >  5:262432/1‑149 (MQ=255)
                         aGTATGGGGTGACCCTTGAGGACCTGGCTGAGCGCTATGAAGCAAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTtgt                                                                                                                           <  5:38649/148‑1 (MQ=255)
                         aGTATGGGGTGACCCTTGAGGACCTGGCTGAGCGCTATGAAGAAAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTtgt                                                                                                                           <  5:70763/148‑1 (MQ=255)
                          gTAGGGGGTGACCCTTGAGGACCTGGCTGAGCGCTATGAAGAAAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTgg                                                                                                                         <  2:122346/149‑1 (MQ=255)
                           tATGGGGTGACCCTTGAGGACCTGGCTGAGCGCTATGAAGAAAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGa                                                                                                                        >  5:249972/1‑149 (MQ=255)
                           tATGGGGTGACCCTTGAGGACCTGGCTGAGCGCTATGAAGAAAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGa                                                                                                                        >  6:221944/1‑149 (MQ=255)
                            aTGGGGTGACCCTTGAGGACCTGGCTGAGCGCTATGAAGAAAGGCTGGCTAGCAACTTCCAGC‑AGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGAtt                                                                                                                      >  8:202517/1‑149 (MQ=255)
                               gggTGACCCTTGAGGACCTGGCTGAGCGCTATGAAGAAAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGAttt                                                                                                                     >  6:349834/1‑148 (MQ=255)
                                ggTGACCCTTGAGGACCTGGCTGAGCGCTATGAAGAAAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGa                                                                                                                   <  3:320909/149‑1 (MQ=255)
                                ggTGACCCTTGAGGACCTGGCTGAGCGCTATGAAGAAAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGa                                                                                                                   <  6:249972/149‑1 (MQ=255)
                                 gTGACCCTTGAGGACCTGGCTGAGCGCTATGAAGAAAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAg                                                                                                                  >  2:248254/1‑149 (MQ=255)
                                      ccTTGAGGACCTGGCTGAGCGCTATGAAGAAAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGTCCTTGTGGATTTGAG‑CCAg                                                                                                             >  6:118169/1‑148 (MQ=255)
                                          gAGGACCTGGCTGAGCGCTATGAAGAAAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGaa                                                                                                          <  5:221944/147‑1 (MQ=255)
                                            ggACCTGGCTGAGCGCTATGAAGAAAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTc                                                                                                       <  5:77363/148‑1 (MQ=255)
                                             gACCTGGCTGAGCGCTATGAAGAAAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAg                                                                                                                                      >  4:371140/1‑117 (MQ=255)
                                             gACCTGGCTGAGCGCTATGAAGAAAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAg                                                                                                                                      <  3:371140/117‑1 (MQ=255)
                                                     tGAGCGCTATGAAGAAAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTCGAGCGTCATc                                                                                             >  4:100018/1‑149 (MQ=255)
                                                       aGCGCTATGAAGAAAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTCGAGCGTCATcgc                                                                                           >  1:89417/1‑149 (MQ=255)
                                                        gcgcTATGAAGAAAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTAGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTAGTGGATTTAAC‑CCAGGAAATCGATCGtactagct                                                                                          >  8:121139/1‑142 (MQ=255)
                                                             aTGAAGAAAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTCGa                                                                                                     <  7:11053/133‑1 (MQ=255)
                                                             aTGAAGAAAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTCGa                                                                                                     >  8:11053/1‑133 (MQ=255)
                                                              tGAAGAAAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTAGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTCGAGc                                                                                                   >  8:25475/1‑134 (MQ=255)
                                                              tGAAGAAAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTAGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTCGAGc                                                                                                   <  7:25475/134‑1 (MQ=255)
                                                                aagaaAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTCGAGCGTCATCGCGCCGAACt                                                                                   <  5:70996/148‑1 (MQ=255)
                                                                aagaaAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTCGAGCGTCATCGCGCCGAACTg                                                                                  >  6:174829/1‑149 (MQ=255)
                                                                aagaaAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTCGAGCGTCATCGCGCCGAACTg                                                                                  >  4:185727/1‑149 (MQ=255)
                                                                aagaaAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTCGAGCGTCATCGCGCCGAACTg                                                                                  <  6:262432/149‑1 (MQ=255)
                                                                 agaaAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTCGAGCGTCATCGCGCCGAACTgg                                                                                 <  3:369579/149‑1 (MQ=255)
                                                                 agaaAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTCGAGCGTCATCGCGCCGAACTg                                                                                  >  2:197409/1‑148 (MQ=255)
                                                                  gaaAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTCGAGCGTCATCGCGCCGAACTg                                                                                  >  2:293382/1‑147 (MQ=255)
                                                                    aaGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTCGAGCGTCATCGCGCCGAACTGgaga                                                                              <  5:349834/149‑1 (MQ=255)
                                                                       gctggctAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTCGAGCGTCATCGCGCCGAACTGGAGAGcc                                                                           <  6:9056/149‑1 (MQ=255)
                                                                       gctggctAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTCGAGCGTCATCGCGCCGAACTGGAGAGcc                                                                           <  6:14287/149‑1 (MQ=255)
                                                                       gctggctAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑‑CAGGAAGTCGAGCGTCATCGCGCCGAACTGGAGAGCCt                                                                          >  1:203251/1‑149 (MQ=255)
                                                                           gctAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTCGAGCGTCATCGCGCCGAACTGGAGAGCCTcaa                                                                       <  7:301726/149‑1 (MQ=255)
                                                                                 aaCTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTCGAGCGTCATCGCGCCGAACTGGAGAGCCTCAAGCAAgt                                                                 <  3:100018/149‑1 (MQ=255)
                                                                                  agttCCAGCAAGTGGGCGACAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGGTGGGTCCGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTCGAGCGTCATCGCGCCGAACTGGAGAGCCTCAAGCAAgtg                                                                <  5:210825/147‑1 (MQ=255)
                                                                                   cTTCCATCAAGTTGTCGAAAGCGTCGAATGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCTGGCAGGCGGTGTGGATTTGAG‑CCAGGAAGTCGAGCGTCATCGCGCCGAACTGGAGATCCTCAAGCAAgtgg                                                               <  5:11253/149‑1 (MQ=255)
                                                                                      ccAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTCGAGCGTCATCGCGCCGAACTGGAGAGCCTCAAGCAAgtggt                                                              <  5:120504/147‑1 (MQ=255)
                                                                                            aGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAGCCCAGGAAGTCGAGCGTCATCGCGCCGAACTGGAGAGCCTCAAGCAAGTGGTGAACGGc                                                       <  1:293382/149‑1 (MQ=255)
                                                                                                 gCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTCGAGCGTCATCGCGCCGAACTGGAGAGCCTCAAGCAAGTGGTGAACGGCCTGAac                                                 <  2:163297/149‑1 (MQ=255)
                                                                                                   gAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTCGAGCGTCATCGCGCCGAACTGGAGAGCCTCAAGCAAGTGGTGAACGGCCTGAacca                                               <  1:242050/149‑1 (MQ=255)
                                                                                                   gAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTCGAGCGTCATCGCGCCGAACTGGAGAGCCTCAAGCAAGTGGTGAACGGCCTGAacca                                               <  3:233850/149‑1 (MQ=255)
                                                                                                     aaGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTCGAGCGTCATCGCGCCGAACTGGAGAGCCTCAAGCAAGTGGTGAACGGCCTGAaccacc                                             <  1:248254/149‑1 (MQ=255)
                                                                                                          tCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTCGAGCGTCATCGCGCCGAACTGGAGAGCCTCAAGCAAGTGGTGAACGGCCTGAACCACCATCTg                                        >  4:365645/1‑149 (MQ=255)
                                                                                                                gCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTCGAGCGTCATCGCGCCGAACTGGAGAGCCTCAAGCAAGTGGTGAACGGCCTGAACCACCATCTGTCCGAt                                  <  7:216810/149‑1 (MQ=255)
                                                                                                                 cTGACATTCTGCATTCTCGAAAGTTAACCTTGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTCGAGCGTCATCGCGCCGAACTGGAGAGCCTCAAGCAAGTGGTGAACGGCCTGAACCACCATCTGTCCGATa                                 <  7:121139/149‑1 (MQ=255)
                                                                                                                  tGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTCGAGCGTCATCGCGCCGAACTGGAGAGCCTCAAGCAAGTGGTGAACGGCCTGAACCACCATCTGTCCGATAc                                >  7:313026/1‑149 (MQ=255)
                                                                                                                     cATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑‑CAGGAAGTCGAGCGTCATCGCGCCGAACTGGAGAGCCTCAAGCAAGTGGTGAACGGCCTGAACCACCATCTGTCCGATACCCA‑g                           <  2:203251/149‑1 (MQ=255)
                                                                                                                      aTTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTCGAGCGTCATCGCGCCGAACTGGAGAGCCTCAAGCAAGTGGTGAACGGCCTGAACCACCATCTGTCCGATACCCA‑g                           <  5:174829/149‑1 (MQ=255)
                                                                                                                           gCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTCGAGCGTCATCGCGCCGAACTGGAGAGCCTCAAGCAAGTGGTGAACGGCCTGAACCACCATCTGTCCGATACCCA‑GCAACg                      >  1:169567/1‑149 (MQ=255)
                                                                                                                           gCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTCGAGCGTCATCGCGCCGAACTGGAGAGCCTCAAGCAAGTGGTGAACGGCCTGAACCACCATCTGTCCGATACCCA‑GCAACg                      >  7:21103/1‑149 (MQ=255)
                                                                                                                                   gCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTCGAGCGTCATCGCGCCGAACTGGAGAGCCTCAAGCAAGTGGTGAACGGCCTGAACCACCATCTGTCCGATACCCA‑GCAACGGCTGGctc              >  2:323737/1‑149 (MQ=255)
                                                                                                                                       gTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTCGAGCGTCATCGCGCCGAACTGGAGAGCCTCAAGCAAGTGGTGAACGGCCTGAACCACCATCTGTCCGATACCCA‑GCAACGGCTGGCTCTTGc          >  2:48481/1‑149 (MQ=255)
                                                                                                                                          tgccTGGTGCTGACCCAGCTCCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTCGAGCGTCATCGCGCCGAACTGGAGAGCCTCAAGCAAGTGGTGAACGGCCTGAACCACCATCTGTCCGATACCCA‑GCAACGGCTGGCTCTTGCtga       >  7:199847/3‑149 (MQ=255)
                                                                                                                                             cTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CAAGGAAGTCGAGCGTCATCGCGCCGAACTGGAGAGCCTCAAGCAAGTGGTGAACGGCCTGAACCCCCATCTGTCCGATACCCA‑GCAACGGATGGCTCTTGCtgat      >  6:253865/1‑147 (MQ=255)
                                                                                                                                                gAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTCGAGCGTCATCGCGCCGAACTGGAGAGCCTCAAGCAAGTGGTGAACGGCCTGAACCACCATCTGTCCGATACCCACGCAACGGCTGGCTCTTGCTGATGAGc  <  1:197409/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GCCTTCGCCAGCGCCTTCGATGCTGAGTATGGGGTGACCCTTGAGGACCTGGCTGAGCGCTATGAAGAAAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTCGAGCGTCATCGCGCCGAACTGGAGAGCCTCAAGCAAGTGGTGAACGGCCTGAACCACCATCTGTCCGATACCCA‑GCAACGGCTGGCTCTTGCTGATGAGC  >  NC_002947/1989420‑1989711

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