Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,270,641 C→A 14.1% S419S (TCC→TCA 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,6410CA14.1% 456.4 / 54.7 163S419S (TCC→TCAPP_4643xanthine/uracil permease family protein
Reads supporting (aligned to +/- strand):  ref base C (82/57);  new base A (11/12);  total (94/69)
Fisher's exact test for biased strand distribution p-value = 3.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 167 total aligned reads displayed.

AGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACATTTCCGCGCTGCGTGATGGG  >  NC_002947/5270492‑5270788
                                                                                                                                                     |                                                                                                                                                   
aGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTcc                                                                                                                                                     >  2:2254966/1‑150 (MQ=255)
aGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTcc                                                                                                                                                     <  1:1291679/150‑1 (MQ=255)
aGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTcc                                                                                                                                                     >  1:4835077/1‑150 (MQ=255)
aGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTcc                                                                                                                                                     >  1:913605/1‑150 (MQ=255)
    gTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCa                                                                                                                                                 <  2:1668128/150‑1 (MQ=255)
      aTCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCa                                                                                                                                              <  1:5072919/150‑1 (MQ=255)
      aTCCGCACACTTTCAAAGGTCGACTATCAGAACAACATCAATCTGATCATTGTGGCTACCTCCATCGGCTTCGGCATGATCCCGATTGCAGCACCAGAATTCTACCATCAGTTCCCGGCCTGGTTTGAAACCATTTTTCATTCAGGCATc                                                                                                                                               >  2:2140693/1‑150 (MQ=255)
      aTCCGCACACTTTCAAAGGTCGACTATCAGAACAACATCAATCTGATCATTGTGGCTACCTCCATCGGCTTCGGCATGATCCCGATTGCAGCACCAGAATTCTACCATCAGTTCCCGGCCTGGTTTGAAACCATTTTTCATTCAGGCATc                                                                                                                                               >  2:2140708/1‑150 (MQ=255)
       tCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCa                                                                                                                                              <  2:181953/150‑1 (MQ=255)
        ccGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAg                                                                                                                                             >  2:1881236/1‑150 (MQ=255)
         cGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGc                                                                                                                                            >  1:4007818/1‑150 (MQ=255)
           cACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTc                                                                                                                                          <  1:552977/150‑1 (MQ=255)
              ccTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTcggc                                                                                                                                       >  1:3746517/1‑150 (MQ=255)
              ccTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGACCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTcggc                                                                                                                                       <  1:713279/150‑1 (MQ=255)
               ctttcaaaGGTCGACTATCAGAACAACATCAATCTGATCATTGTGGCTACCTCCATCGGCTTCGGCATGATCCCGATTGCAGCACCAGAATTCTACCATCAGTTCCCGGCCTGGTTTGAAACCATTTTTCATTCAGGCATCAGCTccgct                                                                                                                                      <  1:2532048/144‑2 (MQ=255)
                 ttcaaaGGTCGACTATCAGAACAACATCAATCTGATCATTGTGGCTACCTCCATCGGCTTCGGCATGATCCCGATTGCAGCACCAGAATTCTACCATCAGTTCCCGGCCTGGTTTGAAACCATTTTTCATTCAGGCATCAGCTCcgctgc                                                                                                                                    >  2:1139876/5‑150 (MQ=255)
                    cAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGag                                                                                                                                 >  2:3028709/1‑149 (MQ=255)
                    cAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGAt                                                                                                                                 >  1:1903880/1‑150 (MQ=255)
                    aaaGGTCGACTATCAGAACAACATCAATCTGATCATTGTGGCTACCTCCATCGGCTTCGGCATGATCCCGATTGCAGCACCAGAATTCTACCATCAGTTCCCGGCCTGGTTTGAAACCATTTTTCATTCAGGCATCAGCTCcgctgcatt                                                                                                                                 >  1:3846870/2‑147 (MQ=255)
                     aaGGTCGACTATCAGAACAACATCAATCTGATCATTGTGGCTACCTCCATCGGCTTCGGCATGATCCCGATTGCAGCACCAGAATTCTACCATCAGTTCCCGGCCTGGTTTGAAACCATTTTTCATTCAGGCATCAGCTCcgctgcattg                                                                                                                                <  2:164160/150‑5 (MQ=255)
                     aaGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATc                                                                                                                                <  2:4856913/150‑1 (MQ=255)
                        gTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATg                                                                                                                             <  1:114966/150‑1 (MQ=255)
                         tCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATgg                                                                                                                            <  1:2012364/150‑1 (MQ=255)
                         tCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATgg                                                                                                                            <  2:533379/150‑1 (MQ=255)
                         tCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATgg                                                                                                                            <  2:2029408/150‑1 (MQ=255)
                             cTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCAt                                                                                                                        >  1:4167392/1‑150 (MQ=255)
                               aCAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATcc                                                                                                                      >  2:3615295/1‑150 (MQ=255)
                               aCAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATcc                                                                                                                      >  1:2057203/1‑150 (MQ=255)
                                  aGAACAACATCAATCTGATCATTGTGGCTACCTCCATCGGCTTCGGCATGATCCCGATTGCAGCACCAGAATTCTACCATCAGTTCCCGGCCTGGTTTGAAACCATTTTTCATTCAGGCATCAGCTCCGCTGCATTGATGGCAATcgcgc                                                                                                                   <  1:2140693/150‑5 (MQ=255)
                                  aGAACAACATCAATCTGATCATTGTGGCTACCTCCATCGGCTTCGGCATGATCCCGATTGCAGCACCAGAATTCTACCATCAGTTCCCGGCCTGGTTTGAAACCATTTTTCATTCAGGCATCAGCTCCGCTGCATTGATGGCAATcgcgc                                                                                                                   <  1:2140708/150‑5 (MQ=255)
                                     acaacGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGa                                                                                                                <  1:828671/150‑1 (MQ=255)
                                        acGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAAcc                                                                                                             <  2:3957203/150‑1 (MQ=255)
                                        acGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAAcc                                                                                                             <  1:4840785/150‑1 (MQ=255)
                                           tCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGa                                                                                                          >  2:1683453/1‑150 (MQ=255)
                                            cAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGAt                                                                                                         >  2:2317122/1‑150 (MQ=255)
                                             aaCCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATc                                                                                                        >  1:765057/1‑150 (MQ=255)
                                                 tGATCATTGTGGCTACCTCCATCGGCTTCGGCATGATCCCGATTGCAGCACCAGAATTCTACCATCAGTTCCCGGCCTGGTTTGAAACCATTTTTCATTCAGGCATCAGCTCCGCTGCATTGATGGCAATCGCGCTCAATCTCATGTTCa                                                                                                    >  1:524326/1‑150 (MQ=255)
                                                 tGATCATTGTGGCTACCTCCATCGGCTTCGGCATGATCCCGATTGCAGCACCAGAATTCTACCATCAGTTCCCGGCCTGGTTTGAAACCATTTTTCATTCAGGCATCAGCTCCGCTGCATTGATGGCAATCGCGCTCAATCTCATGTTCa                                                                                                    >  2:2904352/1‑150 (MQ=255)
                                                 tGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCa                                                                                                    >  2:1386746/1‑150 (MQ=255)
                                                 tGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGACCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCa                                                                                                    <  2:722892/150‑1 (MQ=255)
                                                 tGATCATCGCCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCa                                                                                                    <  1:3078681/150‑1 (MQ=255)
                                                 agATCATTGTGGCTACCTCCATCGGCTTCGGCATGATCCCGATTGCAGCACCAGAATTCTACCATCAGTTCCCGGCCTGGTTTGAAACCATTTTTCATTCAGGCATCAGCTCCGCTGCATTGATGGCAATCGCGCTCAATCTCATGTTCa                                                                                                    <  2:4961541/149‑1 (MQ=255)
                                                  gATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCaa                                                                                                   <  1:4526020/150‑1 (MQ=255)
                                                  gATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGAGTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCaa                                                                                                   <  2:3713892/150‑1 (MQ=255)
                                                      atcGTCGCGGCCTCGCTGGGCTTCGGCATGATACCGATTGCCGCACCGACCTTCAACCACCACATCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCAt                                                                                               >  1:4325821/1‑150 (MQ=255)
                                                         gtcgCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTc                                                                                            >  1:2751369/1‑150 (MQ=255)
                                                         gtcgCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTc                                                                                            <  2:1997780/150‑1 (MQ=255)
                                                              ggCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTcaccac                                                                                       <  2:1202518/150‑1 (MQ=255)
                                                                ccTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCg                                                                                     <  1:43450/150‑1 (MQ=255)
                                                                  tCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGc                                                                                   >  2:3830751/1‑150 (MQ=255)
                                                                      tGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCAGAACTGGTTCGAGACCATCTTCCACACCGGCATCAGATCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACt                                                                               >  1:1909179/1‑150 (MQ=255)
                                                                       gggCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTc                                                                              <  2:2106666/150‑1 (MQ=255)
                                                                        ggCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCg                                                                             <  1:4661375/150‑1 (MQ=255)
                                                                          cTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGa                                                                           <  1:1159488/150‑1 (MQ=255)
                                                                          cTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGa                                                                           >  1:3430165/1‑150 (MQ=255)
                                                                          cTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGa                                                                           <  1:3818725/150‑1 (MQ=255)
                                                                          cTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAAACTGATCTTCAACCATTTCACCACCGGCAACTCGGa                                                                           >  2:4534776/1‑150 (MQ=255)
                                                                               gCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAAcc                                                                      >  1:2896391/1‑150 (MQ=255)
                                                                                cATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACca                                                                     <  2:4635784/150‑1 (MQ=255)
                                                                                  tGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACcagt                                                                   <  2:1651855/150‑1 (MQ=255)
                                                                                    aTCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACcagtca                                                                 <  1:4488339/150‑1 (MQ=255)
                                                                                              ccGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTg                                                       >  2:2800126/1‑150 (MQ=255)
                                                                                                 cACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCCTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAg                                                    >  2:155635/1‑150 (MQ=255)
                                                                                                        cTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGa                                             <  2:3890207/150‑1 (MQ=255)
                                                                                                        cTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGa                                             <  2:2460541/150‑1 (MQ=255)
                                                                                                           cTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAgcg                                          >  1:4888087/1‑150 (MQ=255)
                                                                                                            tcgaaccagTTCGGGAAGTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAgcgc                                         >  2:5034191/5‑150 (MQ=255)
                                                                                                            tACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAgcgc                                         <  1:3188399/150‑1 (MQ=255)
                                                                                                              ccaccaCTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCAc                                       >  2:1898582/1‑150 (MQ=255)
                                                                                                                   acTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATcc                                  >  1:2117387/1‑150 (MQ=255)
                                                                                                                    cTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCa                                 <  1:4709020/150‑1 (MQ=255)
                                                                                                                    cTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCCGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCa                                 <  2:4835077/150‑1 (MQ=255)
                                                                                                                    cTTCCAGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCa                                 <  2:1633209/150‑1 (MQ=255)
                                                                                                                     ttCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAg                                <  2:913605/150‑1 (MQ=255)
                                                                                                                        ccGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACcagtcagt                                                               <  2:1946595/116‑1 (MQ=255)
                                                                                                                              tGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCg                                                          <  2:76648/115‑1 (MQ=255)
                                                                                                                              tGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGAc                       >  1:3287275/1‑150 (MQ=255)
                                                                                                                              tGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGAc                       >  1:3948989/1‑150 (MQ=255)
                                                                                                                              tGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGAc                       >  2:4512286/1‑150 (MQ=255)
                                                                                                                              tGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGAc                       >  1:721411/1‑150 (MQ=255)
                                                                                                                                 ttCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACAtt                    >  1:3367172/1‑150 (MQ=255)
                                                                                                                                  tCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACAttt                   >  2:5008032/1‑150 (MQ=255)
                                                                                                                                    gagaCCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACATTTcc                 <  2:3173480/150‑1 (MQ=255)
                                                                                                                                    gagaCCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACATTTcc                 >  1:1732209/1‑150 (MQ=255)
                                                                                                                                     agaCCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACATTTCcg                >  2:1198114/1‑150 (MQ=255)
                                                                                                                                      gaCCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACATTTCcgc               >  1:2963022/1‑150 (MQ=255)
                                                                                                                                        ccATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACATTTCCgcgc             <  2:1339580/150‑1 (MQ=255)
                                                                                                                                         cATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCACCTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACATTTCCGCGCt            >  2:3299761/1‑150 (MQ=255)
                                                                                                                                         cATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACATTTCCGCGCt            >  1:165774/1‑150 (MQ=255)
                                                                                                                                         cATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACATTTCCGCGCt            >  2:82855/1‑150 (MQ=255)
                                                                                                                                         cATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACATTTCCGCGCt            >  2:2589930/1‑150 (MQ=255)
                                                                                                                                         cATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACATTTCCGCGCt            >  1:165773/1‑150 (MQ=255)
                                                                                                                                         cATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACATTTCCGCGCt            >  2:2589892/1‑150 (MQ=255)
                                                                                                                                           tCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACATTTCCGCGCTGc          >  2:1369354/1‑150 (MQ=255)
                                                                                                                                             ttCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACATTTCCGCGCTGCGt        >  1:17467/1‑150 (MQ=255)
                                                                                                                                             ttCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACATTTCCGCGCTGCGt        <  2:3110636/150‑1 (MQ=255)
                                                                                                                                               ccACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACATTTCCGCGCTGCGTGa      >  2:3615515/1‑150 (MQ=255)
                                                                                                                                                   tCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACATTTCCGCGCTGCGTGATggg  >  1:2829021/1‑150 (MQ=255)
                                                                                                                                                   tCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACATTTCCGCGCTGCGTGATggg  <  2:2285501/150‑1 (MQ=255)
                                                                                                                                                   tCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACATTTCCGCGCTGCGTGATggg  <  2:2896391/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
AGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACATTTCCGCGCTGCGTGATGGG  >  NC_002947/5270492‑5270788

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: