Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,270,525 A→C 22.0% K381Q (AAG→CAG)  PP_4643 → xanthine/uracil permease family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,270,5250AC22.0% 691.5 / 207.2 322K381Q (AAG→CAG) PP_4643xanthine/uracil permease family protein
Reads supporting (aligned to +/- strand):  ref base A (110/140);  new base C (38/33);  total (148/174)
Fisher's exact test for biased strand distribution p-value = 1.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.39e-01

Only 100 of 333 total aligned reads displayed.

GCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTG  >  NC_002947/5270376‑5270674
                                                                                                                                                     |                                                                                                                                                     
 cggcgTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACaa                                                                                                                                                      <  2:1094007/150‑1 (MQ=255)
 cggcgTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACaa                                                                                                                                                      <  2:5388153/150‑1 (MQ=255)
   gcgTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGa                                                                                                                                                    <  2:8121534/150‑1 (MQ=255)
     gTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGaac                                                                                                                                                  <  2:84718/150‑1 (MQ=255)
       gATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGaacaa                                                                                                                                                <  2:2666334/150‑1 (MQ=255)
        tttCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCTAGGTCAGCTTCAAGaacaac                                                                                                                                               <  2:5275879/149‑1 (MQ=255)
        aTTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGTGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGaacaac                                                                                                                                               <  2:1241314/150‑1 (MQ=255)
         ttCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACg                                                                                                                                              >  1:6373684/1‑150 (MQ=255)
         ttCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACg                                                                                                                                              >  2:4934609/1‑150 (MQ=255)
            tggtggTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCa                                                                                                                                           >  1:4079208/1‑150 (MQ=255)
             ggtggtACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCaa                                                                                                                                          >  2:6543840/1‑150 (MQ=255)
              gtggtACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAAc                                                                                                                                         <  2:3122580/150‑1 (MQ=255)
               tggtACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAAcc                                                                                                                                        >  2:1455369/1‑150 (MQ=255)
               tggtACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAAcc                                                                                                                                        >  2:4175678/1‑150 (MQ=255)
                   aCTCGGTCTTCTGCCTGTAATGGGACGGGTGATTGCCTGTGTTCCGCCAGCAGTGCTGGGTGGTGCCGGTATCGTGCTGTTCGGAACGGTTGCCGCAAGCGGCATCCGCACACTTTCAAAGGTCGACTATCAGAACAACATCAATCTGat                                                                                                                                    >  1:1110305/1‑150 (MQ=255)
                   aCTCGGTCTTCTGCCTGTAATGGGACGGGTGATTGCCTGTGTTCCGCCAGCAGTGCTGGGTGGAGCCGGTATCGTGCTGTTCGGAACGGTTGCCGCAAGCGGCATCCGCACACTTTCAAAGGTCGACTATCAGAACAACATCAATCTGat                                                                                                                                    >  2:3036681/1‑150 (MQ=255)
                   aCTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGat                                                                                                                                    <  2:4606985/150‑1 (MQ=255)
                    cTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGatc                                                                                                                                   <  2:7161641/150‑1 (MQ=255)
                     tCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGatca                                                                                                                                  >  2:1684991/1‑150 (MQ=255)
                       ggtcttctgcCTGTAATGGGACGGGTGATTGCCTGTGTTCCGCCAGCAGTGCTGGGTGGTGCCGGTATCGTGCTGTTCGGAACGGTTGCCGCAAGCGGCATCCGCACACTTTCAAAGGTCGACTATCAGAACAACATCAATCTGatcatt                                                                                                                                >  1:3492775/7‑149 (MQ=255)
                       ggtcttctgcCTGTAATGGGACGGGTGATTGCCTGTGTTCCGCCAGCAGTGCTGGGTGGTGCCGGTATCGTGCTGTTCGGAACGGTTGCCGCAAGCGGCATCCGCACACTTTCAAAGGTCGACTATCAGAACAACATCAATCTGatcatt                                                                                                                                >  1:2900881/7‑149 (MQ=255)
                       ggtcttctgcCTGTAATGGGACGGGTGATTGCCTGTGTTCCGCCAGCAGTGCTGGGTGGTGCCGGTATCGTGCTGTTCGGAACGGTTGCCGCAAGCGGCATCCGCACACTTTCAAAGGTCGACTATCAGAACAACATCAATCTGatcatt                                                                                                                                >  1:7519495/7‑149 (MQ=255)
                       ggtcttctgcCTGTAATGGGACGGGTGATTGCCTGTGTTCCGCCAGCAGTGCTGGGTGGTGCCGGTATCGTGCTGTTCGGAACGGTTGCCGCAAGCGGCATCCGCACACTTTCAAAGGTCGACTATCAGAACAACATCAATCTGatcatt                                                                                                                                >  1:8282139/7‑149 (MQ=255)
                       ggCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGatcatc                                                                                                                                >  1:5433470/1‑150 (MQ=255)
                            tctgcCTGTAATGGGACGGGTGATTGCCTGTGTTCCGCCAGCAGTGCTGGGTGGTGCCGGTATCGTGCTGTTCGGAACGGTTGCCGCAAGCGGCATCCGCACACTTTCAAAGGTCGACTATCAGAACAACATCAATCTGATCATTGTggc                                                                                                                           <  1:3316004/149‑1 (MQ=255)
                              tgcCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGaacaac                                                                                                                                               >  2:691725/1‑128 (MQ=255)
                              tgcCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCgg                                                                                                                         >  1:6029282/1‑150 (MQ=255)
                               gcCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGGGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTCTCCAAGGTCAGCGACAAGAACAACGTCAACCTGATCCTCGTCGCGGc                                                                                                                        >  1:4197389/1‑150 (MQ=255)
                                ccAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGcc                                                                                                                       >  1:883191/1‑150 (MQ=255)
                                  aGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTc                                                                                                                     <  1:6894418/150‑1 (MQ=255)
                                    tCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGc                                                                                                                   <  1:2161561/150‑1 (MQ=255)
                                     cATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCt                                                                                                                  <  1:6467813/150‑1 (MQ=255)
                                      aTGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTg                                                                                                                 <  1:5568202/150‑1 (MQ=255)
                                       tGGGACGGGTGATTGCCTGTGTTCCGCCAGCAGTGCTGGGTGGTGCCGGTATCGTGCTGTTCGGAACGGTTGCCGCAAGCGGCATCCGCACACTTTCAAAGGTCGACTATCAGAACAACATCAATCTGATCATTGTGGCTACCTccatcg                                                                                                                >  1:1753892/1‑145 (MQ=255)
                                        gggCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTggg                                                                                                               <  1:1529976/150‑1 (MQ=255)
                                              ggTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCgg                                                                                                         >  1:1854326/1‑150 (MQ=255)
                                              ggTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCgg                                                                                                         <  2:2896982/150‑1 (MQ=255)
                                              ggTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCgg                                                                                                         <  1:8418522/150‑1 (MQ=255)
                                               gTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGc                                                                                                        <  2:4801020/150‑1 (MQ=255)
                                               gTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGc                                                                                                        <  1:9077737/150‑1 (MQ=255)
                                                tGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCa                                                                                                       >  2:8845261/1‑150 (MQ=255)
                                                     gccgccGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATc                                                                                                  <  2:4779383/150‑1 (MQ=255)
                                                        gcagTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCg                                                                                               >  2:7857095/4‑150 (MQ=255)
                                                           gTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGAtt                                                                                            <  1:1277926/150‑1 (MQ=255)
                                                              ccTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCTGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGcc                                                                                         <  2:1715845/150‑1 (MQ=255)
                                                                tACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGc                                                                                       <  2:698665/150‑1 (MQ=255)
                                                                   cccGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCAcc                                                                                    >  2:1186969/1‑150 (MQ=255)
                                                                    gcagtgcTGGGTGGTGCCGGTATCGTGCTGTTCGGAACGGTTGCCGCAAGCGGCATCCGCACACTTTCAAAGGTCGACTATCAGAACAACATCAATCTGATCATTGTGGCTACCTCCATCGGCTTCGGCATGATCCCGATTGCAGCAcca                                                                                   <  2:6791837/144‑2 (MQ=255)
                                                                      ggTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGAc                                                                                 >  2:8921404/1‑150 (MQ=255)
                                                                         aCTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCtt                                                                              <  2:2006002/150‑1 (MQ=255)
                                                                                 gtgcCGGTATCGTGCTGTTCGGAACGGTTGCCGCAAGCGGCATCCGCACACTTTCAAAGGTCGACTATCAGAACAACATCAATCTGATCATTGTGGCTACCTCCATCGGCTTCGGCATGATCCCGATTGCAGCACCAGAATTCTaccatc                                                                      <  2:30944/148‑3 (MQ=255)
                                                                                   gcCGGTATCGTGCTGTTCGGAACGGTTGCCGCAAGCGGCATCCGCACACTTTCAAAGGTCGACTATCAGAACAACATCAATCTGATCATTGTGGCTACCTCCATCGGCTTCGGCATGATCCCGa                                                                                              >  2:1633484/1‑124 (MQ=255)
                                                                                     cGGTATCGTGCTGTTCGGAACGGTTGCCGCAAGCGGCATCCGCACACTTTCAAAGGTCGACTATCAGAACAACATCAATCTGATCATTGTGGCTACCTCCATCGGCTTCGGCATGATCCCGATTGCAGCACCAGAATTCTACCATCAGtt                                                                  >  2:3734231/1‑150 (MQ=255)
                                                                                      ggTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATccc                                                                                                >  1:5097214/1‑119 (MQ=255)
                                                                                      ggTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTc                                                                 <  2:7428412/150‑1 (MQ=255)
                                                                                        tATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTccc                                                               >  1:1022431/1‑150 (MQ=255)
                                                                                            gTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAAc                                                           >  1:3065410/1‑150 (MQ=255)
                                                                                            gTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAAc                                                           <  2:3270143/150‑1 (MQ=255)
                                                                                             tGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACt                                                          >  2:5879327/1‑150 (MQ=255)
                                                                                               cTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTgg                                                        <  1:228151/150‑1 (MQ=255)
                                                                                               cTGTTCGGAACGGTTGCCGCAAGCGGCATCCGCACACTTTCAAAGGTCGACTATCAGAACAACATCAATCTGATCATTGTGGCTACCTCCATCGGCTTCGGCATGATCCCGATTGCAGCACCAGAATTCTACCATCAGTTCCCGGCCTgg                                                        <  2:5802206/150‑1 (MQ=255)
                                                                                                  ttCGGAACGGTTGCCGCAAGCGGCATCCGCACACTTTCAAAGGTCGACTATCAGAACAACATCAATCTGATCATTGTGGCTACCTCCATCGGCTTCGGCATGATCCCGATTGCAGCACCAGAATTCTACCATCAGTTCCCGGCCTGGttt                                                     >  1:7474528/1‑149 (MQ=255)
                                                                                                  ttCGGAACGGTTGCCGCAAGCGGCATCCGCACACTTTCAAAGGTCGACTATCAGAACAACATCAATCTGATCATTGTGGCTACCTCCATCGGCTTCGGCATGATCCCGATTGCAGCACCAGAATTCTACCATCAGTTCCCGGCCTGGttt                                                     >  2:6034721/1‑149 (MQ=255)
                                                                                                   tCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTaccac                                                                       <  2:2962258/131‑1 (MQ=255)
                                                                                                    cGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCga                                                   <  1:7799586/150‑1 (MQ=255)
                                                                                                    cGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCga                                                   <  1:1423975/150‑1 (MQ=255)
                                                                                                     ggCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCgag                                                  <  1:5664986/150‑1 (MQ=255)
                                                                                                      gCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCgaga                                                 <  1:1656010/150‑1 (MQ=255)
                                                                                                      gCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCgaga                                                 <  2:871375/150‑1 (MQ=255)
                                                                                                        aCCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGAcc                                               <  1:6164937/150‑1 (MQ=255)
                                                                                                         ccGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCa                                              >  2:1555730/1‑150 (MQ=255)
                                                                                                              gcggcCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTc                                         <  1:4939586/150‑1 (MQ=255)
                                                                                                              gccgcaaGCGGCATCCGCACACTTTCAAAGGTCGACTATCAGAACAACATCAATCTGATCATTGTGGCTACCTCCATCGGCTTCGGCATGATCCCGATTGCAGCACCAGAATTCTACCATCAGTTCCCGGCCTGGTTTGAAACCAttttt                                         >  1:2720661/7‑149 (MQ=255)
                                                                                                              gccgcaaGCGGCATCCGCACACTTTCAAAGGTCGACTATCAGAACAACATCAATCTGATCATTGTGGCTACCTCCATCGGCTTCGGCATGATCCCGATTGCAGCACCAGAATTCTACCATCAGTTCCCGGCCTGGTTTGAAACCAttttt                                         <  2:1891572/144‑2 (MQ=255)
                                                                                                                ggcCAGTGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCa                                       >  2:83275/1‑150 (MQ=255)
                                                                                                                ggcCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCa                                       >  2:7958175/1‑150 (MQ=255)
                                                                                                                 gcaaGCGGCATCCGCACACTTTCAAAGGTCGACTATCAGAACAACATCAATCTGATCATTGTGGCTACCTCCATCGGCTTCGGCATGATCCCGATTGCAGCACCAGAATTCTACCATCAGTTCCCGGCCTGGTTTGAAACCATTTTTCat                                      <  1:9165274/147‑2 (MQ=255)
                                                                                                                  ccAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACt                                     >  2:5232367/1‑150 (MQ=255)
                                                                                                                  ccAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACt                                     >  2:2225321/1‑150 (MQ=255)
                                                                                                                   ttgCGGCATCCGCACACTTTCAAAGGTCGACTATCAGAACAACATCAATCTGATCATTGTGGCTACCTCCATCGGCTTCGGCATGATCCCGATTGCAGCACCAGAATTCTACCATCAGTTCCCGGCCTGGTTTGAAACCATTTTTCATTc                                    >  2:8260144/3‑150 (MQ=255)
                                                                                                                   cAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTc                                    >  1:8252739/1‑150 (MQ=255)
                                                                                                                    aGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTcc                                   >  2:703986/1‑150 (MQ=255)
                                                                                                                    aGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTcc                                   >  2:5094220/1‑150 (MQ=255)
                                                                                                                     gCGGCATCCGCACACTTTCAAAGGTCGACTATCAGAACAACATCAATCTGATCATTGTGGCTACCTCCATCGGCTTCGGCATGATCCCGATTGCAGCACCAGAATTCTACCATCAGTTCCCGGCCTGGTTTGAAACCATTTTTCATTcag                                  >  1:5779883/1‑148 (MQ=255)
                                                                                                                     gCGGCATCCGCACACTTTCAAAGGTCGACTATCAGAACAACATCAACCTGATCATTGTGGCTACCTCCATCGGCTTCGGCATGATCCCGATTGCAGCACCAGAATTCTACCATCAGTTCCCGGCCTGGTTTGAAACCATTTTTCATTcag                                  <  2:4054548/150‑3 (MQ=255)
                                                                                                                          aTCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATc                             <  1:4809202/150‑1 (MQ=255)
                                                                                                                           tCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCa                            <  1:2009448/150‑1 (MQ=255)
                                                                                                                             cGCACACTTTCAAAGGTCGACTATCAGAACAACATCAATCTGATCATTGTGGCTACCTCCATCGGCTTCGGCATGATCCCGATTGCAGCACCAGAATTCTACCATCAGTTCCCGGCCTGGTTTGAAACCATTTTTCATTCAGGCATCAGc                          >  2:5545170/1‑150 (MQ=255)
                                                                                                                              gCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCt                         >  1:6565526/1‑150 (MQ=255)
                                                                                                                                  ccTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTcggc                     >  1:3852929/1‑150 (MQ=255)
                                                                                                                                   ctttcaaaGGTCGACTATCAGAACAACATCAATCTGATCATTGTGGCTACCTCCATCGGCTTCGGCATGATCCCGATTGCAGCACCAGAATTCTACCATCAGTTCCCGGCCTGGTTTGAAACCATTTTTCATTCAGGCATCAGCTccgct                    <  2:3492775/144‑2 (MQ=255)
                                                                                                                                    tGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTcggcgg                   >  1:2414016/1‑150 (MQ=255)
                                                                                                                                     ttcaaaGGTCGACTATCAGAACAACATCAATCTGATCATTGTGGCTACCTCCATCGTCTTCGGCATGATCCCGATTGCAGCACCAGAATTCTACCATCAGTTCCCGGCCTGGTTTGAAACCATTTTTCATTCAGGCATCAGCTCcgctgc                  <  1:2984157/146‑1 (MQ=255)
                                                                                                                                        cAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCGTGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGAt               >  1:4363519/1‑150 (MQ=255)
                                                                                                                                        aaaGGTCGACTATCAGAACAACATCAATCTGATCATTGTGGCTACCTCCATCGGCTTCGGCATGATCCCGATTGCAGCACCAGAATTCTACCATCAGTTCCCGGCCTGGTTTGAAACCATTTTTCATTCAGGCATCAGCTCcgctgcatt               >  2:615886/2‑147 (MQ=255)
                                                                                                                                          aGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCa             >  2:2378212/1‑150 (MQ=255)
                                                                                                                                                   aCAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATcc    >  2:6323860/1‑150 (MQ=255)
                                                                                                                                                    cAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCct   <  1:600631/150‑1 (MQ=255)
                                                                                                                                                     aaGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCctg  <  1:7036640/150‑1 (MQ=255)
                                                                                                                                                     aaGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCctg  >  2:4367020/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTG  >  NC_002947/5270376‑5270674

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: