Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,270,653 G→C 7.8% S423S (TCG→TCC 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,6530GC 7.8% 457.6 / 22.7 141S423S (TCG→TCCPP_4643xanthine/uracil permease family protein
Reads supporting (aligned to +/- strand):  ref base G (79/51);  new base C (5/6);  total (84/57)
Fisher's exact test for biased strand distribution p-value = 3.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01

Only 100 of 151 total aligned reads displayed.

CCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACATTTCCGCGCTGCGTGATGGGGACTACTTCAAGG  >  NC_002947/5270506‑5270801
                                                                                                                                                   |                                                                                                                                                    
actttcaaaGGTCGACTATCAGAACAACATCAATCTGATCATTGTGGCTACCTCCATCGGCTTCGGCATGATCCCGATTGCAGCACCAGAATTCTACCATCAGTTCCCGGCCTGGTTTGAAACCATTTTTCATTCAGGCATCAGCTccgc                                                                                                                                                    <  2:3908106/143‑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:391452/2‑147 (MQ=255)
       aaGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATc                                                                                                                                             <  2:4856913/150‑1 (MQ=255)
       aaGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATc                                                                                                                                             <  2:2798037/150‑1 (MQ=255)
          gTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATg                                                                                                                                          <  1:114966/150‑1 (MQ=255)
           tCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATgg                                                                                                                                         <  1:2012364/150‑1 (MQ=255)
           tCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATgg                                                                                                                                         <  2:2029408/150‑1 (MQ=255)
           tCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATgg                                                                                                                                         <  2:533379/150‑1 (MQ=255)
            cAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGc                                                                                                                                        >  2:4997670/1‑150 (MQ=255)
               cTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCAt                                                                                                                                     >  1:4167392/1‑150 (MQ=255)
                tACAAGTACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGAGCATGGCCATc                                                                                                                                    >  2:2486294/1‑150 (MQ=255)
                 aCAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATcc                                                                                                                                   >  2:3615295/1‑150 (MQ=255)
                  tcaGAACAACATCAATCTGATCATTGTGGCTACCTCCATCGGCTTCGGCATGATCCCGATTGCAGCACCAGAATTCTACCATCAGTTCCCGGCCTGGTTTGAAACCATTTTTCATTCAGGCATCAGCTCCGCTGCATTGATGGCAATcgc                                                                                                                                  >  1:732761/3‑148 (MQ=255)
                  tcaGAACAACATCAATCTGATCATTGTGGCTACCTCCATCGGCTTCGGCATGATCCCGATTGCAGCACCAGAATTCTACCATCAGTTCCCGGCCTGGTTTGAAACCATTTTTCATTCAGGCATCAGCTCCGCTGCATTGATGGCAATcgc                                                                                                                                  >  2:2916940/3‑148 (MQ=255)
                  cAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCct                                                                                                                                  >  2:2894661/1‑150 (MQ=255)
                   aaGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCctg                                                                                                                                 >  1:2097230/1‑150 (MQ=255)
                    aGAACAACATCAATCTGATCATTGTGGCTACCTCCATCGGCTTCGGCATGATCCCGATTGCAGCACCAGAATTCTACCATCAGTTCCCGGCCTGGTTTGAAACCATTTTTCATTCAGGCATCAGCTCCGCTGCATTGATGGCAATcgcgc                                                                                                                                <  1:2140708/150‑5 (MQ=255)
                    aGAACAACATCAATCTGATCATTGTGGCTACCTCCATCGGCTTCGGCATGATCCCGATTGCAGCACCAGAATTCTACCATCAGTTCCCGGCCTGGTTTGAAACCATTTTTCATTCAGGCATCAGCTCCGCTGCATTGATGGCAATcgcgc                                                                                                                                <  1:2140693/150‑5 (MQ=255)
                    aGAACAACATCAATCTGATCATTGTGGCTACCTCCATCGGCTTCGGCATGATCCCGATTGCAGCACCAGAATTATACCATCAGTTCCCGGCCTGGTTTGAAACCATTTTTCATTCAGGCATCAGCTCCGCTGCATTGATGGCAATcgcgc                                                                                                                                <  1:4743937/150‑5 (MQ=255)
                       acaacGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGa                                                                                                                             <  1:828671/150‑1 (MQ=255)
                          acGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAAcc                                                                                                                          <  2:1468110/150‑1 (MQ=255)
                          acGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAAcc                                                                                                                          <  2:3957203/150‑1 (MQ=255)
                              cAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGAt                                                                                                                      >  2:2317122/1‑150 (MQ=255)
                                   tGATCATTGTGGCTACCTCCATCGGCTTCGGCATGATCCCGATTGCAGCACCAGAATTCTACCATCAGTTCCCGGCCTGGTTTGAAACCATTTTTCATTCAGGCATCAGCTCCGCTGCATTGATGGCAATCGCGCTCAATCTCATGTTCa                                                                                                                 <  2:1800592/150‑1 (MQ=255)
                                   tGATCATTGTGGCTACCTCCATCGGCTTCGGCATGATCCCGATTGCAGCACCAGAATTCTACCATCAGTTCCCGGCCTGGTTTGAAACCATTTTTCATTCAGGCATCAGCTCCGCTGCATTGATGGCAATCGCGCTCAATCTCATGTTCa                                                                                                                 >  1:524326/1‑150 (MQ=255)
                                   tGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCa                                                                                                                 >  2:1386746/1‑150 (MQ=255)
                                   tGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGACCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCa                                                                                                                 <  2:722892/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                                                                                                         <  2:1997780/150‑1 (MQ=255)
                                           gtcgCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTc                                                                                                         >  1:2751369/1‑150 (MQ=255)
                                                ggCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTcaccac                                                                                                    <  2:1202518/150‑1 (MQ=255)
                                                    tCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGc                                                                                                >  2:3830751/1‑150 (MQ=255)
                                                        tGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACt                                                                                            >  1:1908925/1‑150 (MQ=255)
                                                         gggCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTc                                                                                           <  2:2106666/150‑1 (MQ=255)
                                                          ggCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCg                                                                                          <  2:1667635/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)
                                                                ggCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAAc                                                                                    >  1:4881736/1‑150 (MQ=255)
                                                                ggCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAAc                                                                                    <  2:129916/150‑1 (MQ=255)
                                                                 gCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAAcc                                                                                   >  1:2896391/1‑150 (MQ=255)
                                                                 gCATGATCCCGATTACCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAAcc                                                                                   >  2:4898481/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)
                                                                         ccGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAgtg                                                                           >  1:3638220/1‑150 (MQ=255)
                                                                         ccGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAgtg                                                                           >  2:1930112/1‑150 (MQ=255)
                                                                                  gCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCa                                                                  >  2:1841068/1‑150 (MQ=255)
                                                                                   cACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAg                                                                 >  2:155847/1‑150 (MQ=255)
                                                                                        aCCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTAc                                                            >  2:2194702/1‑150 (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)
                                                                                                     acTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATcc                                               >  1:2117387/1‑150 (MQ=255)
                                                                                                      cTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCa                                              <  2:1181783/150‑1 (MQ=255)
                                                                                                      cTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCa                                              >  1:440689/1‑150 (MQ=255)
                                                                                                      cTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCa                                              <  1:4709020/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)
                                                                                                          ccGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACcagtcagt                                                                            >  1:1946595/1‑116 (MQ=255)
                                                                                                                tGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCg                                                                       >  1:76648/1‑115 (MQ=255)
                                                                                                                tGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCg                                                                       <  2:76648/115‑1 (MQ=255)
                                                                                                                tGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGAc                                    >  1:3287275/1‑150 (MQ=255)
                                                                                                                tGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGAc                                    >  2:4512286/1‑150 (MQ=255)
                                                                                                                tGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGAc                                    >  1:3948989/1‑150 (MQ=255)
                                                                                                                tGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGAc                                    >  1:721411/1‑150 (MQ=255)
                                                                                                                   ttCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACAtt                                 >  1:3367172/1‑150 (MQ=255)
                                                                                                                   ttCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACAtt                                 >  1:3819062/1‑150 (MQ=255)
                                                                                                                    tCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACAttt                                >  2:4710028/1‑150 (MQ=255)
                                                                                                                      gagaCCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACATTTcc                              <  2:3173480/150‑1 (MQ=255)
                                                                                                                      gagaCCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACATTTcc                              >  1:1732209/1‑150 (MQ=255)
                                                                                                                      gagaCCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACATTTcc                              <  1:3171008/150‑1 (MQ=255)
                                                                                                                      gagaCCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTGCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTGTGCTGCAGCCTACGAGCGCACCATCCAGGATTCGGACATTTcc                              >  1:1740548/1‑150 (MQ=255)
                                                                                                                        gaCCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACATTTCcgc                            >  1:2963022/1‑150 (MQ=255)
                                                                                                                          ccATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACATTTCCgcgc                          <  2:1339670/150‑1 (MQ=255)
                                                                                                                           cATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCACCTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACATTTCCGCGCt                         >  2:3299761/1‑150 (MQ=255)
                                                                                                                           cATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACATTTCCGCGCt                         >  1:165773/1‑150 (MQ=255)
                                                                                                                           cATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACATTTCCGCGCt                         >  2:2589930/1‑150 (MQ=255)
                                                                                                                           cATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACATTTCCGCGCt                         >  2:82855/1‑150 (MQ=255)
                                                                                                                               ttCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACATTTCCGCGCTGCGt                     <  2:3110636/150‑1 (MQ=255)
                                                                                                                                 ccACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACATTTCCGCGCTGCGTGa                   >  2:3615515/1‑150 (MQ=255)
                                                                                                                                   aCTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACATTTCCGCGCTGCGTGATg                 >  1:613249/1‑150 (MQ=255)
                                                                                                                                     tCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACATTTCCGCGCTGCGTGATggg               <  2:2285501/150‑1 (MQ=255)
                                                                                                                                     tCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACATTTCCGCGCTGCGTGATggg               >  1:2829021/1‑150 (MQ=255)
                                                                                                                                     tCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACATTTCCGCGCTGCGTGATggg               <  2:2896391/150‑1 (MQ=255)
                                                                                                                                               gCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACATTTCCGCGCTGCGTGATGGGGACTACTTCa     <  2:3454581/150‑1 (MQ=255)
                                                                                                                                                cgcggcggCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACATTTCCGCGCTGCGTGATGGGGACTACTTCaa    <  2:840835/148‑1 (MQ=255)
                                                                                                                                                 tCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACATTTCCGCGCTGCGTGATGGGGACTACTTCAAg   <  1:2264448/150‑1 (MQ=255)
                                                                                                                                                  cggcggCGATCATGGCCATCCTGCTGATCCTGATCTTCAACCATTTCCCCACCGGCAACTCGGTAAACCTGTGGGTGTTCGCTGCAGCCTACGAGCTCACCCGCCAGTAGTCGGACATTTCCGCGGTGCGTGCGGGGGACTCGTTGAAgg  >  1:1278151/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACATTTCCGCGCTGCGTGATGGGGACTACTTCAAGG  >  NC_002947/5270506‑5270801

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: