Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,129,594 T→C 5.3% Q90R (CAG→CGG)  PP_5375 ← LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,129,5940TC 5.3% 478.4 / 9.7 152Q90R (CAG→CGG) PP_5375LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (66/76);  new base C (2/6);  total (68/84)
Fisher's exact test for biased strand distribution p-value = 2.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.63e-01

Only 100 of 153 total aligned reads displayed.

CCTGGCGTAATCCTTCACATGATCGTGGGCGATGATCCGCAGGTCGATGTTTTCATGGGCACTCTGGAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACTGGCGTATATCGCCGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCGGTGGGGGTCAGGGACATGCTACGGGTGGTACGGGTGAATAGTGCCTTGCCCAGGTAGTCCTCAAGAAGGCGAACCTGA  >  NC_002947/6129446‑6129743
                                                                                                                                                    |                                                                                                                                                      
 cTGGCGTAATCCTTCACATGATCGTGGGCGATGATCCGCAGGTCGATGTTTTCATGGGCACTCTGGAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACTgg                                                                                                                                                      <  1:727625/150‑1 (MQ=255)
  tGGCGTAATCCTTCACATGATCGTGGGCGATGATCCGCAGGTCGATGTTTTCATGGGCACTCTGGAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACTGGc                                                                                                                                                     <  1:1082091/150‑1 (MQ=255)
    gCGTAATCCTTCACATGATCGTGGGCGATGATCCGCAGGTCGATGTTTTCATGGGCACTCTGGAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACTGGCGt                                                                                                                                                   >  1:42569/1‑150 (MQ=255)
      gTAATCCTTCACATGATCGTGGGCGATGATCCGCAGGTCGATGTTTTCATGGGCACTCTGGAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACTGGCGtat                                                                                                                                                 >  1:537106/1‑150 (MQ=255)
       tAATCCTTCACATGATCGTGGGCGATGATCCGCAGGTCGATGTTTTCATGGGCACTCTGGAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACCGGCgatat                                                                                                                                                <  1:4280817/150‑5 (MQ=255)
         aTCCTTCACATGATCGTGGGCGATGATCCGCAGGTCGATGTTTTCATGGGCACTCTGGAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACTGGCGTATATc                                                                                                                                             <  1:1202886/150‑1 (MQ=255)
         aTCCTTCACATGATCGTGGGCGATGATCCGCAGGTCGATGTTTTCATGGGCACTCTGGAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACTGGCGTATATc                                                                                                                                             <  2:1389380/150‑1 (MQ=255)
          tCCTTCACATGATCGTGGGCGATGATCCGCAGGTCGATGTTTTCATGGGCACTCTGGAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACTGGCGTATATCg                                                                                                                                            >  1:4666689/1‑150 (MQ=255)
            cTTCACATGATCGTGGGCGATGATCCGCAGGTCGATGTTTTCATGGGCACTCTGGAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACTGGCGTATATCGcc                                                                                                                                          <  2:1092426/150‑1 (MQ=255)
             ttCACATGATCGTGGGCGATGATCCGCAGGTCGATGTTTTCATGGGCACTCTGGAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACTGGCGTATATCGCCg                                                                                                                                         >  2:1764851/1‑150 (MQ=255)
             ttCACATGATCGTGGGCGATGATCCGCAGGTCGATGTTTTCATGGGCACTCTGGAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACTGGCGTATATCGCCg                                                                                                                                         >  1:3639049/1‑150 (MQ=255)
             ttCACATGATCGTGGGCGATGATCCGCAGGTCGATGTTTTCATGGGCACTCTGGAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACTGGCGTATATCGCCg                                                                                                                                         >  1:2964876/1‑150 (MQ=255)
             ttCACATGATCGTGGGCGATGATCCGCAGGTCGATGTTTTCATGGGCACTCTGGAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGTGGCCTTGCCACTGGCGTATATCGCCg                                                                                                                                         >  1:407333/1‑150 (MQ=255)
               cacaTGATCGTGGGCGATGATCCGCAGGTCGATGTTTTCATGGGCACTCTGGAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACTGGCGTATATCGCCGGt                                                                                                                                       <  2:602999/150‑1 (MQ=255)
                 caTGATCGTGGGCGATGATCCGCAGGTCGATGTTTTCATGGGCACTCTGGAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACTGGCGTATATCGCCGGTgg                                                                                                                                     >  2:1553408/1‑150 (MQ=255)
                  aTGATCGTGGGCGATGATCCGCAGGTCGATGTTTTCATGGGCACTCTGGAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACTGGCGTATATCGCCGGTGgc                                                                                                                                    >  1:166669/1‑150 (MQ=255)
                  aTGATCGTGGGCGATGATCCGCAGGTCGATGTTTTCATGGGCACTCTGGAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACTGGCGTATATCGCCGGTGgc                                                                                                                                    >  2:2224911/1‑150 (MQ=255)
                  aTGATCGTGGGCGATGATCCGCAGGTCGATGTTTTCATGGGCACTCTGGAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACCGGCGATATACGCCAGTGgc                                                                                                                                    <  1:658566/150‑1 (MQ=255)
                  aTGATCGTGGGCGATGATCCGCAGGTCGATGTTTTCATGGGCACTCTGGAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACCGGCGATATACGCCAGTGgc                                                                                                                                    <  1:2100203/150‑1 (MQ=255)
                     aTCGTGGGCGATGATCCGCAGGTCGATGTTTTCATGGGCACTCTGGAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACTGGCGTATATCGCCGGTGGcgcg                                                                                                                                 >  1:2882589/1‑150 (MQ=255)
                      tCGTGGGCGATGATCCGCAGGTCGATGTTTTCATGGGCACTCTGGAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACTGGCGTATATCGCCGGTGGCGCgg                                                                                                                                <  2:3217845/150‑1 (MQ=255)
                             cGATGATCCGCAGGTCGATGTTTTCATGGGCACTCTGGAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACTGGCGTATATCGCCGGTGGCGCGGGCCAGGt                                                                                                                         <  1:2132480/150‑1 (MQ=255)
                                 gatCCGCAGGTCGATGTTTTCATGGGCACTCTGGAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACTGGCGTATATCGCCGGTGGCGCGGGCCAGGTGAAg                                                                                                                     >  1:4800342/1‑150 (MQ=255)
                                  atCCGCAGGTCGATGTTTTCATGGGCACTCTGGAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACTGGCGTATATCGCCGGTGGCGCGGGCCAGGTGAAgc                                                                                                                    <  1:4679509/150‑1 (MQ=255)
                                   tCCGCAGGTCGATGTTTTCATGGGCACTCTGGAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACTGGCGTATATCGCCGGTGGCGCGGGCCAGGTGAAgca                                                                                                                   >  2:1408727/1‑150 (MQ=255)
                                    ccGCAGGTCGATGTTTTCATGGGCACTCTGGAATTCCGACACCTTGGGTAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACTGGCGTATATCGCCGGTGGCGCGGGCCAGGTGAAgcat                                                                                                                  <  1:3281834/150‑1 (MQ=255)
                                     cGCAGGTCGATGTTTTCATGGGCACTCTGGAAGTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACTGGCGTATATCGCCGGTGGCGCGGGCCAGGTGAAgcatg                                                                                                                 <  2:2329462/150‑1 (MQ=255)
                                      gCAGGTCGATGTTTTCATGGGCACTCTGGAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACTGGCGTATATCGCCGGTGGCGCGGGCCAGGTGAAgcatgc                                                                                                                <  1:4442938/150‑1 (MQ=255)
                                          gTCGATGTTTTCATGGGCACTCTGGAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACTGGCGTATATCGCCGGTGGCGCGGGCCAGGTGAAGCATGCATTc                                                                                                            <  2:2533608/150‑1 (MQ=255)
                                           tCGATGTTTTCATGGGCACTCTGGAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACTGGCGTATATCGCCGGTGGCGCGGGCCAGGTGAAGCATGCATTCa                                                                                                           <  2:37707/150‑1 (MQ=255)
                                            cgttGTTTTCATGGGCACTCTGGAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGTTTGCCTTGCCACTGGCGTATATCGCCGGTGGCGCGGGCCAGGTGAAGCATGCATTCAc                                                                                                          <  1:3984524/147‑1 (MQ=255)
                                              aTGTTTTCATGGGCACTCTGGAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACTGGCGTATATCGCCGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGc                                                                                                        >  2:1232355/1‑150 (MQ=255)
                                                 ttttCATGGGCACTCTGGAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACTGGCGTATATCGCCGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCAcc                                                                                                     >  2:4684371/1‑150 (MQ=255)
                                                     cATGGGCACTCTGGAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACTGGCGTATATCGCCGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGttt                                                                                                 >  2:293405/1‑150 (MQ=255)
                                                     cATGGGCACTCTGGAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACTGGCGTATATCGCCGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGttt                                                                                                 >  1:1785269/1‑150 (MQ=255)
                                                     cATGGGCACTCTGGAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACTGGCGTATATCGCCGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGttt                                                                                                 >  2:4037494/1‑150 (MQ=255)
                                                      aTGGGCACTCTGGAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACTGGCGTATATCGCCGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTc                                                                                                >  1:2632853/1‑150 (MQ=255)
                                                       tGGGCACTCTGGAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACTGGCGTATATCGCCGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCg                                                                                               <  2:4595806/150‑1 (MQ=255)
                                                        gggCACTCTGGAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACTGGCGTATATCGCCGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGt                                                                                              >  1:2084512/1‑150 (MQ=255)
                                                          gCACTCTGGACTTCCGACACCGTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGGTGCCTTGCCACTGGCGTATAGCGCCGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCggac                                                                                            <  2:3673776/150‑4 (MQ=255)
                                                          gCACTCTGGAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACTGGCGTATATCGCCGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTaa                                                                                            <  2:38424/150‑1 (MQ=255)
                                                              tctGGAATTCCGACACCTTGGGCAGTAACCACGGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACTGGCGTATATCGCCGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACt                                                                                        >  1:654753/1‑150 (MQ=255)
                                                              tctGGAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACTGGCGTATATCGCCGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACt                                                                                        >  1:1895305/1‑150 (MQ=255)
                                                              tctGGAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACTGGCGTATATCGCCGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACt                                                                                        <  2:1879990/150‑1 (MQ=255)
                                                              tctGGAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACTGGCGTATATCGCCGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACt                                                                                        >  1:1235501/1‑150 (MQ=255)
                                                                tGGAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACTGGCGTATATCGCCGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTgc                                                                                      >  2:389685/1‑150 (MQ=255)
                                                                tGGAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACTGGCGTATATCGCCGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTgc                                                                                      >  1:3902215/1‑150 (MQ=255)
                                                                 ggAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACTGGCGTATATCGCCGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTgcg                                                                                     >  1:2670002/1‑150 (MQ=255)
                                                                  gAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACTGGCGTATATCGCCGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTgcgc                                                                                    <  1:2505615/150‑1 (MQ=255)
                                                                   aaTTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACTGGCGTATATCGCCGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTgcgcg                                                                                   >  1:684551/1‑150 (MQ=255)
                                                                        cGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACTGGCGTATATCGCCTGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCGGt                                                                              >  1:479183/1‑150 (MQ=255)
                                                                        cGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACTGGCGTATATCGCCGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCGGt                                                                              >  1:3745222/1‑150 (MQ=255)
                                                                        cGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACTGGCGTATATCGCCGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCGGt                                                                              >  1:3745331/1‑150 (MQ=255)
                                                                         gACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACTGGCGTATATCGCCGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCGGTg                                                                             >  1:3230291/1‑150 (MQ=255)
                                                                          acacCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACTGGCGTATATCGCCGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCGGTgg                                                                            >  1:152455/1‑150 (MQ=255)
                                                                            acCTTGGGCAGTAACCACAGCGACGCCATGGCGGTAGTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACTGGCGTATATCGCCGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCGGTgggg                                                                          <  2:4422839/150‑1 (MQ=255)
                                                                            acCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACTGGCGTATATCGCCGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCGGTgggg                                                                          >  1:453155/1‑150 (MQ=255)
                                                                            acCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTGGTTGGTTGCCTTGCCACTGGCGTATATCGCCGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCGGTgggg                                                                          <  2:4611042/150‑1 (MQ=255)
                                                                            acCGTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACTGGCGTATATCGCCGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCGGTgggg                                                                          <  2:2035898/150‑1 (MQ=255)
                                                                              cTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACTGGCGTATATCGCCGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCGGTGGGGGt                                                                        <  2:4226138/150‑1 (MQ=255)
                                                                              cTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACTGGCGTATATCGCCGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCGGTGGGGGt                                                                        >  1:3374386/1‑150 (MQ=255)
                                                                                 gggCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACTGGCGTATATCGCCGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCGGTGGGGGTCAg                                                                     <  2:1208941/150‑1 (MQ=255)
                                                                                     aGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACTGGCGTATATCGCCGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCGGTGGGGGTCAGGGAc                                                                 <  1:4523384/150‑1 (MQ=255)
                                                                                     aGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACTGGCGTATATCGCCGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCGGTGGGGGTCAGGGAc                                                                 >  2:4724407/1‑150 (MQ=255)
                                                                                       tAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACTGGCGTATATCGCCGGTGGCGCGGGCCAGGTGACGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCGGTGGGGGTCAGGGACAt                                                               <  2:2776991/150‑1 (MQ=255)
                                                                                       tAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACTGGCGTATATCGCCGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCGGTGGGGGTCAGGGACAt                                                               <  2:3038361/150‑1 (MQ=255)
                                                                                        aaCCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACTGGCGTATATCGCCGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCGGTGGGGGTCAGGGACATg                                                              <  2:3445680/150‑1 (MQ=255)
                                                                                         aCCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACTGGCTTATATCGCCGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCTCCGGTGGGGGTCAGGGACATGc                                                             <  2:1885098/150‑1 (MQ=255)
                                                                                         aCCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACTGGCGTATATCGCCGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCGGTGGGGGTCAGGGACATGc                                                             >  1:2362830/1‑150 (MQ=255)
                                                                                          ccACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACTGGCGTATATCGCCGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCGGTGGGGGTCAGGGACATGCt                                                            <  2:668540/150‑1 (MQ=255)
                                                                                           cacaGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACTGGCGTATATCGCCGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCGGTGGGGGTCAGGGACATGCTa                                                           >  1:3866427/1‑150 (MQ=255)
                                                                                           cacaGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACTGGCGTATATCGCCGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCGGTGGGGGTCAGGGACATGCTa                                                           >  2:2269261/1‑150 (MQ=255)
                                                                                            acaGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACTGGCGTATATCGCCGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCGGTGGGGGTCAGGGACATGCTAc                                                          >  1:2162616/1‑150 (MQ=255)
                                                                                              aGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACTGGCGTATATCGCCGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCGGTGGGGGTCAGGGACATGCTACgg                                                        <  1:1495923/150‑1 (MQ=255)
                                                                                               gCGACGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACTGGCGTATATCGCCGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCGGTGGGGGTCAGGGACATGCTACggg                                                       >  1:3023779/1‑150 (MQ=255)
                                                                                                  aCGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACTGGCGTATATCGCCGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCGGTGGGGGTCAGGGACATGCTACGggtgg                                                    >  2:3986355/1‑150 (MQ=255)
                                                                                                    gCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACTGGCGTATATCGCCGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCAGTGGGGGTCAGGGACATGCTACGGGTGGTa                                                  <  1:2463489/150‑1 (MQ=255)
                                                                                                        tGGCGGTACTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACTGGCGTATATCGCCGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCGGTGGGGGTCAGGGACATGCTACGGGTGGTACggg                                              <  1:4887523/150‑1 (MQ=255)
                                                                                                         ggCGGTACGGGTGGCAACGGGGACGTGGTGGTTGCCTTGCCACGGGCGTAGATCGCCGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCGGGGGGGGTCAGGGACATGCTACGGGTGGTACGGGt                                             <  2:2603895/150‑1 (MQ=255)
                                                                                                                  ggtggCAACGGTGACTTGTTGGTTGCCTTGCCACTGGCGTATATCGCCGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCGGTGGGGGTCAGGGACATGCTACGGGTGGTACGGGTGAATAGTGc                                    <  1:2944910/150‑1 (MQ=255)
                                                                                                                   gtggCAACGGTGACTTGTTGGTTGCCTTGCCACTGGCGTATATCGCCGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCGGTGGGGGTCAGGGACATGCTACGGGTGGTACGGGTGAATAGTGcc                                   <  2:2882589/150‑1 (MQ=255)
                                                                                                                      gCAACGGTGACTTGTTGGTTGCCTTGCCACTGGCGTATATCGCCGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCGGTGGGGGTCAGGGGCATGCTACGGGTGGTACGGGTGAATAGTGCCTTg                                >  2:4438818/1‑150 (MQ=255)
                                                                                                                      gCAACGGTGACTTGTTGGTTGCCTTGCCACTGGCGTATATCGCCGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCGGTGGGGGTCAGGGACATGCTACGGGTGGTACGGGTGAATAGTGCCTTg                                <  2:648635/150‑1 (MQ=255)
                                                                                                                            gTGACTTGTTGGTTGCCTTGCCACTGGCGTATATCGCCGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCGGTGGGGGTCAGGGACATGCTACGGGTGGTACGGGTGAATAGTGCCTTGCCCAgg                          <  2:1589230/150‑1 (MQ=255)
                                                                                                                            gTGACTTGTTGGTTGCCTTGCCACGGGCGTATATCGCCGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTGTCGTAATACTGCGCGCCGGTGGGGGTCAGGGACATGCTACGGGTGGTACGGGTGAATAGTGCCTTGCCCAgg                          <  2:3674288/150‑1 (MQ=255)
                                                                                                                              gACTTGTTGGTTGCCTTGCCACTGGCGTATATCGCCGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCGGTGGGGGTCAGGGACATGCTACGGGTGGTACGGGTGAATAGTGCCTTGCCCAGGTa                        <  2:42569/150‑1 (MQ=255)
                                                                                                                                    ttggttgCCTTGCCACTGGCGTATATCGCCGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCTCCGTTTCGTAATACTGCGCGCCGGTGGGGGTCAGGGACATGCTACGGGTGGTACGGGTGAATAGTGCCTTGCCCAGGTAGTCCTc                  <  1:1716629/150‑1 (MQ=255)
                                                                                                                                     tggttgCCTTGCCACTGGCGTATATCGCCGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCGGTGGGGGTCAGGGACATGCTACGGGTGGTACGGGTGAATAGTGCCTTGCCCAGGTAGTCCTCa                 <  2:4666689/150‑1 (MQ=255)
                                                                                                                                        ttgccttgccACTGGCGTATATCGCCGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCGGTGGGGGTCAGGGACATGCTACGGGTGGTACGGGTGAATAGTGCCTTGCCCAGGTAGTCCTCaaga              >  2:2427444/1‑150 (MQ=255)
                                                                                                                                         gcccgcgccACCGGCGATATACGCCAGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCGGTGGGGGTCAGGGACATGCTACGGGTGGTACGGGTGAATAGTGCCTTGCCCAGGTAGTCCTCaagaa             >  2:3830378/7‑150 (MQ=255)
                                                                                                                                            cgcgccACCGGCGATATACGCCAGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCGGTGGGGGTCAGGGACATGCTACGGGTGGTACGGGTGAATAGTGCCTTGCCCAGGTAGTCCTCAAGAAGGc          >  1:1006274/4‑150 (MQ=255)
                                                                                                                                              tgccACTGGCGTATATCGCCGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCGGTGGGGGTCAGGGACATGCTACGGGTGGTACGGGTGAATAGTGCCTTGCCCAGGTAGTCCTCAAGAAGGCGa        <  2:859362/150‑1 (MQ=255)
                                                                                                                                               gccACTGGCGTTTATCGCCGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCGGTGGGGGTCAGGGACATGCTACGGGTGGTACGGGTGTATAGTGCCTTGCCCAGGTAGTCGTCATGAAGGCGaa       <  2:407333/150‑1 (MQ=255)
                                                                                                                                               gccACTGGCGTATATCGCCGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCGGTGGGGGTCAGGGACATGCTACGGGTGGTACGGGTGAATAGTGCCTTGCCCAGGTAGTCCTCAAGAAGGCGaa       <  2:2964876/150‑1 (MQ=255)
                                                                                                                                               gccACTGGCGTATATCGCCGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCGGTGGGGGTCAGGGACATGCTACGGGTGGTACGGGTGAATAGTGCCTTGCCCAGGTAGTCCTCAAGAAGGCGaa       <  2:407345/150‑1 (MQ=255)
                                                                                                                                               gccACTGGCGTATATCGCCGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCGGTGGGGGTCAGGGACATGCTACGGGTGGTACGGGTGAATAGTGCCTTGCCCAGGTAGTCCTCAAGAAGGCGaa       <  1:4449844/150‑1 (MQ=255)
                                                                                                                                                 cACTGGCGTATATCGCCGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCGGTGGGGGTCAGGGACATGCTACGGGTGGTACGGGTGAATAGTGCCTTGCCCAGGTAGTCCTCAAGAAGGCGAAcc     >  2:4711812/1‑150 (MQ=255)
                                                                                                                                                 cACTGGCGTATATCGCCGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCGGTGGGGGTCAGGGACATGCTACGGGTGGTACGGGTGAATAGTGCCTTGCCCAGGTAGTCCTCAAGAAGGCGAAcc     >  1:4094558/1‑150 (MQ=255)
                                                                                                                                                  aCTGGCGTATATCGCCGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCGGTGGGGGTCAGGGACATGCTACGGGTGGTACGGGTGAATAGTGCCTTGCCCAGGTAGTCCTCAAGAAGGCGAACCt    >  2:4082672/1‑150 (MQ=255)
                                                                                                                                                    tGGCGTATATCGCCGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCGGTGGGGGTCAGGGACATGCTACGGGTGGTACGGGTGAATAGTGCCTTGCCCAGGTAGTCCTCAAGAAGGCGAACCTGa  <  2:4800342/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
CCTGGCGTAATCCTTCACATGATCGTGGGCGATGATCCGCAGGTCGATGTTTTCATGGGCACTCTGGAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTTGTTGGTTGCCTTGCCACTGGCGTATATCGCCGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCGGTGGGGGTCAGGGACATGCTACGGGTGGTACGGGTGAATAGTGCCTTGCCCAGGTAGTCCTCAAGAAGGCGAACCTGA  >  NC_002947/6129446‑6129743

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: